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ABSTRACT 



Almost none of the r- modes ordinarily found in rotating stars exist, if the star 
and its perturbations obey the same one-parameter equation of state; and rotating 
relativistic stars with one-parameter equations of state have no pure r-modes at 
all, no modes whose limit, for a star with zero angular velocity, is an axial-parity 
oscillation. Rotating stars of this kind similarly have no pure g-modes, no modes 
whose spherical limit is a perturbation with polar parity and vanishing perturbed 
pressure and density. Where have these modes gone? 

In spherical stars of this kind, r-modes and g-modes form a degenerate zero- 
frequency subspace. We find that rotation splits the degeneracy to zeroth order in 
the star's angular velocity Q, and the resulting modes are generically hybrids, whose 
limit as Q — > is a stationary current with axial and polar parts. Lindblom and 
Ipser have recently found these hybrid modes in an analytic study of the Maclaurin 
spheroids. We present the first calculation of these modes in relativistic stars. 
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Because each mode has definite parity, its axial and polar parts have alternating 
values of I. We show that each mode belongs to one of two classes, axial- led or 
polar-led, depending on whether the spherical harmonic with lowest value of I that 
contributes to its velocity field is axial or polar. We numerically compute these modes 
for slowly rotating newtonian polytropes and Maclaurin spheroids, and for slowly 
rotating relativistic stars with uniform density. Timescales for the gravitational-wave 
driven instability and for viscous damping are computed for the hybrid modes of the 
newtonian models using assumptions appropriate to neutron stars. The instability 
to nonaxisymmetric modes is, as expected, dominated by the / = m r-modes with 
simplest radial dependence, the only modes which retain their axial character in 
newtonian isentropic models. For relativistic isentropic stars, these I = m modes 
are replaced for I > 2 by axial-led hybrids. We find analytically the post-newtonian 
corrections to these modes for uniform density stars. 
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Chapter 1 



Introduction 



1.1 Background and Motivation 

This dissertation examines a new class of oscillation modes of rotating stars. The 
work has been motivated by the recently discovered r-mode instability (see below), 
and answers a number of previously unresolved questions concerning the nature of 
the r-mode spectrum in newtonian and relativistic stellar models. 

The structure and stability of rotating relativistic stars has recently been reviewed 
in detail (Stergioulas [[75], Friedman |22|, |23|], Friedman and Ipser f24|), and a general 



discussion of the small oscillations of relativistic stars may be found in a recent review 
article by Kokkotas |4(| (see also Kokkotas and Schmidt In this work, we will 

focus our attention on non-radial oscillations, which were first studied in relativistic 
stars by Thorne and collaborators (Thorne and Campolattaro [f79H , Price and Thorne 



3j, Thorne |76, T7\, Campolattaro and Thorne fill , Ipser and Thorne |3 



The spherical symmetry of a non-rotating star implies that its perturbations can 
be divided into two classes, polar or axial, according to their behaviour under parity. 
Where polar tensor fields on a 2-sphere can be constructed from the scalars Y™ and 
their gradients VY™ (and the metric on a 2-sphere), axial fields involve the pseudo- 
vector f x VV;" 1 , and their behavior under parity is opposite to that of Y™. That is, 
axial perturbations of odd / are invariant under parity, and axial perturbations with 
even / change sign. 

It is useful to further divide stellar perturbations into subclasses according to 
the physics dominating their behaviour. This classification was first developed by 
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Cowling [jnj for the polar perturbations of newtonian polytropic models. The f- and 
p-modes are polar-parity modes having pressure as their dominant restoring force. 
They typically have large pressure and density perturbations and high frequencies 
(higher than a few kilohertz for neutron stars). The other class of polar-parity modes 
are the g-modes, which are chiefly restored by gravity. They typically have very small 
pressure and density perturbations and low frequencies. Indeed, for isentropic stars, 
which are marginally stable to convection, the g-modes are all zero-frequency and 
have vanishing perturbed pressure and density (see Sect. 2.1). Similarly, all axial- 
parity perturbations of newtonian perfect fluid models have zero frequency in a non- 
rotating star. The perturbed pressure and density as well as the radial component 
of the fluid velocity are all rotational scalars and must have polar parity. Thus, the 
axial perturbations of a spherical star are simply stationary horizontal fluid currents 
(see Sect. 2.1). 

The analogues of these modes in relativistic models of neutron stars have been 
studied by many authors. More recently, an additional class of outgoing modes has 
been identified that exist only in relativistic stars. Like the modes of black holes, 
these are essentially associated with the dynamical spacetime geometry and have 
been termed w-modes, or gravitational wave modes. Their existence was first argued 
by Kokkotas and Schutz [fl2"| |. The polar w-modes were first found by Kojima 37 



as 



rapidly damped modes of weakly relativistic models, while the axial w-modes were 



first studied by Chandrasekhar and Ferrari |14[] as scattering resonances of highly 
relativistic models. (See the reviews by Kokkotas 0] and Kokkotas and Schmidt 

In general, this classification of modes also describes the oscillations of rotating 
stars, although the character of the modes may be significantly affected by rotation. 
Because a rotating star is also invariant under parity, its perturbations can be classi- 
fied according to their behaviour under parity. If a mode varies continuously along a 
sequence of equilibrium configurations that starts with a spherical star and continues 
along a path of increasing rotation, the mode will be called axial if it is axial for 
the spherical star. Its parity cannot change along the sequence, but I is well-defined 
only for modes of the spherical configuration. 

Rotation imparts a finite frequency to the axial-parity perturbations of newtonian 
models. Because these modes are restored by the Coriolis force, their frequencies are 
proportional to the star's angular velocity, Q. These rotationally restored axial modes 



3 



were first studied by Papaloizou and Pringle [p0[ , who called them r-modes because 
of their similarity to the Rossby waves of terrestrial meteorology. For a normal mode 
of the form e t ( <jt + m v) ^ Papaloizou and Pringle found the r-mode frequency to be, 

„ 2mVL 

It is only rather recently that the oscillation modes of rotating relativistic stars 
have begun to be accessible to numerical study (see below). Early work on the 
perturbations of such stars focused mainly on the criteria for their stability, and 
led to the surprising discovery that all rotating perfect fluid stars are subject to a 
non-axisymmetric instability driven by gravitational radiation. The instability was 



discovered by Chandrasekhar [|12| in the I = m = 2 polar mode of the uniform- 
density, uniformly rotating Maclaurin spheroids. Although this mode is unstable 
only for rapidly rotating models, by looking at the canonical energy of arbitrary 
initial data sets, Friedman and Schutz [^SJ and Friedman [2T| showed the instability 
to be a generic feature of rotating perfect fluid stars. 

In essence, the CFS (Chandrasekhar-Friedman-Schutz) instability operates by 
converting the rotational energy of the star partly into the oscillation energy of the 
perturbation and partly into gravitational waves. For a normal mode of the form 
e i(at+mip) nonaxisymmetric instability acts in the following manner. 

In a non-rotating star, gravitational radiation removes positive angular momen- 
tum from a forward moving mode and negative angular momentum from a backward 
moving mode, thereby damping all time- dependent, non-axisymmetric modes. In a 
star rotating sufficiently fast, however, a backward moving mode can be dragged 
forward as seen by an inertial observer; and it will then radiate positive angular 
momentum. The mode continues to carry negative angular momentum because the 
perturbed star has lower total angular momentum than the unperturbed star. As 
positive angular momentum is removed from a mode with negative angular momen- 
tum, the angular momentum of the mode becomes increasingly negative, implying 
that the amplitude of the mode increases. Thus, the mode is driven by gravitational 
radiation. 

Since the instability acts on modes that are retrograde with respect to the star, 
but prograde as seen by an inertial observer, a mode will be unstable if and only if 
its frequency satisfies the condition, 

(x(> + mft)<0. (2) 
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For the polar f- and p-modes, the frequency is large and approximately real. Condi- 
tion (fj) will be met only if \mQ\ is of order |er|, so that for a given angular velocity 
the instability will set in first through modes with large m. 

The CFS instability spins a star down by allowing it to radiate away its angular 
momentum in gravitational waves. However, to determine whether this mechanism 
may be responsible for limiting the rotation rates of actual neutron stars, one must 
also consider the effects of viscous damping on the perturbations. Detweiler and 



Lindblom j[J| suggested that viscosity would stabilize any mode whose growth time 
was longer than the viscous damping time, and this was confirmed by Lindblom 
and Hiscock [fi~9|| . Recent work has indicated that the gravitational- wave- driven 
instability can only limit the rotation rate of hot neutron stars, with temperatures 
above the superfluid transition point, T ~ 10 9 K, but below the temperature at 
which bulk viscosity apparently damps all modes, T ~ 10 10 K. (Ipser and Lindblom 
[33|; Lindblom ||7] and Lindblom and Mendell [|51[]) Because of uncertainties in 
the temperature of the superfluid phase transition and in our understanding of the 
dominant mechanisms for effective viscosity, even this brief temperature window is 
not guaranteed. 

To calculate the timescales associated with viscous and radiative dissipation it is 
necessary to compute explicitly the normal modes of oscillation. Until recently, the 
polar f- and p-modes were expected to dominate the CFS instability through their 
coupling to mass multipole radiation. As we have already noted, all axial-parity fluid 
oscillations are time-independent in a spherical model, and therefore do not couple 
to gravitational radiation at all. (Thorne and Campolattaro J79[) In a rotating star, 
the rotationally restored r-modes do couple to current multipole radiation. However, 
their low frequencies and negligible perturbed densities in newtonian stars made it 
seem implausible that their contribution to gravitational radiation would compare 
to that of the polar-parity modes. 

Indeed, apart from studies of the axial-parity oscillations of models of the neutron 
star crust (van Horn |SIJ, Schumaker and Thorne []TT|), the axial modes were almost 
universally ignored in the early research on perturbations of relativistic stars. It has 
only been recently that interest in these modes has been revived, following the work 



of Chandrasekhar and Ferrari on the resonant scattering of axial wave modes [14 



and on the coupling between axial and polar modes induced by stellar rotation [13 



(Further recent studies of axial modes are reviewed by Kokkotas JIO] and Kokkotas 
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and Schmidt [0.) 

Thus, the first explicit calculations of the dissipative timescales associated with 
neutron star oscillations focused on the / = m f-modes (Ipser and Lindblom |33 



Lindblom | 47fl , Lindblom and Mendell [[51]]), and until very recently it was only these 
modes that had been studied in connection with the CFS instability. It has long 
been hoped that some neutron stars rotate sufficiently fast to be subject to the CFS 
instability, and that the gravitational radiation produced might be detectable by 
gravitational wave observatories. However, based on the initial studies of dissipa- 
tion in f-mode oscillations, the prospects were rather unpromising. Neutron stars 
formed from stellar collapse would certainly be hot enough to pass through the tem- 
perature window at which viscous damping is apparently suppressed, but there was 
little evidence that neutron stars formed in supernovae (or by the accretion-induced 
collapse of white dwarves) rotate rapidly enough for the onset of instability. Follow- 



ing an early suggestion of Papaloizou and Pringle |5I |, Wagoner |8T] had proposed 
another scenario in which an old, accreting neutron star, spun up past the onset of 
nonaxisymmetric instability would achieve an equilibrium state with angular momen- 
tum acquired by accretion balanced by angular momentum radiated in gravitational 
waves. This scenario, too, appeared to have been ruled out by the strength of damp- 
ing by mutual friction and viscosity at the temperatures expected for such stars, 
T ~ 10 8 K. 

Very recently, however, a series of surprising results have emerged that dramati- 
cally improve these prospects. 

The first surprise was the discovery that the r-modes are CFS unstable in perfect 
fluid models with arbitrarily slow rotation. First indicated in numerical work by An- 
dersson B, the instability is implied in a nearly newtonian context by the newtonian 
expression for the r-mode frequency ([!]), which satisfies the CFS instability criterion, 
(H), for arbitrarily small Q, 

2(1 - l)(l + 2)m 2 Q 2 

A computation by Friedman and Morsink [^| of the canonical energy of initial data 
showed (independent of assumptions on the existence of discrete modes) that the 
instability is a generic feature of axial-parity fluid perturbations of relativistic stars. 

As we have just observed, the generic instability of perfect fluid models will be of 
no astrophysical importance if, in actual stars, the unstable modes are damped by 



6 



viscous dissipation. Studies of the viscous and radiative timescales associated with 
the r-modes (Lindblom et al. pi], Owen et al. p9| , Andersson et al. ||, Kokkotas 
and Stergioulas fl3| , Lindblom et al. [f)3[0 have revealed a second surprising result: 
The growth time of r-modes driven by current-multipole gravitational radiation is 
significantly shorter than had been expected. In fact, it has turned out be so short for 
some of the r-modes that their instability to gravitational radiation reaction easily 
dominates viscous damping in hot, newly formed neutron stars. A neutron star that 
is rapidly rotating at birth now appears likely to spin down by radiating most of 
its angular momentum in gravitational waves. (See, however, the caveats indicated 
below.) 

Hot on the heels of these theoretical surprises was the discovery by Marshall et.al. 
[]57[ of a fast (16ms) pulsar in a supernova remnant (N157B) in the Large Magellanic 
Cloud. Estimates of the initial period put it in the 6-9ms range, thus providing the 
long-sought evidence of a class of neutron stars that are formed rotating rapidly. 
Hence, the newly discovered instability appears to set the upper limit on the spin of 
the newly discovered class of neutron stars! 

The current picture that has emerged of the spin-down of a hot, newly formed 
neutron star can be readily understood in terms of a model of the r-mode instability 
due to Owen, Lindblom, Cutler, Schutz, Vecchio and Andersson (hereafter OLCSVA) 
||59|| . Since one particular mode (with spherical harmonic indices I = m = 2 and 
frequency a = — 4Q/3) is expected to dominate the r-mode instability, the perturbed 
star is treated as a simple system with two degrees of freedom: the uniform angular 
velocity Q of the equilibrium star, and the (dimensionless) amplitude a of the I = 
m = 2 r-mode. Initially, the neutron star forms with a temperature large enough for 
bulk viscosity to damp any unstable modes, T ^ 10 10 K; the star is assumed to be 



rotating close to its maximum (Kepler) velocity, VLk ~ y M/-R 3 , the angular velocity 
at which a particle orbits the star's equator. The star then cools by neutrino emission 
at a rate given by a standard power law cooling formula (Shapiro and Teukolsky 
|[72|). Once it reaches the temperature window at which the I = m = 2 r-mode can 
go unstable, the system is assumed to evolve in three stages. 

First, the amplitude of the r-mode undergoes rapid exponential growth from some 
arbitrary tiny magnitude. Using conservation of energy and angular momentum, 
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OLCSVA derive the following equations for the evolution of the system in this stage. 

dn 2VL a 2 Q 



dt Ty 1 + a 2 Q 
da a a 1 — a 2 Q 



(4) 
(5) 



dt \t gr \ t v 1 + a 2 Q 
Here, tqr and Ty are, respectively, the timescales for the growth of the mode by 
gravitational radiation reaction and the damping of the mode by viscosity (see Sect. 
2.6). (The parameter Q is a constant of order 0.1 related to the initial angular mo- 
mentum and moment of inertia of the equilibrium star.) Since the initial amplitude 
a of the mode is so small, the angular momentum changes very little at first (Eq. 
(01)). That this stage is characterized by the rapid exponential growth of a is the 
statement that the first term in Eq. (H) (the radiation reaction term) dominates over 
the second (viscous damping). 

Eventually the mode will grow to a size at which linear perturbation theory is 
insufficient to describe its behaviour. It is expected that a non-linear saturation will 
occur, halting the growth of the mode at some amplitude of order unity, although the 
details of these non-linear effects are poorly understood at present. When this satu- 
ration occurs, the system enters a second evolutionary stage during which the mode 
amplitude remains essentially unchanged and the angular momentum of the star is 
radiated away. During this stage OLCSVA evolve their model system according to 
the equations 

a 2 = k (6) 
dVt 2VL kQ 
~dt ~ ~\t gr \ 1-kQ 

where k is constant of order unity parameterizing the uncertainty in the degree of 
non-linear saturation. The star spins down by Eq. (|7|), radiating away most of its 
angular momentum while continuing to cool gradually. 

When its temperature and angular velocity are low enough that viscosity again 
dominates the gravitational- wave- driven instability, the mode will be damped. Dur- 
ing this third stage, OLCSVA return to Eqs. (f§)-([5|) to continue the evolution of 
their system. That the mode amplitude decays is the statement that the second 
term in Eq. (|) (the viscous damping term) dominates the first (radiation reaction), 
at this temperature and angular velocity. 



s 



The net effect of this three-stage evolutionary process is that the newly formed 
neutron star is left with an angular velocity small compared with Qk- This final 
angular velocity appears to be fairly insensitive to the initial amplitude of the mode 
and to its degree of non-linear saturation. A final period P ^ 5 — 10ms appar- 
ently rules out accretion-induced collapse of white dwarves as a mechanism for the 
formation of millisecond pulsars with P ^ 3ms. 

The r-mode instability has also revived interest in the Wagoner |5T| mechanism, 
involving old neutron stars spun up by accretion to the point at which the accretion 
torque is balanced by the angular momentum loss in gravitational radiation. Bildsten 
|| and Andersson, Kokkotas and Stergioulas have proposed that the r-mode 
instability might succeed in this regard where the instability to polar modes seems 
to fail. However, the mechanism appears to be highly sensitive to the temperature 
dependence of viscous damping. Levin |4!| has argued that if the r-mode damping 
is a decreasing function of temperature (at the temperatures expected for accreting 
neutron stars, T ~ 10 8 K) then viscous reheating of the unstable neutron star could 
drive the system away from the Wagoner equilibrium state. Instead, the star would 
follow a cyclic evolution pattern. Initially, the runaway reheating would drive the 
star further into the r-mode instability regime and spin it down to a fraction of 
its angular velocity. Once it has slowed to the point at which the r-modes become 
damped, it would again slowly cool and begin to spin up by accretion. Eventually, 
it would again reach the critical angular velocity for the onset of instability and 
repeat the cycle. Since the radiation spin-down time is of order 1 year, while the 
accretion spin-up time is of order 10 6 years, the star spends only a small fraction 
of the cycle emitting gravitational waves via the unstable r-modes. This would 
significantly reduce the likelihood that detectable gravitational radiation is produced 
by such sources. On the other hand, if the r-mode damping is independent of - or 
increases with - temperature (at T ~ 10 8 K) then the Wagoner equilibrium state may 
be allowed (Levin [45]). Work is currently in progress (Lindblom and Mendell p2| ) 
to investigate the r-mode damping by mutual friction in superfluid neutron stars, 
which was the dominant viscous mechanism responsible for ruling out the Wagoner 
scenario in the first place (Lindblom and Mendell |5lfl ). 

Other uncertainties in the scenarios described above are still to be investigated. 
There is substantial uncertainty in the cooling rate of neutron stars, with rapid 
cooling expected if stars have a quark interior or core, or a kaon or pion condensate. 
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Madsen |57J suggests that an observation of a young neutron star with a rotation 
period below 5 — 10ms would be evidence for a quark interior; but even without 
rapid cooling, the uncertainty in the superfTuid transition temperature may allow a 
superffuid to form at about 10 10 K, possibly killing the instability. We noted above 
the expectation that the growth of the unstable r-modes will saturate at an amplitude 
of order unity due to non-linear effects (such as mode- mode couplings); however, this 
limiting amplitude is not yet known with any certainty and could be much smaller. 
In particular, it has been suggested that the non-linear evolution of the r-modes will 
wind up the magnetic field of a neutron star, draining energy away from the mode 
and eventually suppressing the unstable modes entirely (Rezzolla, Lamb and Shapiro 



p6| ; see also Spruit [(75|). 

The excitement over the r-mode instability has generated a large literature. (An- 



dersson 0, Friedman and Morsink [25], Kojima [3£], Lindblom et al. 54], Owen et 
al. |59fl , Andersson, Kokkotas and Schutz ||, Kokkotas and Stergioulas P3f , Ander- 
sson, Kokkotas and Stergioulas ||, Madsen p6| , Hiscock [[H[], Lindblom and Ipser 
[§g, Bildsten [§, Levin [g§, Ferrari et al. [gj, Spruit 0, Brady and Creighton §, 
Lockitch and Friedman |55|] (see Ch. 2), Lindblom et al. [Q, Beyer and Kokkotas 
0, Kojima and Hosonuma [|3£|, Lindblom [48], Schneider et al. JFU], Rezzolla et al. 
]|67f , Yoshida and Lee fl8~3| ) It has also generated a number of questions which have 
not been properly answered, some of which are addressed in this dissertation. 

Despite the sudden interest in the r-modes they are not yet well-understood for 
stellar models appropriate to neutron stars. A neutron star is accurately described 
by a perfect fluid model in which both the equilibrium and perturbed configurations 
obey the same one-parameter equation of state. Hereafter, I will call such models 
isentropic, because isentropic models and their adiabatic perturbations obey the 
same one-parameter equation of state. 

For stars with more general equations of state, the r-modes appear to be complete 
for perturbations that have axial-parity. However, this is not the case for isentropic 
models. Early work on the r-modes focused on newtonian models with general 
equations of state (Papalouizou and Pringle JB0| , Provost et al. p4J, Saio p9" 



Smeyers and Martens |73|) and mentioned only in passing the isentropic case. In 
isentropic newtonian stars, one finds that the only purely axial modes allowed are 
the r-modes with I = m and simplest radial behavior. (Provost et al. f64}|p|; see Sect. 
*An appendix in this paper incorrectly claims that no I = m r-modes exist, based on an incorrect 
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2.3.1) It is these r-modes only that have been studied (and found to be physically 
interesting) in connection with the gravitational-wave driven instability. 

The first part of this dissertation (Ch. 2) addresses the question of the missing 
modes in isentropic newtonian models. (Lockitch and Friedman |]55|| ) The disap- 
pearance of the purely axial modes with I > m occurs for the following reason. 
We have already noted that all axial perturbations of a spherical star are time- 
independent convective currents with vanishing perturbed pressure and density. We 
have also noted that in spherical isentropic stars the gravitational restoring forces 
that give rise to the g-modes vanish and they, too, become time-independent con- 
vective currents with vanishing perturbed pressure and density. Thus, the space of 
zero frequency modes, which generally consists only of the axial r-modes, expands 
for spherical isentropic stars to include the polar g-modes. This large degenerate 
subspace of zero-frequency modes is split by rotation to zeroth order in the star's 
angular velocity, and the corresponding modes of rotating isentropic stars are gener- 
ically hybrids whose spherical limits are mixtures of axial and polar perturbations. 
These hybrid modes have already been found analytically for the uniform- density 
Maclaurin spheroids by Lindblom and Ipser [|5(| in a complementary presentation 
that makes certain features transparent but masks properties (such as their hybrid 
character) that are our primary concern. Lindblom and Ipser point out that since 
these modes are also restored by the Coriolis force, it is natural to refer to them as 
rotation modes, or generalized r-modes. 

Having found the missing modes in isentropic newtonian stars, I then turn to the 
corresponding problem in general relativity. The r-modes of rotating relativistic stars 
have been studied for the first time only recently (Andersson [[|; Kojima |5S|| ; Beyer 
and Kokkotas J7|; Kojima and Hosonuma ||39||), but none of these calculations have 
found the modes in the isentropic stellar models appropriate to neutron stars. As 
in the newtonian spherical isentropic relativistic star has a large degenerate 

subspace of zero-frequency modes consisting of the axial-parity r-modes and the 
polar-parity g-modes. Again, the degeneracy is split by rotation and the generic 
mode of a rotating isentropic star is a hybrid whose spherical limit is a mixture 
of axial and polar perturbations. The second part of this dissertation (Chs. 3- 
4) presents the first calculation finding these modes in isentropic relativistic stars, 
assumption about their radial behavior. 
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Although isentropic newtonian stars retain a vestigial set of purely axial modes 
(those having I = m), rotating relativistic stars of this type have no pure r-modes, 
no modes whose limit for a spherical star is purely axial. Instead, the newtonian 
r-modes with / = m > 2 acquire relativistic corrections with both axial and polar 
parity to become discrete hybrid modes of the corresponding relativistic models (see 
Sects. 3.4.1-3.4.2). 

This dissertation examines the hybrid rotational modes of rotating isentropic 
stars, both newtonian and relativistic. Sect. 1.2 begins with a brief summary of the 
theory of self-gravitating perfect fluids and their linearized perturbations. Ch. 2 con- 
siders the hybrid modes of newtonian stars, first proving that the time-independent 
modes of spherical isentropic stars are the r- and g-modes (Sect. 2.1), and then mov- 
ing to consider rotating stars (Sect. 2.2). Sect 2.3 distinguishes two types of modes, 
axial-led and polar-led, and shows that every mode belongs to one of the two classes. 
Sects. 2.4-2.5 deal with the computation of eigenfunctions and eigenfrequencies for 
modes in each class, adopting what appears to be a method that is both novel and 
robust. For the uniform-density Maclaurin spheroids, these modes have been found 
analytically by Lindblom and Ipser. I find machine precision agreement with their 
eigenfrequencies and corresponding eigenfunctions to lowest nontrivial order in the 
angular velocity Q. I also examine the frequencies and modes of n — 1 polytropes, 
finding that the structure of the modes and their frequencies are very similar for the 
polytropes and the uniform- density configurations. The numerical analysis is com- 
plicated by a curious linear dependence in the Euler equations, detailed in Appendix 
B. The linear dependence appears in a power series expansion of the equations about 
the origin. It may be related to difficulty other groups have encountered in searching 
for these modes. Finally, Sect. 2.6 examines unstable modes, computing their growth 
time and expected viscous damping time. The pure I = m = 2 r-mode retains its 
dominant role, but the 3 < I = m ^ 10 r-modes and some of the fastest growing 
hybrids remain unstable in the presence of viscosity. 

Chs. 3 and 4 are concerned with the hybrid modes of isentropic relativistic stars, 
which turn out be very similar in character to their newtonian counterparts. In Sect. 
3.1, the proof that the time-independent modes of spherical isentropic stars are the 
r- and g-modes is generalized to relativity. In Sect. 3.2 the perturbation equations 
governing the hybrid modes in slowly rotating stars are derived; their structure 
parallels the corresponding newtonian equations of Sect. 2.2. This similarity between 
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the newtonian and relativistic equations leads to an identical structure of the mode 
spectrum and to a parallel theorem in Sect 3.3 that every non-radial mode is either 
an axial-led or polar-led hybrid (the result has so far been proven only for slowly 
rotating relativistic stars). This chapter concludes with a discussion of the boundary 
conditions appropriate to the hybrid modes and the construction of some explicit 
solutions (Sect 3.4). I show that there are no modes in isentropic relativistic stars 
whose limit as Q —>■ is a pure axial perturbation with 1^1. In particular, the 
newtonian r-modes having I = m > 2 do not exist in isentropic relativistic stars and 
must be replaced by axial-led hybrid modes (Sect. 3.4.1). I explicitly construct these 
particular modes to first post-newtonian order in slowly rotating, uniform density 
stars (Sect. 3.4.2). 

Finally, Ch. 4 involves the computation of eigenfunctions and eigenfrequencies, 
applying essentially the same numerical method as was used in the newtonian calcu- 
lation. A set of modes from each parity class is constructed for uniform density stars 
and compared with their newtonian counterparts. The relativistic corrections turn 
out to be small for the modes and stellar models considered. As in the newtonian 
calculation, the numerical analysis is complicated by a curious linear dependence in 
the perturbation equations. The linear dependence, again, appears in a power series 
expansion of the equations about the origin, and is discussed in Appendix D. 

Throughout this dissertation I will work in geometrized units, (G = c = 1), 
except in Sect. 2.6 where G and c are restored to their cgs values for the explicit 
computation of dissipative timescales. I use the conventions of Misner, Thorne and 
Wheeler |5^] for the metric signature (—,+,+,-1-) and the sign of the curvature 
tensors. I adopt the abstract index notation (see, e.g., Wald, Sect 2.4 fl82| ) with 
latin spatial indices and greek spacetime indices; components of tensors will always 
be written with respect to a choice of coordinates (t, r, 9, ip). 

1.2 Self- Gravitating Perfect Fluids and Their Lin- 
earized Perturbations 

We will be considering stationary perfect fluid stellar models in both newtonian grav- 
ity and in general relativity. We construct both rotating and non-rotating equilibrium 
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stars and then study the equations of motion, linearized about these equilibria, gov- 
erning their small oscillations. We will make use of both the eulerian and lagrangian 
perturbation formalisms, which we now briefly review. 

1.2.1 Newtonian Gravity 

In the newtonian theory of gravity, a complete description of an isentropic perfect 
fluid configuration is provided by the fluid density p, the pressure p, the fluid velocity 
v a and the newtonian gravitational potential $. These must satisfy a barotropic 
(one-parameter) equation of state, 

P = p(j>), ( 8 ) 

the equation of mass conservation, 

d t p + V a (pv a ) = 0, (9) 

Euler's equation, 

(d t + v b V b )v a + -\7 a p + V a $ = 0, (10) 
P 

and the newtonian gravitational equation, 

V 2 $ = 4np. (11) 

An equilibrium stellar model is a time-independent solution (p,p, v a , <&) to these 
equations. The small perturbations of such a star may be studied using either the 
eulerian or the lagrangian perturbation formalism (see, e.g., Friedman and Schutz 



pY|). If [p(A),p(A), v a (X), $(A)] is a smooth family of solutions to the exact equations 
(|8[)-(|TT|) that coincides with the equilibrium solution at A = 0, 

[p(0),p(0),t; a (0),<l(0)] = (p,p,<<l>), 

then the eulerian change 5Q in a quantity Q may be defined (to linear order in A) 
as, 

(12) 



ah - dQ 



A=0 

Thus, an eulerian perturbation is simply a change (Sp, Sp, 5v a , 5$) in the equilibrium 
configuration at a particular point in space. 
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In the lagrangian perturbation formalism (Friedman and Schutz |[27||), on the 



other hand, perturbed quantities are described in terms of a lagrangian displacement 
vector £ a that connects fluid elements in the equilibrium and perturbed star. The 
lagrangian change AQ in a quantity Q is related to its eulerian change 5Q by 

AQ = 5Q + £^Q, (13) 

with £^ the Lie derivative along £ a . The fluid perturbation is then entirely deter- 
mined by the displacement £ a : 

Av a = dtC (14) 

^ = * = _ v „r, (i5) 

IP P 

(where 7 is the adiabatic index) and the corresponding Eulerian changes are 

5v a = (d t + £ v )C (16) 

sp = _ 7P v a e a -e°v p (is) 

with the change in the gravitational potential determined by 

V 2 5$ = Ait5p. (19) 



1.2.2 General Relativity 

In general relativity, a complete description of an isentropic perfect fluid configura- 
tion is provided by a spacetime with metric g a p, sourced by an energy- momentum 
tensor, 

T a p = (e + p)u a u p + pg aP , (20) 
where the fluid 4-velocity u a is a unit timelike vector field, 

u a u a = -1, (21) 

and e and p are, respectively, the total energy density and pressure of the fluid as 
measured by an observer moving with 4-velocity u a . The metric and fluid variables 
must, again, satisfy a barotropic equation of state, 



p =p(e), 



(22) 



as well as the Einstein field equations, 



G a/ 3 = 8nT Q/3 . (23) 

An equilibrium stellar model is a stationary solution (g Q/ 3,w a , e,p) to these equa- 
tions. The small perturbations of such star may be studied using either the eulerian 
or the lagrangian perturbation formalism (Friedman |^TJ; see also Friedman and Ipser 
|p4|). As in the newtonian case, an eulerian perturbation may described in terms 
of a smooth family, [g a /3(X), u a (\), e(A),p(A)], of solutions to the exact equations 
(]2~T|)-(|2"3"D that coincides with the equilibrium solution at A = 0, 

[M0),M a (0),e(0),p(0)] = (g a ^u a ,e,p). 



Then the eulerian change 5Q in a quantity Q may be defined (to linear order in A) 

as, 

rid 

(24) 



5Q =lx 



A=0 

Thus, an eulerian perturbation is simply a change (Sg a/ 3, 5u a , 5e, 5p) in the equilib- 
rium configuration at a particular point in spacetime. 

In the lagrangian perturbation formalism (Friedman pi]; see also Friedman and 
Ipser on the other hand, perturbed quantities are expressed in terms of the 

eulerian change in the metric h a p = 5g a/ 3, and a lagrangian displacement vector £ a , 
which connects fluid elements in the equilibrium star to the corresponding elements 
in the perturbed star. The lagrangian change AQ in a quantity Q is related to its 
eulerian change SQ by 

AQ = 6Q + £ ( Q, (25) 

with £^ the Lie derivative along £ Q . 
The identities, 

Ag al3 = h a/3 + 2V(a€p) (26) 
1 

Ae a/3l s = 2 £a P7 s 9^Ag^ u (27) 
then allow one to express the fluid perturbation in terms of h a $ and £ Q , 

Au a = -uWAfl^ (28) 
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where n is the baryon density and q a/3 = g a ^ + u a u@. Using Eqs. (25)-( f29|) , it is 
straightforward to express the corresponding eulerian changes also in terms of h a p 
and £ a . 
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Chapter 2 
Newtonian Stars 



The oscillation modes considered in this dissertation are dominantly restored by the 
Coriolis force and have frequencies that scale with the star's angular velocity, Q. 
Thus, all of these modes will be degenerate at zero frequency in a non-rotating star. 
To study properly the spectrum of these rotationally restored modes, we must first 
examine the perturbations of a non-rotating star, and find all of the modes belonging 
to its degenerate zero-frequency subspace. 



2.1 Stationary Perturbations of Spherical Stars 

We consider a static spherically symmetric, self-gravitating perfect fluid described 
by a gravitational potential $, density p and pressure p. These satisfy an equation 
of state of the form 

p = p(p), (30) 
as well as the newtonian equilibrium equations 

-V aP + V a $ = (31) 
P 

V 2 $ = 47rp. (32) 

We are interested in the space of zero-frequency modes, the linearized time- 
independent perturbations of this static equilibrium. This zero-frequency subspace 
is spanned by two types of perturbations: (i) perturbations with Sv a 7^ and bp = 
5p = 5$ = 0, and (ii) perturbations with Sp, 5p and <5<3> nonzero and Sv a = 0. If one 
assumes that no solution to the linearized equations governing a static equilibrium is 
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spurious, that each corresponds to a family of exact solutions, then the only solutions 
(ii) are spherically symmetric, joining neighboring equilibria. 

The decomposition into classes (i) and (ii) can be seen as follows. The set of 
equations satisfied by (Sp, Sp, Sv a ) are the perturbed mass conservation equation, 



S[d t p + V a (pv a )] = 0, 



(33) 



the perturbed Euler equation 
5 



(d t + £ v )v a + jVaP + V a ($ - \v 2 ' 



0. 



(34) 



the perturbed Poisson equation, <5[Eq. and an equation of state for the per- 

turbed configuration (which may, in general, differ from that of the equilibrium 
configuration). 

For a time-independent perturbation these equations take the form 



and 



V a (pSv a ) = 0, 

- V a 5p - ^V aP + V a 5$ = 0, 
p p 2 



V 2 5$ = AixSp. 



(35) 
(36) 

(37) 



Because Eq. ( |3"oD for Sv a decouples from Eqs. ( PB[ ) and ( |3"7| ) for (Sp, Sp, £$), any 
solution to Eqs. (|35| ) -(|37 D is a superposition of a solution (0, 0, 0, Sv a ) and a solution 
(Sp, Sp, 5$, 0). This is the claimed decomposition. 

The theorem that any static self-gravitating perfect fluid is spherical implies that 
the solution (Sp, Sp, <5$, 0) is spherically symmetric, to within the assumption that the 
static perturbation equations have no spurious solutions ( "linearization stability" )[]. 

Thus, under the assumption of linearization stability we have shown that all 
stationary non-radial (/ > 0) perturbations of a spherical star have Sp = Sp = 5$ = 
and a velocity field Sv a that satisfies Eq. (|3"o]). 

A perturbation with axial parity has the form (see, e.g., Friedman and Morsink 



H) 



Sv a = f/(r)e afec V 6 F z m V c r. 



(3f 



*We are aware of a proof of this linearization stability for relativistic stars under assumptions 
on the equation of state that would not allow polytropes (Kiinzle and Savage H|). 



and automatically satisfies Eq. (j35|). 

A perturbation with polar parity perturbation has the form, 



5v a = ^iy, m V a r + V{r)V a Y, m ; (39) 
r 

and Eq. (|35| ) gives a relation between W and V, 

4-(rpW)-l(l + l)pV = 0. (40) 
dr 

These perturbations must satisfy the boundary conditions of regularity at the center, 
r = and surface, r = R, of the star. Also, the lagrangian change in the pressure 
(defined in the next section) must vanish at the surface of the star. These boundary 
conditions result in the requirement that 

W(0) = W(R) = 0; (41) 

however, apart from this restriction, the radial functions U(r) and W(r) are unde- 
termined. 

Finally, we consider the equation of state of the perturbed star. For an adiabatic 
oscillation of a barotropic star (i.e., a star that satisfies a one-parameter equation of 
state, p = p(p)) Eq. (|l"5"l) implies that the perturbed pressure and energy density are 
related by 

' p__ p_ 

p 7p_ 

for some adiabatic index j(r) which need not be the function 



5p = 6p + 
IP P 



(42) 



TV ^ - P d P fAO\ 

!» = --r (43) 
p dp 

associated with the equilibrium equation of state. Here, £ a is the lagrangian dis- 
placement vector and is related to our perturbation variables by Eq. (|16D , which 
becomes 

Sv a = d t e, (44) 
or (taking the initial displacement (at t = 0) to be zero) 

C =tSv r . (45) 

For the class of perturbations under consideration, we have seen that 5p = 5p = 0, 
thus Eqs. (f42|) and ( |45l) require that 
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For axial-parity perturbations this equation is automatically satisfied, since Sv a has 
no r-component (Eq. fl38|)). Thus, a spherical barotropic star always admits a class 
of zero-frequency r-modes. 

For polar-parity perturbations, 5v r oc W(r) ^ 0, and Eq. (|46| ) will be satisfied if 
and only if 

l(r) = r(r) = P - d J. (47) 
p dp 

Thus, a spherical barotropic star admits a class of zero-frequency g-modes if and 
only if the perturbed star obeys the same one-parameter equation of state as the 
equilibrium star. We call such a star isentropic, because isentropic models and their 
adiabatic perturbations obey the same one-parameter equation of state. 

Summarizing our results, we have shown the following. A spherical barotropic 
star always admits a class of zero-frequency r-modes (stationary fluid currents with 
axial parity); but admits zero-frequency g-modes (stationary fluid currents with polar 
parity) if and only if the star is isentropic. Conversely, the zero-frequency subspace 
of non-radial perturbations of a spherical isentropic star is spanned by the r- and 
g-modes - that is, by convective fluid motions having both axial and polar parity 
and with vanishing perturbed pressure and density.0 Being stationary, these modes 
do not couple to gravitational radiation. One would expect this large subspace of 
modes, which is degenerate at zero-frequency, to be split by rotation, so let us now 
consider the perturbations of rotating stars. 



2.2 Perturbations of Rotating Stars 

We consider perturbations of an isentropic newtonian star, rotating with uniform 
angular velocity Q. No assumption of slow rotation will be made until we turn to 
numerical computations in Sect. 2.4. The equilibrium of an axisymmetric, self- 
gravitating perfect fluid is described by the gravitational potential $, density p, 
pressure p and a 3-velocity 

v a = tt(p a , (48) 

where (p a is the rotational Killing vector field. 

2 Note that for spherical stars, nonlinear couplings invalidate the linear approximation after a 
time t ~ R/Sv, comparable to the time for a fluid element to move once around the star. For 
nonzero angular velocity, the linear approximation is expected to be valid for all times, if the 
amplitude is sufficiently small, roughly, if \5v\ < R£l. 
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We will use the lagrangian perturbation formalism reviewed in Sect. 1.2.1. Since 
the equilibrium spacetime is stationary and axisymmetric, we may decompose our 
perturbations into modes of the formP] e^ <Tt+mip \ In this case, the Eulerian change in 
the 3-velocity (|I"6| ) is related to the lagrangian displacement £ a by, 



5v a = i(a + mQ)C- (49) 

We can expand this perturbed fluid velocity in vector spherical harmonics 

00 r l i 
5v ° = E -WiY^r + V^Y™ - iUie^VbY^Vcr e lu \ (50) 

l=m ^ r J 

and examine the perturbed Euler equation. 

The lagrangian perturbation of Euler's equation is 

= A[(d t + £ v )v a + V a (h-~v 2 + ®)] 

= (d t + £ v )Av a + V a [A(h-^v 2 + <!>)], (51) 

and its curl, which expresses the conservation of circulation for an isentropic star, is 

= q a = i( a + mn)e abc V b Av c , (52) 

or 

= q a = i(a + mn)e abc V b 5v c + tte abc V h (£ 6v y c ) . (53) 

Using the spherical harmonic expansion (EDI) of 5v a we can write the components 
of q a as 



1 00 f 

= q r = — \ [(<?+ mQ)l(l + 1) - 2mn}U l Y l m (54) 



l=m 



-2ft [sin ^YT™ + 1(1 + 1) cos AY]" 1 ] 
+2VLWi [sin 6d e Y^ + 2 cos OY^} \ e iat 



3 We will always choose to > since the complex conjugate of an to < mode with frequency a 
is an m > mode with frequency —a. Note that a is the frequency in an inertial frame. 
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= q b 



r 2 sin 



1 00 f / w \ 



(55) 



-2Qd r Vi cos 9 sm 9d e Y m + 2fim 2 — K m 

r 



-2Qd r Wi sin 2 ^F™ - 2mVtd r Ui cos 0^™ 



+ (a + mn)d r Ui sin 0d e K m + 2mfi— sin 9d Y m \e iat , 

r 



and 

= q ip 



- 2 sin 2 9 



\ m{a + mn)d r UiY l m - 2fi9 r E7i cos 9 sin 03f^ m 

+2tt— [m 2 -1(1 + 1) sin 2 9}Y m - 2mQd r V l cos 9Y T 
r 



(56) 



+ 



/ Wi\ Vi 

[a + mO) (d r Vi + 2mQ— 

V r / r 



x sin^Y,™ f e i<T *. 



These components are not independent. The identity V a q a = 0, which follows from 
equation (152), serves as a check on the right-hand sides of (plp-(|5l)D. 
Let us rewrite these equations making use of the standard identities, 



sinWT" = ZQj+i^-tf + 1)3^1 
cos 9Y t m = Qi+iYftx + QiY™ x 



(57) 
(58) 



where 



Qi 



(I + m)(l — m) 



_(2/-l)(2/ + l) 
Defining a dimensionless comoving frequency 



K 



(er + mQ) 

n 



(59) 



(60) 



we find that the equation q r = becomes 

oo ( 

= E [±kI(1 + 1) - mpiYf 



l=m 



+(W l -lVi)(l + 2)Q l+1 Y% 1 



-[Wt + il + mjil-VQtY^A, 



q 9 = becomes 
= 



l=m 



E { -Q1+1Q1+2 - w7 



J Z+2 



(m — -kI)U[ — ml— 



Y 



i+i 



+ 



Q«m + (Z + 1)Q 2 - ZQ? +1 ) VJ' 



-(l-^-g? +1 )Wi'-i«m^ + m 2 ^ 



yr. 



m+^(Z + l))[/; + m(Z + l)^ 



(i + i)v; + w{ 



r Z-2 



and g v = becomes 



o = E -iQi+iQm c//-(z + i)^ 

Z=m L 



I l+2 



1 . Vj 1 W/' 
-nl — m)V, + ml nl — 

2 11 r 2 r . 



+ 



(^m + (z + i)Qf-zgf +1 )c/; 



+ (m 2 -/(/ + !) (l - Q? - g? +1 )) ^ 
Qk(Z + 1) + m) ^' + m(Z + 1)^ - i/t(Z + 1) 



r . 
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+ (l + i)Qi-iQi 



uui u -i 



1 l-2 



where ' = 4-. 

ar 

Let us rewrite the equations one last time using the orthogonality relation for 
spherical harmonics, 

J Y l 7 l 'Yrdn = 5 ll 'S mm/ , (64) 

where dfl is the usual solid angle element. 

From equation (|61|) we find that J q r Y l * m dfl = gives 

= [±kI(1 + l)-m]U l + (l + - (Z - - IQi+i [W l+1 + (/ + 2)^ +1 ] (65) 

Similarly, / q 9 Y l * m dQ = gives 

= Q l Q l - 1 {(l-2)V{_ 2 -W;_ 2 } + Q l \[m-±K(l-l)}UU-m(l-l) — 

12 r 

i-Qf- Qf+i) w/ - [l^rn + + i)q 2 - zq 2 +1 ] v;' + ^rr3 - m 2 ^ 



+ Q m |[m + ^(/ + 2)}U' l+1 + m(l + 2)^} 

- Q, +2 Q m {(z + 3)v/ +2 + w' } 



(66) 



and / q*YC m dVL = gives 



-(Z-2)g,Q,_ 1 

1 



J-2 



(J-l) 



u, 



1-2 



+ 



Km+(l + 1)01 - l Qi+i 



U! 



(I + 3)Q l+2 Q 



i+i 



U' l+2 + (l + 2) 



U l+2 



m 2 -1(1 + 1) (1-Q 2 -Q 



'l+l 



El 

r 



+ Qi {[^(Z - 1) - rn]VU + m(l-l 
- Qi+i {["«(* + 2) + m}V( +l + m(l + 2 



K(l — 1) 

r 2 v ; r 

_ i K(i + 2) !*» 

r 2 r 



(67) 



2.3 The Character of the Perturbation Modes 

From this last form of the equations it is clear that the rotation of the star mixes 
the axial and polar contributions to 5v a . That is, rotation mixes those terms in ( |50|) 
whose limit as Q — > is axial with those terms in ((5^) whose limit as Q — > is 
polar. It is also evident that the axial contributions to 5v a with I even mix only with 
the odd / polar contributions, and that the axial contributions with I odd mix only 
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with the even / polar contributions. In addition, we prove that for non-axisymmetric 
modes the lowest value of I that appears in the expansion of Sv a is always I = m. 
(When to = this lowest value of I is either or 1.) 

For an equilibrium model that is axisymmetric and invariant under parity, one 
can resolve any degeneracy in the perturbation spectrum to make each discrete mode 
an eigenstate of parity with angular dependence e tm(fi . The following theorem then 
holds. 

Theorem 1 Let (Sp, Sv a ) with Sv a /Oka discrete normal mode of a uniformly 
rotating stellar model obeying a one-parameter equation of state. Then the decom- 
position of the mode into spherical harmonics Y™ (i.e., into (l,m) representations 
of the rotation group about the center of mass) has I = m as the lowest contributing 
value of I, when m ^ 0; and has or 1 as the lowest contributing value of I, when 
m = 0. 

Thus, we find two distinct classes of mixed, or hybrid, modes with definite be- 
havior under parity. Let us call a non-axisymmetric mode an "axial-led hybrid" (or 
simply "axial-hybrid") if Sv a receives contributions only from 

axial terms with / = to, to + 2, to + 4, ... and 
polar terms with / = to + 1, to + 3, to + 5, .... 

Such a mode has parity (— l) m+1 . 

Similarly, we define a non-axisymmetric mode to be a "polar-led hybrid" (or 
"polar-hybrid" ) if Sv a receives contributions only from 

polar terms with / = to, to + 2, to + 4, ... and 
axial terms with / = to + 1, to + 3, m + 5, .... 

Such a mode has parity (— l) m . 

For the case to = 0, there exists a set of axisymmetric modes with parity +1 that 
we call "axial-led hybrids" since Sv a receives contributions only from 

axial terms with / = 1, 3, 5, . . . and 
polar terms with / = 2, 4, 6, . . .; 

and there exist two sets of axisymmetric modes that may be designated as "polar-led 
hybrids." One set has parity —1 and Sv a receives contributions only from 
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polar terms with I = 1, 3, 5, . . . and 
axial terms with / = 2, 4, 6, . . .. 

The other set has parity +1 and 5v a receives contributions only from 

polar terms with I = 0, 2, 4, . . . and 
axial terms with I = 1, 3, 5, . . .. 

Note that the theorem holds for p-modes as well as for the rotational modes 
that are our main concern. A p-mode is determined by its density perturbation and 
is therefore dominantly polar in character regardless of its parity. For a rotational 
mode, however, the lowest / term in its velocity perturbation is at least comparable 
in magnitude to the other contributing terms. 

We prove the theorem separately for each parity class in Appendix A. 



2.3.1 The purely axial solutions 



We have proved that the generic mode of a rotating isentropic newtonian star is 
a hybrid mixture of axial and polar terms. However, it is known that newtonian 
stars of this type do allow a set of purely axial modes (Provost et al. |63]]). To find 
these r-modes, let us assume that the only non- vanishing coefficient in the spherical 
harmonic expansion (|50"D of the perturbed velocity is Ui(r), for some particular value 
of I. Eqs. (|65D-(|67D must be satisfied for all I, but with our ansatz they reduce to 
the following set. 
Eq. QBjD becomes, 

1 Ui = 0, (68) 



+ 1) — m 

and Eq. (0) with I -»• I + 2, I -»• I and I -»• I- 



-IQ1+1Q1+2 



u;-(i + i) 



2 gives the equations 
Ui 



(69) 



- K m+(l + l)Qf 



u; 



(70) 



+ 



m 



1(1 + 1) (1-Qf-QLi 



= (Z + l)Q,_iQ, 



Ul + l 



u. 



(71) 
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respectively. Recall that we need only work with two of the three equations (|65|)-(|67D, 
since they are linearly dependent as a result of the identity, V a q a = 0. 

We see immediately from Eq. (^) that a non-trivial solution to these equations 
exists if and only if 

K ^( a + mfi) = ___ (72) 



which is the r-mode frequency, Eq. (p]), found by Papalouizou and Pringle |KJ 



By Thm. ([]]), we know that a non-axisymmetric (m > 0) mode must have I = m 
as its lowest value of / and that an axisymmetric (m = 0) mode must have I = 1 
as its lowest value of I. Hence, in the present context of pure a spherical harmonic 
these are also the only allowed values of /. An r-mode with I > m (or I > 1 if 
m = 0) cannot exist in isentropic newtonian models. We consider the axisymmetric 
and non-axisymmetric cases separately. 

The case m = and 1 = 1. 

It is well known that uniform rotation is a purely axial perturbation with m = and 
I = 1, and we can see this from our Eqs. (|5IS|)-([72]) as follows. 

With m = and / = 1, Eq. ( |72|) simply becomes kQ = a = 0. The radial 
behaviour of this stationary solution is then determined by the other equations. The 
definition of Qi, Eq. ([59D, gives, 



Q?_i = 0, Ql= l - and Q 2 l+ i = ^- (73) 

These imply that Eq. (|n]) is trivially satisfied, while Eqs. ( |69|) and ( |70|) both 
become, 

= Ul- -U { (74) 
r 

or 

Ui(r) = Kr\ (75) 
for some constant K. If we define the constant 

ci=-i[iy K , (76) 

then our perturbed 3-velocity Q50D simply becomes^, 

Sv a = -iKr 2 e abc V h Y?V c r = (l<p a , (77) 



4 We use the standard normalization for the spherical harmonic Y® = y/S/Air cos 8. (See, e.g. 
Jackson p§, p.99.) 
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which represents a small change Cl in the uniform angular velocity Q of the star, as 
claimed. This perturbed velocity field is displayed in Fig. (|1|). 

The case I — m > 0. 



Eqs. (|69D-(|T2"D also have a simple solution when I = m > 0. The Papalouizou and 
Pringle frequency (J72[) becomes, 

Ktt = (a + mn) = - -. (78) 

v ; (m + 1) v ; 

The definition of Qi, Eq. (p9f), gives, 



Ql = 0, and gk 1 = ^-L_ (79) 

These, again, imply that Eq. (|7T|) is trivially satisfied, while Eqs. fl69|) and (|70|) both 
become, 

= EC - ^^^m (80) 



j, \ m+l 



with solution, 

EV» = • (81) 

We have chosen the normalization of this solution so that U m = 1 at the surface 
of the star, r = R. Images of the perturbed velocity field 5v a for the r-modes with 
1 < I = m < 3 are shown in Fig. ([I]). We note that the mode with I = m = 1 has 
zero-frequency in the inertial frame, 

o- = (re - 1)0 = 0, (82) 

and represents rotation of the star about an axis perpendicular to its original axis of 
rotation. In addition, we note that the r-mode with I = m = 2 is the one expected to 
dominate the gravitational radiation driven instability in hot, young neutron stars. 



2.4 Numerical Method 



In our numerical solution, we restrict consideration to slowly rotating stars, finding 
axial- and polar-led hybrids to lowest order in the angular velocity Q. That is, we 
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assume that perturbed quantities introduced above obey the following ordering in 
Q: 

W,~0(1), V5~0(l), Ut-Oil), 
5p~0(tt), 5p~0(tt), 5$~0(fi), <T~0(fi). 
The Q — > limit of such a perturbation is a sum of the zero-frequency axial and 
polar perturbations considered in Sect. 2.2. Note that, although the relative orders 
of 5p and 5v a are physically meaningful, there is an arbitrariness in their absolute 
order. If (Sp, Sv a ) is a solution to the linearized equations, so is (Q5p,Q5v a ). We 
have chosen the order fl53|) to reflect the existence of well-defined, nontrivial velocity 
perturbations of the spherical model. Other authors (e.g., Lindblom and Ipser H50f ) 
adopt a convention in which 5v a = 0(Q) and 5p = 0(Q 2 ). 

To lowest order, the equations governing these perturbations are the perturbed 
Euler equations fl65|)-(6?D and the perturbed mass conservation equation, (pop, which 
becomes 

rW; + (l+ r^jWi- 1(1 + 1)^ = 0. (84) 

In addition, the perturbations must satisfy the boundary conditions of regularity 
at the center of the star, r = 0, regularity at the surface of the star, r = R, and the 
vanishing of the lagrangian change in the pressure at the surface of the star, 

= Ap = dp + = CP' +0(Q). (85) 

Equations ( |S^) -(|67D and (|4j) are a system of ordinary differential equations for 
Wi'(r), Vi'(r) and Ui>(r) (for all I'). Together with the boundary conditions, these 
equations form a non-linear eigenvalue problem for the parameter k, where kQ is the 
mode frequency in the rotating frame. 

To solve for the eigenvalues we proceed as follows. We first ensure that the bound- 
ary conditions are automatically satisfied by expanding Wp{r), Vj/(r) and £fy(r) (for 
all V) in regular power series about the surface and center of the star. Substituting 
these series into the differential equations results in a set of algebraic equations for 
the expansion coefficients. These algebraic equations may be solved for arbitrary 
values of k using standard matrix inversion methods. For arbitrary values of k, how- 
ever, the series solutions about the center of the star will not agree with those about 
the surface of the star. The requirement that the series agree at some matching 
point, < r < R, then becomes the condition that restricts the possible values of 
the eigenvalue, kq. 



30 



The equilibrium solution (p, $) appears in the perturbation equations only through 
the quantity (p'/p) in equation (|84|). We begin by writing the series expansion for 
this quantity about r = as 

{pi Rf,"\R 



' M 4f>f;i, <«>> 



i odd 



and about r = i? as 
where the 7Tj and 7Tfc are determined from the equilibrium solution. 



f> ' ] ' E , (87) 



Because ( j65|) relates i7j(r) algebraically to Wi±i(r) and VJ±i(r), we may eliminate 
Ui'(r) (all /') from fl66| ) and (|67j). We then need only work with one of equations ( |66| ) 
or (^) since the equations (^) through ( |d7| ) are related by V a g a = 0. 

We next replace p'/p, Wy, and Vy in equations (|^) or (|67|) by their series ex- 
pansions. We eliminate the Uy (r) from either (^) or (|6^) and, again, substitute for 
the Wi'(r) and V//(r). Finally, we write down the matching condition at the point 
?"o equating the series expansions about r = to the series expansions about r = R. 
The result is a linear algebraic system which we may represent schematically as 

Ax = 0. (88) 

In this equation, A is a matrix which depends non-linearly on the parameter k, and 
x is a vector whose components are the unknown coefficients in the series expansions 
for the Wi'(r) and VJ/(r). In Appendix B, we explicitly present the equations making 
up this algebraic system as well as the forms of the regular series expansions for 
W//(r) and V v (r). 

To satisfy equation (^) we must find those values of k for which the matrix A 
is singular, i.e., we must find the zeroes of the determinant of A. We truncate the 
spherical harmonic expansion of 5v a at some maximum index Z max and we truncate 
the radial series expansions about r = and r = R at some maximum powers z max 
and fc max , respectively. 

The resulting finite matrix is band diagonal. To find the zeroes of its determinant 
we use standard root finding techniques combined with routines from the LAPACK 
linear algebra libraries (Anderson et al. [l|]). We find that the eigenvalues, k , 
computed in this manner converge quickly as we increase Z max , z max and fc max . 
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The eigenfunctions associated with these eigenvalues are determined by the per- 
turbation equations only up to normalization. Given a particular eigenvalue, we 



find its eigenfunction by replacing one of the equations in the system ( p8|) with the 
normalization condition that 

V m (r = R) = 1 for polar- hybrids, or that 
V m+ i(r = R) = 1 for axial-hybrids.) 

Since we have eliminated one of the rows of the singular matrix A in favor of this 
condition, the result is an algebraic system of the form 

Ax = b, (90) 

where A is now a non-singular matrix and b is a known column vector. We solve this 
system for the vector x using routines from LAPACK and reconstruct the various 
series expansions from this solution vector of coefficients. 

2.5 Eigenvalues and Eigenfunctions 

We have computed the eigenvalues and eigenfunctions for uniform density stars and 
for n = 1 polytropes, models obeying the polytropic equation of state p = Kp 2 , where 
K is a constant. Our numerical solutions for the uniform density star agree with the 
recent results of Lindblom and Ipser |5(| who find analytic solutions for the hybrid 



modes in rigidly rotating uniform density stars with arbitrary angular velocity - the 
Maclaurin spheroids. Their calculation uses the two-potential formalism (Ipser and 
Managan |j3"l| ; Ipser and Lindblom |j3"2| ) in which the equations for the perturbation 



modes are reformulated as coupled differential equations for a fluid potential, 5U, 
and the gravitational potential, <5$. All of the perturbed fluid variables may be 
expressed in terms of these two potentials. The analysis follows that of Bryan |T0|| 
who found that the equations are separable in a non-standard spheroidal coordinate 
system. 

The Bryan/Lindblom-Ipser eigenfunctions 5Uq and 5&o turn out to be products 
of associated Legendre polynomials of their coordinates. This simple form of their 
solutions leads us to expect that our series solutions might also have a simple form 
- even though their unusual spheroidal coordinates are rather complicated functions 
of r and 9. In fact, we do find that the modes of the uniform density star have a 



32 



particularly simple structure. For any particular mode, both the angular and radial 
series expansions terminate at some finite indices Zo and i$ (or ko). That is, the 
spherical harmonic expansion ( |50"1) of Sv a contains only terms with m < I < Iq for 
this mode, and the coefficients of this expansion - the Wi(r), Vftr) and Ui(r) - are 
polynomials of order m+i . For all Z > m there exist a number of modes terminating 
at Iq. 

In Tables |IJ to |] we present the functions Wi(r), Vi(r) and Ui(r) for all of the 
axial- and polar-led hybrids with m = 1 and m = 2 for a range of values of the 
terminating index Z . (See also Figure 0.) For given values of m > and Z there are 
lo — m + 1 modes. (When m = there are Zo modes. See Eq. fl9~2] ) below.) We also 
find that the last term in the expansion flSO), the term with Z = Zo, is always axial 
for both types of hybrid modes. This fact, together with the fact that the parity of 
the modes is, 

{(— l) m for polar-led hybrids 
(91) 
(— l) m+ for axial-led hybrids, 

(for m > 0) implies that Zo — m + 1 must be even for polar-led modes and odd for 
axial-led modes. 

The fact that the various series terminate at Zo, zq and ko implies that Equations 
(|88|) and ( j9"0D will be exact as long as we truncate the series at Z max > Z , i max > io 
and Zc max > k . 

To find the eigenvalues of these modes we search the k axis for all of the zeroes 



of the determinant of the matrix A in equation (pq) . We begin by fixing m and 
performing the search with Z max = m. We then increase Z max by 1 and repeat the 
search (and so on). At any given value of Z max , the search finds all of the eigenvalues 
associated with the eigenfunctions terminating at Zo < Z max . 

In Table ||], we present the eigenvalues kq found by this method for the axial- and 
polar-led hybrid modes of uniform density stars for a range of values of Zo and m. 
Observe that many of the eigenvalues, (marked with a *) satisfy the CFS instability 
condition a (a + mfl) < (see Sect. 1.1). The modes whose frequencies satisfy 
this condition are subject to the non-axisymmetric gravitational radiation driven 
instability in the absence of viscosity. The modes having Zo = m > (or Zo = 1 for 
m = 0) are the purely axial r- modes with eigenvalues k = 2/(m+ 1) (or n Q = 
for m = 0) discussed in Sect. 2.3.1. We find that there are no purely polar modes 
satisfying our assumptions in these stellar models. 
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We have compared these eigenvalues with those of Lindblom and Ipser ||50|| . To 



lowest non-trivial order in Q their equation for the eigenvalue, kq, can be expressed 
in terms of associated Legendre polynomials^ (see Lindblom and Ipser's equation 
6.4), as 



( 4 " 4)-K + i{^o) ~ 2mP% +1 (-K ) = 0. (92) 



d 

For given values of m > and Iq this equation has Iq — m + 1 roots (corresponding to 
the number of distinct modes), which can easily be found numerically. (For m = 
there are Iq roots.) For the range of values of m and Iq checked our eigenvalues agree 
with these to machine precision. (Compare our Table § with Table 1 in Lindblom 
and Ipser 0].) 

We have also compared our eigenfunctions with those of Lindblom and Ipser. For 
a uniformly rotating, isentropic star, the fluid velocity perturbation, 5v a , is related 



(Ipser and Lindblom f32| ) to their fluid potential 8U by 

VJU = - [iKttg ab + 2V b v a ] 5v b . (93) 



Since the ip component of this equation is simply 



im5U = —VLr 2 sin 2 9 



-5v r + 2 cot 95v e + ik5v< 
r 



(94) 



it is straightforward numerically to construct this quantity from the components of 
our 5v a and compare it with the analytic solutions for SU given by Lindblom and 
Ipser (see their Eq. 7.2). We have compared these solutions on a 20 x 40 grid in the 
(r — 9) plane and found that they agree (up to normalization) to better than 1 part 
in 10 9 for all cases checked. 

Because of the use of the two-potential formalism and the unusual coordinate sys- 
tem used in their analysis, the axial- or polar-hybrid character of the Bryan/Lindblom- 
Ipser solutions is not obvious. Nor is it evident that these solutions have, as their 
Q — > limit, the zero-frequency convective modes described in Sect. 2.2. The com- 
parison of their analytic results with our numerical work has served the dual purpose 
of clarifying these properties of the solutions and of testing the accuracy of our code. 
The computational differences are minor between the uniform density calculation 
and one in which the star obeys a more realistic equation of state. Thus, this testing 



5 The index I used by Lindblom and Ipser is related to our Iq by I = Iq + 1. Our convention 
agrees with the usual labelling of the Iq — m pure axial modes. 



34 



gives us confidence in the validity of our code for the polytrope calculation. As a 
further check, we have written two independent codes and compared the eigenvalues 
computed from each. One of these codes is based on the set of equations described 
in Appendix B. The other is based on the set of second order equations that results 
from using the mass conservation equation, (0), to substitute for all the Vi(r) in 
favor of the Wi(r). 

For the n — 1 polytrope we will consider and, more generally, for any isentropic 
equation of state, the purely axial r-modes are independent of the equation of state. 
In both isentropic and non-isentropic stars, pure r-modes exist whose velocity field 
is, to lowest order in Q, an axial vector field belonging to a single angular harmonic 
(and restricted to harmonics with I = m in the isentropic case). The frequency of 
such a mode is given (to order Q) by the Papalouizou and Pringle pD| expression, 
Eq. (D), and is independent of the equation of state. As we saw in Sect. 2.3.1, only 
those modes having I = m (or / = 1 for m = 0) exist in isentropic stars, and for 
these modes the eigenfunctions are also independent of the (isentropic) equation of 
state. This independence of the equation of state occurs for the r-modes because (to 
lowest order in Q) fluid elements move in surfaces of constant r (and thus in surfaces 
of constant density and pressure). For the hybrid modes, however, fluid elements 
are not confined to surfaces of constant r and one would expect the eigenfrequencies 
and eigenfunctions to depend on the equation of state. 

Indeed, we find such a dependence. The hybrid modes of the n = 1 polytrope 
are not identical to those of the uniform density star. On the other hand, the modes 
do not appear to be very sensitive to the equation of state. We have found that 
the character of the polytropic modes is similar to the modes of the uniform density 
star, except that the radial and angular series expansions do not terminate. For 
each eigenfunction in the uniform density star there is a corresponding eigenfunction 
in the polytrope with a slightly different eigenfrequency (See Table [|) F° r a given 
mode of the uniform density star, the series expansion ( ^0|) terminates at I — Iq. 
For the corresponding polytrope mode, the expansion ( pT?D does not terminate, but 
it does converge quickly. The largest terms in (|5Up with I > Iq are more than 
an order of magnitude smaller than those with I < Iq and they decrease rapidly 
as I increases. Thus, the terms that dominate the polytrope eigenfunctions are 
those that correspond to the non-zero terms in the corresponding uniform density 
eigenfunctions. 
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In Figures and|| we display the coefficients Wi(r), Vi(r) and Ui(r) of the expan- 
sion ( p0[ ) for the same m = 2 axial-led hybrid mode in each stellar model. For the 
uniform density star (Figure 0) the only non-zero coefficients for this mode are those 
with I < l = 4. These coefficients are presented explicitly in Table |^ and are low 
order polynomials in r. For the corresponding mode in the polytrope, we present in 
Figure |3] the first seven coefficients of the expansion fl50|). Observe that those coef- 
ficients with I < 4 are similar to the corresponding functions in the uniform density 
mode and dominate the polytrope eigenf unction. The coefficients with 4 < / < 6 are 
an order of magnitude smaller than the dominant coefficients and those with I > 6 
are smaller still. (Since they would be indistinguishable from the (r/R) axis, we do 
not display the coefficients having I > 6 for this mode.) 

Just as the angular series expansion fails to terminate for the polytrope modes, so 
too do the radial series expansions for the functions Wi(r), Vftr) and Ui(r). We have 
seen that in the uniform density star these functions are polynomials in r (Tables 
[I] through In the polytropic star, the radial series do not terminate and we are 
required to work with both sets of radial series expansions - those about the center of 
the star and those about its surface - in order to represent the functions accurately 
everywhere inside the star. 



In Figures |] through [12] we compare corresponding functions from each type of 
star. For example, Figures f|, [| and ^| show the functions Wi(r), Vi(r) and Ui(r) 
(respectively) for I < 6 for a particular m = 1 polar-led hybrid mode. In the uni- 
form density star this mode has eigenvalue kq = 1.509941, and in the polytrope it 
has eigenvalue kq = 1.412999. The only non-zero functions in the uniform density 
mode are those with I < Iq — 2 and they are simple polynomials in r (see Table 0). 
Observe that these functions are similar, but not identical to, their counterparts in 
the polytrope mode, which have been constructed from their radial series expansions 
about r = and r = R (with matching point tq = 0.5R). Again, note the conver- 
gence with increasing I of the polytrope eigenf unction. The mode is dominated by 
the terms with I < 2 and those with / > 2 decrease rapidly with I. (The / = 5 and 
I = 6 coefficients are virtually indistinguishable from the {r/R) axis.) 

Because the polytrope eigenf unctions are dominated by their I < Iq terms, the 
eigenvalue search with Z max = l Q will find the associated eigenvalues approximately. 
We compute these approximate eigenvalues of the polytrope modes using the same 
search technique as for the uniform density star. We then increase / max and search 
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near one of the approximate eigenvalues for a corrected value, iterating this procedure 
until the eigenvalue converges to the desired accuracy. We present the eigenvalues 
found by this method in Table |6]. 

As a further comparison between the mode eigenvalues in the polytropic star and 
those in the uniform density star we have modelled a sequence of "intermediate" 
stars. By multiplying the expansions (|86| ) and flS?]) for (p'/p) by a scaling factor, 



e G [0, 1], we can simulate a continuous sequence of stellar ]models connecting the 
uniform density star (e = 0) to the polytrope (e = 1). We find that an eigenvalue 
in the uniform density star varies smoothly as function of e to the corresponding 
eigenvalue in the polytrope. 

2.6 Dissipation 

The effects of gravitational radiation and viscosity on the pure l Q = m r-modes 
discussed in Sect. 2.3.1 have already been studied by a number of authors. (Lindblom 
et al. [[54]], Owen et al. HS9(| , Andersson et al. ||, Kokkotas and Stergioulas f43fl , 



Lindblom et al. p3"|) All of these modes are unstable to gravitational radiation 
reaction, and for some of them this instability strongly dominates viscous damping. 
We now consider the effects of dissipation on the axial- and polar-hybrid modes. 

To estimate the timescales associated with viscous damping and gravitational 
radiation reaction we follow the methods used for the Iq = m modes (Lindblom et al. 



]51| |, see also Ipser and Lindblom [Q). When the energy radiated per cycle is small 
compared to the energy of the mode, the imaginary part of the mode frequency is 
accurately approximated by the expression 

■-ST 

where E is the energy of the mode as measured in the rotating frame, 

d 3 x. (96) 



2 



The rate of change of this energy due to dissipation by viscosity and gravitational 
radiation is, 



dE 



- J {2r]5a ab 5a* ab + (5656* 



dt 

-<j(a + mfi) £ N t a 21 (\5D lrn \ 2 + \5J lrn \ 2 ) . (97) 



l>2 
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The first term in ( |97l) represents dissipation due to shear viscosity, where the shear, 
5a ab , of the perturbation is 

fo ab = ~ (v a 5v b + V b 5v a - ^g ab V c 5v c ^j , (98) 

and the coefficient of shear viscosity for hot neutron-star matter is (Cutler and 
Lindblom |7|; Sawyer §§) 



9 



^ = 2xlO l8 f in1 / X (^) g-cm^-s- 1 . (99) 



s 2 

10 15 g-cm _3 y \^ T 

The second term in (|97|) represents dissipation due to bulk viscosity, where the 
expansion, 56, of the perturbation is 

56 = V c 5v c (100) 
and the bulk viscosity coefficient for hot neutron star matter is (Cutler and Lindblom 



T7|| ; Sawyer |68| 



1Hz \ 2 / p \ 2 / T 



C = 6xl0 25 ( — tt - ^ ? ^T77 g-cm- 1 ^- 1 . (101) 

s \cr + mttj V10 15 g-cm- 3 / \W 9 KJ 6 v ; 

The third term in fl9~T|) represents dissipation due to gravitational radiation, with 
coupling constant 

n = AnG ( / + 1 )( / + 2 ) (l()2) 

1 c 2 '+ 1 /(/-l)[(2/ + l)!!] 2 ' 1 J 

The mass, 5Di m , and current, 5Ji m , multipole moments of the perturbation are given 
by (Thorne [78|, Lindblom et al. |5^| ) 



5 D lm = J 5pr l Y; m d 3 x, (103) 



and 



si 2 I 1 

O-Jlrn 



T jr l (p5v a + 5pv a )Y l a m B *a i x (104) 



c\l + 

where Y^ B is the magnetic type vector spherical harmonic (Thorne |78| ) given by, 

YL* = ~ e ahc V b YrV c r. (105) 

# + 1) 



To lowest order in Vt, the energy (|9T)| ) of the hybrid modes is positive definite. 
Their stability is therefore determined by the sign of the right hand side of equation 
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(|97D . We have seen that many of the hybrid modes have frequencies satisfying the 
CFS instability criterion a (a + mQ) < 0. It is now clear that this makes the third 
term in Eq. (|97|) positive, implying that gravitational radiation reaction tends to 
drive these modes unstable. As discussed in Sect. 1.1, however, to determine the 
actual stability of these modes, we must evaluate the various dissipative terms in 



m 

We first substitute for 5v a the spherical harmonic expansion flSUp and use the 
orthogonality relations for vector spherical harmonics (Thorne f78| ) to perform the 
angular integrals. The energy of the modes in the rotating frame then becomes 

1 r R 

l=m 

To calculate the dissipation due to gravitational radiation reaction we must evalu- 
ate the multipole moments (|103|) and ( |104|) . To lowest order in Q the mass multipole 
moments (|103|) vanish and the current multipole moments are given by 



00 1 r R r i 
£=EW p\w? + l{l + l)V? + l(l + l)U?]dr. (106) 
ZZ, 2 Jo 1 J 



SJim = — / pr^Urfr. (107) 
c Jo 

To calculate the dissipation due to bulk viscosity we must evaluate the expansion, 
59 = V c 5f c , of the perturbation. For uniform density stars this quantity vanishes 
identically by the mass conservation equation (|35|) . For the Iq = m, pure axial 
modes the expansion, again, vanishes identically, regardless of the equation of state. 
To compute the bulk viscosity of these modes it is necessary to work to higher order 
in Q (Andersson et al. P], Lindblom et al. El). On the other hand, for the new 



hybrid modes in which we are interested, the expansion of the fluid perturbation is 
non-zero in the slowly rotating polytropic stars. After substituting for Sv a its series 
expansion and performing the angular integrals, the bulk viscosity contribution to 
(P?D becomes 



'dE\ ^ f R C 

~~dl 



B l=m 



E / \ Wi + w i - l ( l + dr ( 108 ) 

f „ Jo r 
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In a similar manner, the contribution to from shear viscosity becomes 



'dE\ 
, dt )s 



00 rR 

Jo r 2 1 3 

(=771 ^ 



277/2 



r — ~ 



1(1 + l)W 2 



+ -/(Z + l)(2Z 2 + 2/-3)Vf 



(109) 



+1(1 + l)Wi 



Vi 



5 1 y i 



+ -l(l + l)V l (rWi)' 



--1(1 
2 v 



1 2 



+ h(l + l)(l 2 + l-2)uAdr. 



Given a numerical solution for one of the hybrid mode eigenf unctions, these 
radial integrals can be performed numerically The resulting contributions to ( |97| ) 
also depend on the angular velocity and temperature of the star. Let us express the 
imaginary part of the hybrid mode frequency (|)5|) as, 



1 

T 



1 



'l0 9 iT 
T 



1 

+ — 



T 



10 9 fT 



nGp / 



1+1 



(110) 



where p is average density of the star. (Compare this expression to the corresponding 
expression in Lindblom et al. j54] - their equation (22) - for the Iq = m pure axial 
modes.) 

The bulk viscosity term in this equation is stronger by a factor Q~ A than that for 
the lo = m pure axial modes. This is because the expansion 58 of the hybrid mode 
is nonzero to lowest order in Q for the polytropic star, whereas it is order Q 2 for the 
pure axial modes. This implies that the damping due to bulk viscosity will be much 
stronger for the hybrid modes than for the pure axial modes in slowly rotating stars. 

Note that the contribution to ( |110| ) from gravitational radiation reaction consists 
of a sum over all the values of / with a non-vanishing current multipole. This sum 
is, of course, dominated by the lowest contributing multipole. 

In Tables [7] to |^ we present the timescales for these various dissipative effects 
in the uniform density and polytropic stellar models that we have been considering 
with R = 12.57km and M = 1AM Q . For the reasons discussed above, we do not 
present bulk viscosity timescales for the uniform density star. 



40 



Given the form of their eigenfunctions, it seems reasonable to expect that some 
of the unstable hybrid modes might grow on a timescale which is comparable to 
that of the pure Iq = m r-modes. For example, the m = 2 axial-led hybrids all 
have U 2 (r) ^ (see, for example, Figures and |3|). By equation (|107|) , this leads 



one to expect a non-zero current quadrupole moment 6 J22, and this is the multipole 
moment that dominates the gravitational radiation in the r-modes. Upon closer 
inspection, however, one finds that this is not the case. In fact, we find that all of 
the multipoles 5Ji m vanish (or nearly vanish) for I < l , where Iq is the largest value 
of I contributing a dominant term to the expansion ([50]) of 5v a . 

In the uniform density star, these multipoles vanish identically. Consider, for 
example, the m = 2, Iq = 4 axial-hybrid with eigenvalue k = 0.466901. (See Table 
[3] and Figure [7]) For this mode, f/ 2 oc (7x 3 — 9x 5 ), where x = (r/R). By equation 



(|107|) , we then find that 



SJ22 oc / x 6 (7x A - 9x b )dx = 0, (111) 
Jo 

and that dJ^ is the only non-zero radiation multipole. In general, the only non-zero 
multipole for an axial- or polar-hybrid mode in the uniform density star is SJi om . 

That this should be the case is not obvious from the form of our eigenfunctions. 
However, Lindblom and Ipser's ]50| analytic solutions provide an explanation. Their 



equations (7.1) and (7.3) reveal that the perturbed gravitational potential, 5$, is a 
pure spherical harmonic to lowest order in Q. In particular, 

M> oc Y™ +1 . (112) 

This implies that the only non-zero current multipole is SJi om . 

We find a similar result for the polytropic star. Because of the similarity between 
the modes in the polytrope and the modes in the uniform density star, we find that 
although the lower I current multipoles do not vanish identically, they very nearly 
vanish and the radiation is dominated by higher / multipoles. 

The fastest growth times we find in the hybrid modes are of order 10 4 seconds 
(at lO 9 ^ and Q = VnGp). Thus, the spin-down of a newly formed neutron star will 
be dominated by the Iq = m = 2 mode with small contributions from the Iq = m 
pure axial modes with 2 < m ^ 10 and from the fastest growing hybrid modes. 
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-7a; 4 ) 


-0.646064a; 4 





-3(a- 2 - a- 4 ) 





-1.5a; 2 + 2.5a; 4 
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1.728851(5ai 2 


-7a; 4 ) 


1.431460x 4 





-3(a- 2 - a- 4 ) 





-1.5a; 2 + 2.5a; 4 





1.708024 
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-7a; 4 ) 
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-3(a- 2 - a- 4 ) 





-1.5a; 2 + 2.5a; 4 
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- 31.5a; 4 
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5.25a; 2 - 13.5a; 4 
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2.625a; 2 - 11.25a; 4 


-0.490203(5ai 4 




+24. 75a; 6 ; 


1 
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+8. 25a; 6 


-20a- 6 ) 


+9. 625a; 6 


-7a- 6 ) 
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a The eigenfunctions are normalized so that V2 = 1 at the surface of the star, x = 1. Here x = (r/R). 

b lo — m + 1 = Iq is the maximum value of / in the spherical harmonic expansion of Sv a . Observe that this I = Iq term is always axial. 
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a The eigenfunctions are normalized so that Vi — 1 at the surface of the star, x — 1. Here x — (r/R). 

h lo — m + 1 = Iq is the maximum value of / in the spherical harmonic expansion of 5v a . Observe that this / = Zq term is always axial. 
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a The eigenfunctions are normalized so that V3 — 1 at the surface of the star, x — 1. Here x — (r/R). 

h lo — m + 1 = /q — 1, where Zq is the maximum value of / in the spherical harmonic expansion of 8v a . Observe that this I = lo term is always axial. 
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a The eigenfunctions are normalized so that V2 — 1 at the surface of the star, x — 1. Here x — (r/R). 

b /o — m + 1 = Iq — 1, where Zo is the maximum value of / in the spherical harmonic expansion of 5v a . Observe that this I = Iq term is 
always axial. 
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TABLE 5 
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a Kofl = (a + mfl) is the mode frequency in the rotating frame to lowest order in fi. The 
modes whose frequencies are marked with a * satisfy the condition a (a + mft) < and are 
subject to a gravitational radiation driven instability in the absence of viscous dissipation. 

b For m = 0, this is simply l . For the uniform density star, l is the maximum value of I 
appearing in the spherical harmonic expansion of Sv a . 

c This denotes the parity class of the mode; a meaning axial-led hybrids, and p meaning 
polar-led hybrids. 

d These are the eigenvalues of the pure lo = m r-modes. For isentropic stars they are 
independent of the equation of state and have the value /to = 2/(m + 1) (or kq = for 
to = 0) to lowest order in ft. (See Sect. 2.3.1) 



Table 5: Eigenvalues k for Uniform Density Stars. 
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TABLE 6 



(lo — m + l) b parity m = m = 1 m = 2 m = 3 m = 4 



l d 


a 





,000000 


1 


,000000 


0. 


.666667* 


0. 


.500000* 





.400000* 


2 


P 


-1 


.028189 


-0 


.401371 


-0 


.556592 


-0 


.631637 


-0 


.672385 




P 


1 


.028189 


1 


.412999 


1. 


.100026* 


0. 


.904910* 





.771078* 


3 


a 


-1 


.358128 


-1 


.032380 


-1 


.025883 


-1 


.014866 


-1 


.002175 




a 


0. 


,000000 





.690586* 





,517337* 


0. 


.412646* 





.342817* 




a 


1 


,358128 


1 


.613725 


1 


,357781* 


1. 


.176745* 


1 


.041683* 


4 


P 


-1 


.542065 


-1 


.312267 


-1 


.272885 


-1 


.238631 


-1 


.208390 




P 


-0 


.701821 


-0 


.178792 


-0. 


.275335 


-0 


,333267 


-0 


.370450 




P 


0. 


.701821 


1 


.051525 





.862948* 





,734297* 





.640592* 




P 


1. 


.542065 


1 


.726257 


1 


.519573* 


1. 


.360560* 


1 


.234698* 


5 


a 


-1 


,656481 


-1 


.483402 


-1. 


.433916 


-1 


.391305 


-1 


.354057 




a 


-1 


.013703 


-0 


.705182 


-0 


.703898 


-0. 


.699942 


-0 


.694498 




a 


0. 


.000000 





.528102* 





.421678* 





.350192* 





.299055* 




a 


1 


.013703 


1 


.281962 


1 


.104402* 


0. 


.974192* 





.874124* 




a 


1. 


.656481 


1 


.795734 


1. 


.627215* 


1 


.489441* 


1 


.375406* 



a Kofi = (<r + mfi) is the mode frequency in the rotating frame to lowest order in fi. The 
modes whose frequencies are marked with a * satisfy the condition a (a + mVt) < and are 
subject to a gravitational radiation driven instability in the absence of viscous dissipation. 

b For m = 0, this is simply l . For the n = 1 polytrope, Iq is the largest value of I that 
contributes a dominant term to the spherical harmonic expansion of Sv a . 

c This denotes the parity class of the mode; a meaning axial-led hybrids, and p meaning 
polar-led hybrids. 

d These are the eigenvalues of the pure lo = m r-modes. For isentropic stars they are 
independent of the equation of state and have the value /to = 2/(m + 1) (or k = for 
m = 0) to lowest order in fl (See Sect. 2.3.1). 



Table 6: Eigenvalues k for the p = Kp 2 Polytrope. 



TABLE 7 



lo 


n b 


K 


f B c 


fs 


h 




3 





0.611985 




7.67 x 10 7 


-9.79 x 10 6 






1 


0.690586 


5.86 x 10 9 


9.29 x 10 7 


-1.25 x 10 8 


-1.22 x 10 20 


5 





0.440454 




2.04 x 10 7 


— 00 


-2.07 x 10 13 




1 


0.528102 


2.57 x 10 9 


3.87 x 10 7 


-2.17 x 10 10 


-5.75 x 10 14 



a We present dissipative timescales only for those modes that are unstable to 
gravitational radiation reaction. None of the m = 1 polar-hybrid modes are 
unstable for low values of lo. 

b The polytropic index, n, where p = Kp 1+1 / n . The n=0 case represents the 
uniform density equilibrium star. 

c Dissipation due to bulk viscosity is not meaningful for uniform density stars. 



Table 7: Dissipative timescales (in seconds) for m — 1 axial-hybrid modes at T 
10 9 K and n = ^Gp. 
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TABLE 8 



Co - 1) 


n & 














T 2 


T4 




f 




l c 








.666667 




4 


.46 


X 


10 s 


-1.56 x 10° 














1 





.666667 


2.0 x 10 11 


2 


.52 


X 


10 s 


-3.26 x 10° 












3 








.466901 




4 


.10 


X 


10 7 


— oo 


-3.88 x 10 5 












1 





.517337 


6.43 x 10 9 


6 


.21 


X 


10 7 


< -10 18 


-1.85 x 10 6 


-4. 


.97 


X 


10 16 







1 


.496436 




3 


.92 


X 


10 7 


— oo 


-5.85 x 10 9 












1 


1 


.357781 


4.10 x 10 9 


7 


.18 


X 


10 7 


< -10 19 


-1.60 x 10 9 


-4. 


.35 


X 


10 19 


5 








.359536 




1 


.34 


X 


10 7 


— oo 


— oo 


-1. 


.28 


X 


10" 




1 





.421678 


2.65 x 10 9 


3 


.01 


X 


10 7 


< -10 16 


-2.01 x 10 9 


-1 


.15 


X 


10 12 







1 


.153058 




1 


.32 


X 


10 7 


— oo 


— oo 


-3 


.11 


X 


10 14 




1 


1 


.104402 


2.45 x 10 9 


3 


.65 


X 


10 7 


< -10 12 


-1.37 x 10 11 


-4. 


.89 


X 


10 14 







1 


.733971 




1 


.31 


X 


10 7 


— oo 


— oo 


-1. 


.92 


X 


10 21 




1 


1 


.627215 


5.32 x 10 9 


3 


.44 


X 


10 7 


< -10 19 


-2.30 x 10 15 


-8. 


.33 


X 


10 19 



a The polytropic index, n, where p = Kp 1+1 /". The n=0 case represents the uniform density 
equilibrium star. 

b Dissipation due to bulk viscosity is not meaningful for uniform density stars. 

c This is the Iq = m = 2 r-mode already studied by Lindblom et al. (1998), Owen et al. (1998), 
Andersson et al. (1998), Kokkotas and Stergioulas (1998) and Lindblom et al. (1999). The value 
of the bulk viscosity timescale for this mode is taken from Lindblom et al. (1999) who calculate it 
self-consistently using an order CI 2 calculation. 



Table 8: Dissipative timescales (in seconds) for m = 2 axial-hybrid modes at T = 
10 9 K and n = y/WGp. 



TABLE 9 



(*o - 1) 


n a 


K 




fs 




f 3 




2 





1.231925 




9.03 x 


10 7 


-4.77 x 10 4 






1 


1.100026 


3.32 x 10 9 


1.24 x 


10 8 


-3.37 x 10 4 


-3.13 x 10 14 


4 





0.884249 




2.17 x 


10 7 


— 00 


-5.64 x 10 9 




1 


0.862948 


1.93 x 10 9 


4.94 x 


10 7 


-1.10 x 10 7 


-1.45 x 10 10 







1.643443 




2.13 x 


10 7 


— 00 


-2.12 x 10 15 




1 


1.519573 


4.79 x 10 9 


4.77 x 


10 7 


-1.92 x 10 11 


-2.31 x 10 14 



a The polytropic index, n, where p = Kp 1+1 l n . The n=0 case represents the uniform 
density equilibrium star. 

b Dissipation due to bulk viscosity is not meaningful for uniform density stars. 



Table 9: Dissipative timescales (in seconds) for m = 2 polar-hybrid modes at T 
10 9 K and n = ^Gp. 
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Figure 1: Images of the perturbed velocity field 5v a for a few of the newtonian r- 
modes at a fixed time. The velocity profile shown rotates forward in the inertial 
frame with angular velocity (/ — + 2)Q/l(l + 1) and backward in the rotating 
frame with angular velocity 2Q/l(l + 1). Fluid elements oscillate in small circles in 
the rotating frame as their velocity changes in accordance with this rotating profile. 
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-1 - 



G- 



B U 2 (r) 
W 3 (r) 
V 3 (r) 
U 4 (r) 



'G . ,G 



0.2 



0.4 0.6 

(r/R) 



0.8 



Figure 2: All of the non-zero coefficients Wi(r), VJ(r), t/j(r) of the spherical harmonic 
expansion ([50]) for a particular m = 2 axial-led hybrid mode of the uniform density 
star. The mode has eigenvalue kq = —0.763337. Note that the largest value of / that 
appears in the expansion floTf ) is Iq = 4 and that the functions Wi(r), V/(r) and Ui(r) 
are low order polynomials in (r/R). (See Table |3|.) The mode is normalized so that 
V 3 (r = R) = l. 
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_ 2 I , 1 , 1 

0.5 1 

(r/R) 

Figure 3: The coefficients Wi(r), Vftr), Ui(r) with I < 6 of the spherical harmonic 
expansion ([50]) for a particular m = 2 axial-led hybrid mode of the polytropic star. 
This is the polytrope mode that corresponds to the uniform density mode displayed 
in Figure 0. For the polytrope the mode has eigenvalue k = —1.025883. The 
expansion (50) converges rapidly with increasing I and is dominated by the terms 
with 2 < I < 4, i.e., by the terms corresponding to those which are non-zero for 
the uniform density mode. Observe that the coefficients shown with 4 < / < 6 are 
an order of magnitude smaller than those with 2 < I < 4. Those with / > 6 are 
smaller still and are not displayed here. The mode is, again, normalized so that 
V 3 (r = R) = l. 
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0.5 



Wi(r) for the uniform density star 
Wi(r) for the n = 1 polytrope 
W3(r) for the n = 1 polytrope 
W5(r) for the n = 1 polytrope 




-0.5 



0.5 

(r/R) 



Figure 4: The functions Wi(r) with Z < 6 for a particular m = 1 polar- led hybrid 
mode. For the uniform density star this mode has eigenvalue = 1.509941 and 
W\ = —x + x 3 (x = r / R) is the only non- vanishing Wi(r) (see Table El). The cor- 
responding mode of the polytropic star has eigenvalue kq = 1.412999. Observe that 
Wi(r) for the polytrope, which has been constructed from its power series expansions 
about r = and r = R, is similar, though not identical, to the corresponding Wi(r) 
for the uniform density star. Observe also that the functions Wi(r) with Z > 1 for 
the polytrope are more than an order of magnitude smaller than Wi(r) and become 
smaller with increasing Z. {W^(r) is virtually indistinguishable from the {r/R) axis.) 



54 



0.5 







1 1 




/ 


Vi(r) 


for the uniform density star 
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for the n = 1 polytrope 
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for the n = 1 polytrope 




V 5 (r) 
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Figure 5: The functions V/(r) with / < 6 for the same mode as in Figure |j. 
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U2(r) for the uniform density star 
U%{r) for the n = 1 polytrope 
U^r) for the n = 1 polytrope 
Ue(r) for the n = 1 polytrope 




-0.1 
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Figure 6: The functions Ui(r) with Z < 6 for the same mode as in Figure £§ 



U2(r) for the uniform density star 

£72 (r) for the n = 1 polytrope 

3 £74 (r) for the uniform density star 

£74 (r) for the n = 1 polytrope 

£76 (r) for the n = 1 polytrope 



0.2 0.4 0.6 0.8 1 

(r/R) 

Figure 7: The functions £7;(r) with Z < 7 for a particular m = 2 axial-led hybrid 
mode. For the uniform density star this mode has eigenvalue k = 0.466901 and 
£72 (r) and £74(r) are the only non- vanishing Ui(r). (See Table |3] for their explicit 
forms.) The corresponding mode of the polytropic star has eigenvalue k,q = 0.517337. 
Observe that £72 (r) and £74 (r) for the polytrope, which have been constructed from 
their power series expansions about r = and r = R, are similar, though not 
identical, to the corresponding functions for the uniform density star. Observe also 
that £76 (r) is more than an order of magnitude smaller than £72 (r) and £74 (r). 



-6 



-8 
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1.5 



W3(r) for the uniform density star 
W3 (r) for the n = 1 polytrope 
W$(r) for the n = 1 polytrope 
Wyir) for the n = 1 polytrope 
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Figure 8: The functions Wi(r) with / < 7 for the same mode as in Figure 0. 
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(r/R) 



Figure 9: The functions Vj(r) with I < 7 for the same mode as in Figure [7|. 
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Figure 10: The functions Ui(r) with / < 8 for a particular m = 2 axial-led hybrid 
mode. For the uniform density star this mode has eigenvalue k = 0.359536 and 
U2(r), U±(r) and Uq{t) are the only non- vanishing Ui(r). (See Table ^ for their explicit 
forms.) The corresponding mode of the polytropic star has eigenvalue k,q = 0.421678. 
Observe that U^r), U^{r) and Uq{t) for the polytrope, which have been constructed 
from their power series expansions about r = and r = R, are similar, though not 
identical, to the corresponding functions for the uniform density star. Observe also 
that U$(r) is more than an order of magnitude smaller than U2{r), U^(r) and Uq{t). 
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Figure 11: The functions Wi(r) with I < 8 for the same mode as in Figure [T(| 
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Figure 12: The functions Vi(r) with I < 8 for the same mode as in Figure |10 
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Chapter 3 

Relativistic Stars: Analytic 
Results 

In Ch. 2 we examined the rotationally restored hybrid modes of isentropic newtonian 
stars. We now consider the corresponding modes in relativistic stars. As in the 
newtonian case, we must begin by examining the perturbations of the non-rotating 
star, and finding all of the modes belonging to its degenerate zero-frequency subspace. 

3.1 Stationary Perturbations of Spherical Stars 

The equilibrium of a spherical perfect fluid star is described by a static, spherically 
symmetric spacetime with metric g af3 of the form, 

rfs 2 = _ e 2v{r) dt 2 + e 2X { r) dr 2 + ^2 + ^2^ (m) 

and with energy-momentum tensor, 

T a( 5 = (e + p)u a Ufs + pg afS , (114) 
where e(r) is the total fluid energy density, p(r) is the fluid pressure and 

u a = e -u t a ( n5 ) 

is the fluid 4-velocity - with t a = (d t ) a the timelike Killing vector of the spacetime. 
These satisfy an equation of state of the form 

p = p(e) (116) 
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as well as the Einstein equations, G Q/ g = 87rT aj g, which are equivalent to 

dp _ (e + p)(M + 4irr 3 p) 
dr ~ r(r - 2M) ' 



(117) 



dr 
and 

dv 1 dp 

dr (e + p) dr 

where 



4vrr 2 e (118) 



(119) 



M(r) = \r(l- e~ 2X ). (120) 



(See, e.g., Wald @, Ch.6.) 

We are, again, interested in the space of zero-frequency modes, the linearized, 
time-independent perturbations of this static equilibrium. As in the newtonian case, 
we find that this zero-frequency subspace is spanned by two classes of perturbations. 
To identify these classes explicitly, we must examine the equations governing the 
perturbed configuration. 

Writing the change in the metric as h a p = 5g a/ 3, we express the perturbed con- 
figuration in terms of the set (h a / 3 ,5u a ,5e,5p). These must satisfy the perturbed 
Einstein equations 5G£ = SttSTJ 3 , together with an equation of state (which may, 
in general, differ from that of the equilibrium configuration). 

The perturbed Einstein tensor is given by 

6Gj> = ~l V 7 VV - V 7 V^V - VW a /i/ + V a V?h (121) 

+2R2hf + (V a V^h Ql3 - V 7 VU - R ap h ap )g£ 

where h = /i Q a , V Q is the covariant derivative associated with the equilibrium metric 
and 

RJ? = 8vr(T/ - l -Tg£) = 8vr [(e + p)u a u? + ±(e- p)gj\ (122) 
is the equilibrium Ricci tensor. The perturbed energy-momentum tensor is given by 

6Tf = (5e + 6p)u aU P + 6p6f + (e + p)5u a vP + (e + p)u a 6u . (123) 



Following Thorne and Campolattaro |[79|| , we expand our perturbed variables in 
scalar, vector and tensor spherical harmonics. The perturbed energy density and 
pressure are scalars and therefore must have polar parity, (—1)'. 

5e = 5e(r)Yr, (124) 
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5p = 6p{r)XT. 

The perturbed 4-velocity for a polar-parity mode can be written 



5u% 



hl^Y^e + ^-W{r)Y™r a + \/(r)V Q Y, m \< " 



while that of an axial-parity mode can be written 

-?7(r)e (A - I/) e a/3 ^V /3 y; m M 7 V 5 r. 



5u a A 



(125) 



(126) 



(127) 



(We have chosen the exact form of these expressions for later convenience.) 

To simplify the form of the metric perturbation we will, again, follow Thorne and 
Campolattaro [J7j| and work in the Regge- Wheeler [B5| gauge. The metric perturba- 
tion for a polar-parity mode can be written 



hi/ 



H (r)e 
Hx(r 





2;/ 



#i(r) 
H 2 (r)e 2X 








r 2 K(r) 







r 2 sin 2 6K(r) 



:i28i 



while that of an axial-parity mode can be written 






symm symm 

symm symm 








h (r) sin Qd e Y{ 
hi(r) sin 6deY{ 







(129) 



The Regge- Wheeler gauge is unique for perturbations having spherical harmonic 
index / > 2. However, when Z = 1 or Z = 0, there remain additional gauge degrees of 
freedom^. In addition, the components of the perturbed Einstein equation acquire 
a slightly different form in each of these three cases (cf. Campolattaro and Thorne 
[ p"Tf ) and will be presented separately below. 

We have derived these components using the Maple tensor package by substi- 
tuting expressions (|124|) -( |i~29| ) into Eqs. (|121|) and (|123|) (making liberal use of the 
equilibrium equations (|117|) through ( |120| ) to simplify the resulting expressions). 

1 Letting eAB be the metric on a two-sphere with cab and Da the associated volume element and 
covariant derivative, respectively, one finds the following. When I > 2 the polar tensors DaDbY" 1 
and cabY™ are linearly independent, but when I = 1, they coincide. In addition, the axial tensor 

' Dq)DbY™ vanishes identically for I = 1 and, of course, DaY™ vanishes for I = 0. 
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The resulting set of equations for the case I > 2 are equivalent to those presented in 
Thorne and Campolattaro W9\ upon specializing their equations to the case of sta- 
tionary perturbations and making the necessary changes of notation^. Similarly, the 
set of equations for the case / = 1 are equivalent to those presented in Campolattaro 



and Thorne 11 



3.1.1 The case / > 2 

The non-vanishing components of the perturbed Einstein equation for / > 2 are as 
follows. We will use Eq. ( |133| ) below, to replace H 2 by H . From 5G t l = 87r5T t * we 
have 

1 







e- 2 \ 2 K" + e" 2A (3 - r\')rK' - 



-1(1 + 1) 



K 



-e~ 2X rH' Q 



-1(1 + 1) + l-8nr 2 



H + 8nr 2 5e. 



From 6G r r = 8ir5T r r we similarly have 
= e _2A (l +rv')rK' - 



-1(1 + 1) 



K 



(130) 



-e- 2X rH 



-1(1 + 1) - 1 - 8nr 2 p 



H — 8nr 2 5p. 



From SGf + 5G* = 8n (5T e e + 5T/) we have 

= e' 2X r 2 K" + e~ 2X [r(u - A') + 2] rK' - 16nr 2 5p 



e- 2X r 2 H' ' - e" 2A (3rz/' - r\' + 2)rH' - lQ7rr 2 pH . 



From SG e ° - 5G* = 8vr (5T e ° - ST/) we have 



Ho = Hn- 



From 5Gf = 8it5Tj we have 



(131) 



(132) 



(133) 



K' = e 



~2u 



e li 



(134) 



2 In particular, their fluid variables (denoted by the subscript tc) are related to ours as follows: 
U TC (r,t) = Ut, W T c(r,t) = ~re x Wt, V TC (r,t) = Vt and Se/(e+p) = Sp/jp = -(K + \H ). Their 
equilibrium metric has the opposite signature and differs in the definitions of the metric potentials 
v tc — \ v and Atc = 5 A. 
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From 8G t r = 8n8T t r we have 

From 8G t = 8n5T t 9 we have 

= e~ {u ~ x) [e^-^H^ + 16vr(e + p)e 2X V. 

From 8G t 9 = 8ir5T t ^ we have 

"(2-/ 2 -/) , A 2, , 



/i - [y 1 + A')/i' + 



From SGjf = 8n6T* we have 



-e— (v' + A') - 



(/- l)(Z + 2)/i 1 = 0. 



h = -(*/ + A')C/. 
r 



Finally, from 5G e ' p = 8tt5T ip we have 



3 -("-A) 



(135) 



(136) 



(137) 



;i38) 



(139) 



3.1.2 The case / = 1 

The I = 1 case differs from / > 2 in two respects (Campolattaro and Thorne [|TT|]). 
Firstly, H 2 (r) ^ H (r), because the equation 8G e 6 - 8G* = 8ir (dT e e - van- 
ishes identically. Secondly, we may exploit the aforementioned gauge freedom for 
this case to eliminate the metric functions K(r) and hi(r). (We note that Eq.( |138| ) 
implies hi(r) = for / > 2 anyway.) With these two differences taken into account 
the non-vanishing components of the perturbed Einstein equation for I = 1 are as 
follows. From 5G t l = 87r<5T t * we have 



= e~ 2X rH' 2 + (2 - 87rr 2 e) H 2 - 8nr 2 5e. 
From 5G r r = 87i5T r r we have 

= e~ 2X rH - H + (l + 8nr 2 p) H 2 + 8ixr 2 5p. 



(140) 



(141) 



From 5Gf + 5G„ 



ir [8TJ + 8TJ) we have 



e- 2X r 2 H' Q ' + e~ 2X {2rv' - r\' + l)rH' - H 
+e~ 2X (l + rv')rH 2 + (1 + lQ7rr 2 p)H 2 + 16irr 2 8p 



(142) 



From 8GJ? = 8n5Tj we have 
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= rH' + [rv' - 1)H + (rV + l)H 2 



From SG t r = 8ir5T t r we, again, have 



(143) 



(144) 



= H 1 + 87r(e + p)e 2A rW. 
From 5G t e = 8n5T t 6 we, again, have 

= e^ u ~ x) [e (l/ - A) #i] ' + 167r(e + p)e 2X V. (145) 
Finally, from SG^ = SttST^ we have 

h'o - (i/ + A') /id - [J(i/ + A') + ^Uo = -y + A')C/. (146) 

3.1.3 The case / = 

The / = case differs yet again from the previous two, being the case of stationary, 
spherically symmetric perturbations of a static, spherical equilibrium. To maximize 
the similarity to the preceding two cases we will use the same form for the perturbed 
metric except that we may now exploit the gauge freedom for this case to elimi- 
nate the functions K(r), Hi(r) and hi(r). The non- vanishing components of the 
perturbed Einstein equation for I — are as follows. From SG^ = 8rr5T t 1 we have 



= e~ 2X rH' 2 + (l - 87rr 2 e) H 2 - 8nr 2 5e. 
From 5G r r = 8n5T r r we have 

= e - 2X rH' + (l + 8vrr 2 p) H 2 + 8nr 2 5p. 

From 5G e e + 8G* = 8tt (6T e + 5T/) we have 

= e - 2 V#o + e- 2X (2ru' - rX' + l)rH' 

+e" 2A (l + rv')rH' 2 + lQirr 2 pH 2 + lQirr 2 5p 



(147) 



(148) 



(149) 



Finally, from 5G t r = 8n5T t r we have 



= I6n(e + p)W. 



(150) 
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3.1.4 Decomposition of the zero-frequency subspace. 

By inspection of the above three sets of equations, it is evident that they decouple 
into two independent classes. For I > 2 Eqs. (|130|) -( |T34] ) involve only the variables 
(H , H 2 , K, 5e, 5p) while Eqs. (p5|) - (p7|) involve only the variables (H u h , W, V, U) 
(with hi(r) = implied by Eq. (p8j)). Similarly for I = 1 Eqs. ([140[)-([143D involve 
only (H ,H 2 ,5e,5p) while Eqs. (P^)-(|HBD involve only (H u h , W, V, U) and, in 
fact, are identical to Eqs. (p^-fpTp. Finally, for I = Eqs. (g47D-(p9|) involve 
only (H , H 2 , Se, dp) while Eq. ( p.50| ) simply implies that W{r) = 0. 
Thus, any solution, 

(Ho, H h H 2 , K, h , W, V, U, 5e, dp), (151) 

to the equations governing the time-independent perturbations of a static, spherical 
star is a superposition of (i) a solution 

(o,ifi,o,o,/io,w;v;tf,o,o) (152) 

to Eqs. (P5D-(P?1) or (pop and (ii) a solution 

(^0,0,^2,^,0,0,0,0,56,^) (153) 

to Eqs. O-O, or 0-(H). 

For the solutions of type (ii), the vanishing of the (tr), (t9) and (tip) components 
of the perturbed metric in our coordinate system implies that these solutions are 
static. If, as in the newtonian case, one assumes the linearization stability^ of these 
solutions, i.e., that any solution to the static perturbation equations is tangent to 
a family of exact static solutions, then the theorem that any static self-gravitating 
perfect fluid is spherical implies that any solution of type (ii) is simply a neighboring 
spherical equilibrium. 

Thus, under the assumption of linearization stability we have shown that all 
stationary non-radial (/ > 0) perturbations of a spherical star have 

H = H 2 = K = 6e = 5p = 

3 Again, we are aware of a proof of this linearization stability property under assumptions on 
the equation of state that are satisfied by uniform density stars, but would not allow polytropes 
(Kunzle and Savage Q). 
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and satisfy Eqs. ( |135| )-( [L371) ; that is 
= Ht 







h'o - (v 1 + X')h' + 



>-A) 



1(1 + 1) 



e^-^HA' + 16vr(e + j9)e 2 V, 



(154) 



(155) 



(2-Z 2 



0^ 



2 2 

V+Ao-4 



4 



/U, = -(*/ + A') U. (156) 



Observe that if we use Eq. ( |154| ) to eliminate H\{r) from Eq. (|155 ) we obtain 



-(i/+A) 



[e+p)e v+x rW 



(157) 



J(J + l)(e + p) 

This equation is clearly the generalization to relativistic stars of the conservation 
of mass equation in newtonian gravity, Eq. (|40|) from Sect. 2.1. The other two 
equations relate the perturbation of the spacetime metric to the perturbation of the 
fluid 4- velocity and vanish in the newtonian limit. 

These perturbations must be regular everywhere and satisfy the boundary con- 
dition that the lagrangian change in the pressure vanish at the surface of the star. 
(See Sect. 3.4 below.) As with newtonian stars, this boundary condition requires 
only that 

W(R) = 0. (158) 

leaving W(r) and U(r) otherwise undetermined. If W(r) and U(r) are specified, 
then the functions H\{r), ho(r) and V(r) are determined by the above equations. 
(These solutions are subject to matching conditions to the solutions in the exterior 
spacetime, which must also be regular at infinity. See Sect. 3.4.) 

Finally, we consider the equation of state of the perturbed star. For an adiabatic 
oscillation of a barotropic star (i.e., a star that satisfies a one-parameter equation of 
state, p = p(e)) Eq. ( |2"9"D implies that the perturbed pressure and energy density are 
related by 



dp 

IP 



Se 



(e + p) 



(e + p) 



p 
IP 



for some adiabatic index j(r) which need not be the function 

(e + p) dp 



r r 



V 



de 



(159) 



(160) 
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associated with the equilibrium equation of state. Here, £ Q is the lagrangian dis- 
placement vector and is related to our perturbation variables by 

q a p£ u f = 5u a - \u a u p u^h Pl (161) 

Thus, we have 

e- v dtC = Su r (162) 
or (taking the initial displacement (at t — 0) to be zero) 

f = te u 5u T . (163) 

For the class of perturbations under consideration, we have seen that Sp = 5e = 0, 
thus Eqs. (|159|) and ( |163|) require that 



5u r 



P 



(e + p) 7p 



0. (164) 



For axial-parity perturbations this equation is automatically satisfied, since Su'X has 
no r-component (Eq. ( |127| )). Thus, a spherical barotropic star always admits a class 
of zero-frequency r-modes. 

For polar-parity perturbations, 5u r P = e~ u W(r)/r ^ 0, and Eq. ( |164| ) will be 
satisfied if and only if 

7r)=T(r) = — . (165) 

p de 

Thus, a spherical barotropic star admits a class of zero-frequency g-modes if and only 
if the perturbed star obeys the same one-parameter equation of state as the equilib- 
rium star. Once again, we will call such a star isentropic, because isentropic models 
and their adiabatic perturbations obey the same one-parameter equation of state. 
(That all axial-parity fluid perturbations of a spherical relativistic star are time- 
independent was shown by Thorne and Campolattaro [f79|l . The time-independent 
g-modes in spherical, isentropic, relativistic stars were found by Thorne [f76| .) 

Summarizing our results, we have shown the following. A spherical barotropic 
star always admits a class of zero-frequency r-modes (stationary fluid currents with 
axial parity); but admits zero-frequency g-modes (stationary fluid currents with polar 
parity) if and only if the star is isentropic. Conversely, the zero-frequency subspace 
of non-radial perturbations of a spherical isentropic star is spanned by the r- and 
g-modes - that is, by convective fluid motions having both axial and polar parity 



71 



and with vanishing perturbed pressure and density. Being stationary, these r- and 
g-modes do not couple to gravitational radiation, although the r-modes do induce 
a nontrivial metric perturbation (h t g, h tip 7^ 0) in the spacetime exterior to the star 
(frame-dragging). One would expect this large subspace of modes, which is degen- 
erate at zero-frequency, to be split by rotation, as it is in newtonian stars, so let us 
now consider the perturbations of slowly rotating relativistic stars. 



3.2 Perturbations of Slowly Rotating Stars 

The equilibrium of an isentropic, perfect fluid star that is rotating slowly with uni- 
form angular velocity Q is described (Hartle P9| , Chandrasekhar and Miller ]T5||) by 
a stationary, axisymmetric spacetime with metric, g a p, of the form 

ds 2 = -e 2u{r) dt 2 + e 2X{r) dr 2 + r 2 d6 2 + r 2 sin 2 9dp 2 - 2uo{r)r 2 sin 2 6dtd V (166) 

(accurate to order Q) and with energy-momentum tensor 

T a/3 = (e + p)u a up + pg a/3 , (167) 

where e(r) is the total fluid energy density, p(r) is the fluid pressure and 

is the fluid 4-velocity to order Q - with t a = {dt) a and (p a = (d ip ) a , respectively the 
timelike and rotational Killing vectors of the spacetime. 

That the star is rotating slowly is the assumption that Q is small compared to 



the Kepler velocity, VLk ~ y M/R 3 , the angular velocity at which the star is unstable 
to mass shedding at its equator. In particular, we neglect all quantities of order fl 2 
or higher. To this order, the star retains its spherical shape, because the centrifugal 
deformation of its figure is an order Q 2 effect (Hartle pU[). 



In constructing such an equilibrium configuration, the equations (|116|) -( |T20| ) gov- 
erning a spherical star, 

p = p(e), (169) 



dp _ (e + p)(M + Anr 3 p) 
dr ~ r(r - 2M) 



(170) 
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dM , 



dr 

and 



47rr 2 e, (171) 



du 1 dp , 

= i_ ? (172) 

dr (e + p) dr ' 

with 

M ( r ) = I r (l_ e - 2 *), (173) 

are joined by an equation (Hartle [^|) that determines the new metric function u(r) 
in terms of the spherical metric functions z/(r) and A(r), 

e(u+X) d ( r 4 e -(»+V dcu \ 4 ( du + dX \oj - (174) 
r 4 \ dr J r \dr dr J 

where 

w(r) = fi-w. (175) 
Outside the star, Eq. (|174|) has the solution, 

2 7 

cu = fi-— (176) 

where J is the angular momentum of the spacetime. This new metric variable is 
a quantity of order Q and governs the dragging of inertial frames induced by the 
rotation of the star (Hartle |29|). Apart from the frame- dragging effect, however, 
the spacetime is unchanged from the spherical configuration. 

Since the equilibrium spacetime is stationary and axisymmetric, we may decom- 
pose our perturbations into modes of the formQ e K crt + m( p) We will use the lagrangian 
perturbation formalism reviewed in Sect. 1.2.2 and begin by expanding the displace- 
ment vector £ a and metric perturbation h a p in tensor spherical harmonics. 

The lagrangian displacement vector can be written 

1 00 r 1 

C Q = ~q E{ -Wi(r)Yrr a + Vtir^Yr (177) 

l=m 

-iU^Py^VpYrVjtVsAe' 



Sat 



where we have defined, 



pa = e {u+X) _ ^ya^ ^g) 



4 We will, again, always choose m > since the complex conjugate of an m < mode with 
frequency a is an m > mode with frequency —a. Note that a is the frequency measured by an 
inertial observer at infinity. 
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and the comoving frequency kQ = a + mil The exact form of this expression 
has been chosen for later convenience. In particular, we have chosen a gauge in 
which £ t = 0. Note also the choice of phase between the terms in ( |177| ) with polar 
parity (those with coefficients W\ and V/) and the terms with axial parity (those with 
coefficients U{). 

Working again in the Regge- Wheeler gauge, we express our metric perturbation 

as 



h 



in/ 



E 

l=m 



2vym 

Hi,i(r)Yf' 

symm 
symm 



H 2 ,i(r)e 2A Y, 



2\\rm 
I 



symm 
symm 



r 2 ^(r)Y" ; m 




ihi t i(r) sinOdgYf 


r 2 sin 2 0irj(r)Y l Tn 



Sat 



(179) 

Again, note the choice of phase between the polar-parity components (those with 
coefficients H j, Hu, H 2i i and K{) and the axial-parity components (those with 
coefficients h 0i i and hi t i). 

Based on our knowledge of the newtonian spectrum, we expect to find a class 
of hybrid modes whose spherical limit (Q —>■ 0) is a sum of the relativistic zero- 
frequency r- and g-modes found in the preceding section. We will, therefore, assume 
that our perturbation variables obey an ordering in powers of Q that reflects this 
spherical limit. 

W h V h U h H 1 , h h G , l ~ 0(1) 

H 0; i,H 2! i,Ki,h l! i,8e,5p,a ~ 0(D). 

In addition, we assume that the 0(1) quantities obey the 0(1) perturbation equations 
(1I53D , ([156D and for all /. 

Eq. ( |25| ) relates the Eulerian change in the 4- velocity to £ a and h a p, 

(181) 



;iso) 



5u a = q a p£ u f + \u a u p u'h M . 
The ordering ( 180|) implies that 5u a is given, to zeroth order in Q, by 

5u a = i/cfie-"^^, 



(182) 



and the spherical limit of this expression reveals that the mode (£ Q , h a p) is manifestly 
a sum of the axial and polar perturbations considered in Sect. 3.1. (Eqs. (|126|) , ( 127 ), 



1281) and @.) 
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In the newtonian calculation, it was the conservation of circulation in isentropic 
stars that brought about the rotational splitting of the zero-frequency modes at 
zeroth order in the angular velocity Q (see Sect. 2.2). The equation that enforces 
this conservation law is the curl of the perturbed Euler equation (|52|) and in general 
relativity it takes the form (Friedman [p^] ; see also Friedman and Ipser |S3| ) . 



= A£ v u 



or simply 



where 



= 2V[ a 



Auj afS = 0, 
' e + p 



a/3 



-u 



n 



ft 



(183) 

(184) 
(185) 



is the relativistic vorticity tensor. 
We begin by expressing Eq. ( |184j ), 



= Acj, 



aj3 



\ I e+ P 

A I up 



A 



t + p 



-u, 



n 



(186) 



in terms of £° and h a p. 

Making use of Eq. fl29|) we have 



A 



e + p 



-u, 



n 



e + p 



n 



1 

2' 



Au a - =-q aP Ag 



a/3 



IP 



u, 



where 



;i87) 



188) 



Au a = A(5f Q/3 M /3 ) = Ag a pu p + g a pAu p 

The ordering (|180|) implies that u a u^h a p and g a ^h a p vanish to zeroth order in Q 
since the only zeroth order metric components are h tr , h t g and h tip . Therefore, 



-u a u p Ag aP 
\q a ^9^ 



Au a = uVu 7 V^ 

Au a = h Q pu fS + u p Vp£ a + u p V a ip + u a u p u < Vpi 1 . 
From Eqs. ( p9|) and (|172|) and the relation, 

= -^Vpu + OiO), (193) 



(189) 

(190) 
(191) 
(192) 
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we obtain, 



u 



e + p 

v IP J 
Ve' 2A £ r 



v e 



to zeroth order in Q. We will also use the explicit form of u v from Eq. 



—v- 2 • 2/i 

u<p = e cur sin 0, 



and the components of Aw Q to zeroth order in Q, 



Au r 
Aug 

Au,„ 



/ 1 \ c 2u 

h tr + ZKfi£ r + nr 2 d r (— f J + — 8 r {r 2 Ue~ 2u ) £ 



e 

h te + inQ^g + fldgC,^ — 2u cot 9£ 



h tLp + iKft^tp + Q.d^w + 2u sin 9 cos 



+e u d r {r 2 ue- v ) sin 2 #e" 2A £ 



(194) 
(195) 

(196) 

(197) 
(198) 
(199) 



For reference, we explicitly write the components of inQ£ to zeroth order in Q, 



i/d2f = 0(D) 



1 



r 2 sin 9 



IK 



DC = y -WiY, m e ici mVt^ = y - , 3- [mV^Y, m + tf, sin #<9 e Y, m ] e i<T * 



r 2 sin' 9 



infl^t = 



wcfi& = ^ _ [VJ sin #<9 Y, m + m^F™] e iff * 



,2A 



= E — WiYre iat i/cfi^ = ^ z [mV^ m + tf, sin fld^" 1 ] e 



let 



(200) 

By making use of Eqs. ( |187| ) through ( |199| ) and the expressions fl200| ) and ( |179| ) 
for the components of ikQ£ and h a p, we may now write the spatial components of 



Au a p. We will use Eq. (|154j) to eliminate H\ i (for all I) from the resulting expressions 
and drop the "0" subscript on the metric functions hoi, writing h i = hi. 



Auj d u 



e + p 



n 



$yd 6 Au v - dp Au e - d e 



IP 
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e + p\ e v sin 6 



x 



W [1(1 + + C/i) - 27nwC/i] 

-2a)VJ [sin 00 fl y, m + /(/ + 1) cos OYf 1 } 



+ ^d r (r 2 ue- 2v ) Wt [sin Qd e Y{ n + 2 cos dY{"\ \e l 



Act 



(201) 



n J 
e + p 



d r (e-"Auo) - d e (e-"Au r ) 



n J kQ sin 

x m/tft<9 r [e-^(/i, + C/i)] - 2d r (Qe~ 2lJ Ui) cos 9 sin edgYf 



(202) 



i 



+ ^d r (r 2 uje~ 2v ) Ui [m 2 + 1(1 + l)(cos 2 6 - 1)] Y{ 



-2md r (ue- 2v Vi) cos OY^ + ^d r (r 2 ue- 2u ) VJ sin 00 fl 17 



[0 r (e- 2 ^) + e" 2 " (^g^ - I) e 2X Wt] sin^F^ 



Acu^ = 



e + p 



n 



e"| <9 V (e-'Au,.) - <9 r (e^Aw^) 



e + p\ e 



n J inQ 



(203) 
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-2<9 r (we- 2 "^) costfsinfl^Y,™ + ^<9 r (r^e" 2 ^) ViY^ 



+d r 



±& (r 2 ooe- 2u ) W, (cos 2 9 - l)Y, m 



-2md r {uie- 2u Ui) cosdY™ + nQd r [e" 2l/ (^ + 17,)] mnflfyl? 



+ ^|<9 r (r 2 we- 2y ) Ui smOdeY^je^ 

Note that the three spatial components of Au al 3 = are not independent. They 
are related by the identity 

V [a AuJ M = 0, (204) 
which, therefore, serves as a check on the above expressions. 



Let us write these three equations making use of the identities (p)7j)-(|58|) 



sinWr = ZQ, + i*£ a -(Z + 1)3,1^ 
cos BY? = Q^Y^ + QtY^ 



where Qi was defined in Eq. (|59| ) to be 

(I + m)(l — m) 



Qi 



(2Z-1)(2/ + 1) 



From Aw to = we have, 



+ [^d r [r 2 Qe- 2 ») Wi - 2Z^1 (/ + 2)Q l+l Y l 



rn 
l+l 



(205) 
(206) 



(207) 



(208) 



^d r (r 2 ooe- 2 n Wi + 2(1 + 1)^1 (/ - 1)3,1^ 
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From Au r e = we have, 

= E{ [~ 2 ^ i^~ 2v Ui) + ^d r (r 2 coe- 2 n U t ] lQ l+1 Q l+2 Y% 2 
UVld r {e- 2v Vi) - 2md r (cDe" 2 ^) 



+ 



Qi+tXi 



m 



+ 



mnttdr \e- 2u {hi + U x )\ + 2d r [yje^U^ ((/ + 1)Q 2 - IQ 2 1+1 ) 



+ - 



±d r (r 2 ue- 2 n U t [m 2 + 1(1 + 1) (Q 2 +1 + Q 2 - l) 



+ 



(I + l) K nd r (e-^Vi) + 2md r [yje- 2v V x ) 



+ ^d r (r 2 ue' 2 n V l + (l + l^Qe' 2 " (±^g±ri - J) e 2 ^" 
2<9 r (cue" 2 ^) + (r 2 cje~ 2 n Ut] (I + l)Q,_iQ,l^ 2 | 



QiYf: 



(209) 



From Au^ = we have, 
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o = E< 



Or 



1 d r {r 2 uje- 2v )Wi - 2W r {ue- 2u Vi) 



Q1+2Q1+1Y1+2 



+ 



Infldr [e- 2u (hi + Ui)\ - 2md r (ue- 2v Ui) + ^d r (r 2 uje~ 2v ) U t 



Qi+iYi 



m 
2+1 



+ 



mnVLdr [e-^Vx) + 2d r (ue^Vi) ((/ + l)Qf - lQ 2 l+l ) 



(Z + l)/cftd r [e- 2v {h x + tf,)] 



+2md r (cue" 2 ^) + ^^d r (r 2 uje~ 2v ) XJ X 



Qi Y i-i 



+ 



d r 



% (r 2 ue" 2u ) W t + 2(1 + l)d r (uje~ 2v Vx) 



m 
1-2 



(210) 

Let us rewrite the equations one last time using the orthogonality relation for 
spherical harmonics, 

J YfYrdQ = 5 lv 5 mm/ , (211) 

where dfl is the usual solid angle element on the unit 2- sphere. 
From Au dv = we have, for all allowed /, 



= [1(1 + l)/cfi(/i, + Ui) - 2mwU l ) 



(212) 



+(l + 1)Q, e —d r (r 2 iue- 2 -) - 2(1 - l)c^-i 



-iQi+i e —d r (r 2 Lue- 2 ») W l+1 + 2(1 + 2)uV l+1 



From Au r g = we have, for all allowed I, 



= (l-2)Q l _ 1 Q l 



(213) 
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+Qi 



(I - l) K Qd r (e-^V^) - 2md r (ue'^V^) 



+ 



mKVtdr \e- 2v {hi + U t )] + 2d r (we" 2 "!/',) ((/ + 1)Q 2 - lQ 2 l+l 



+±d r (r 2 iue- 2 ») U t [m 2 + 1(1 + 1) (Q 2 +1 + Q 2 - l) 



—Qi+i 



(I + 2)nVtd r (e~ 2v V l+l ) + 2md r (ue~ 2u V l+1 



+ ^d r (r 2 ue- 2 n V l+1 + + 2)Ktle- 2 » (^(0. - i) e 2 *W l+1 



+(l + 3)Q l+1 Q l+2 



2d r (uje- 2 »U l+2 ) + { ^-d r (r 2 u>e- 2 ») U l+2 



From Aui^r = we have, for all allowed I, 



= Qi-iQi 



-d r {r 2 ue' 2u ) W1-2 - 2(1 - 2)d r (cue^VU) 



(214) 



+Qi 



(I - l)Kttd r [e- 2 "(/i,_i + XJi-x, 



+ 



-2md r (iue^U^) + '^d r (r 2 cue- 2 <0 



mKm r (e-^Vi) + 2d r (ue~ 2v V x ) ((/ + 1)Q 2 - IQ 2 +1 ) 
+^-d r (r 2 coe- 2 n V x + d r [i0 r (rfce" 2 ") W,] (Q 2 +1 + Q 2 - l) 



+m^e-- (i^±rf - I) e* Wl 



-Qi+i 



(l + 2) K nd r [e- 2 »(h l+1 + U l+1 )] 



+2md r (ue- 2 »U l+1 ) + r -^d r (r 2 Qe~ 2v ) U l+1 



+Q1+1Q1+2 



Or 



-d r (r 2 ue~ 2 ») W l+2 \ + 2(1 + 3)d r (u}e~ 2v V l+2 
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It is instructive to consider the newtonian limit of these equations, 

u(r),u(r),X(r),hi(r) -> 0. (215) 



We have already seen that Eq. ( 157|) is the relativistic generalization of the new- 



tonian mass conservation equation (fjQ|) (or Eq. fl8"4|)), and that the other zeroth 
order perturbation equations, ( |154| ) and ( |1 56| ) , simply vanish in the newtonian limit. 



Similarly, one can readily observe that the conservation of circulation equations have 
as their newtonian limit the corresponding equations presented in Sect. 2.2, 

Eq. (HI) - Eq. ©, 

Eq. O - Eq. 

Eq. (|T|) -> Eq. (§§ 



(and similarly for the other forms of these equations). 

This correspondence leads us to expect the same structure for the relativistic 
modes as was found in the newtonian case: we expect to find a discrete set of 
axial- and polar-led hybrid modes with opposite behaviour under parity. Further, we 
expect a one-to-one correspondence between these relativistic hybrid modes and the 
newtonian modes, to which the relativistic hybrids should approach in the newtonian 
limit. 



3.3 Character of the Perturbation Modes 

In deriving the components of the curl of the perturbed Euler equation in newtonian 
gravity (|65|)-(|6T|), we required no assumptions about the ordering of the perturbation 
variables (dp, 5v a ) in powers of the angular velocity Q. Thus, our theorem concerning 
the character of the newtonian modes, Thm. (P applied to any discrete normal mode 
of a uniformly rotating barotropic star with arbitrary angular velocity. 

We conjecture that the perturbations of relativistic stars obey the same principle: 
If (£ a , h a p) with £ a ^ is a discrete normal mode of a uniformly rotating stellar 
model obeying a one-parameter equation of state, then the decomposition of the 
mode into spherical harmonics Y™ has I = m as the lowest contributing value of I, 
when m^O; and has or 1 as the lowest contributing value of I, when m = 0. 

However, in deriving the curl of the perturbed Euler equation for relativistic 
models we have imposed assumptions that restrict its generality. We have derived 
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Eqs. (|212 )- (pT4|) in a form that requires a slowly rotating equilibrium model, assumes 
the ordering (|180|) and neglects terms of order Q 2 and higher. Under these more 



restrictive assumptions, the following theorem holds. 

Theorem 2 Let (g a p(Q), T a/ 3(Q)) be a family of stationary, axisymmetric space- 
times describing a sequence of stellar models in uniform rotation with angular ve- 
locity fl and obeying a one-parameter equation of state, where (g Q/ 3(0), T Q/ g(0)) is 
a static spherically symmetric spacetime describing the non-rotating model. Let 
(^ a (Q),h a 0(Q)) with £ a ^ be a family of discrete normal modes of these space- 
times obeying the same one-parameter equation of state, where (£ a (0), /i a /?(0)) is a 
stationary non-radial perturbation of the static spherical model. Let (£ a (f2o), h a p(Qo) 
be a member of this family with Qq -C Qk, the angular velocity of a particle in orbit 
at the star's equator. Then the decomposition of (£ a (f2 ), h a p(Qo) into spherical har- 
monics Y™ (i.e., into (l,m) representations of the rotation group about its center of 
mass) has I = m as the lowest contributing value of I, when m ^ 0; and I = 1 as the 
lowest contributing value of I, when m = 0. 

We designate a non-axisymmetric mode with parity (— l) m+1 an "axial-led hy- 
brid" if £° and h a p receive contributions only from 

axial terms with I = m, m + 2, m + 4, ... and 
polar terms with / = m + 1, m + 3, m + 5, .... 

Similarly, we designate a non-axisymmetric mode with parity (— l) m a "polar-led 
hybrid" if £ a and h a p receive contributions only from 

polar terms with / = m, m + 2, m + 4, ... and 
axial terms with / = m + l, m + 3, m + 5, .... 

For the case m = 0, we designate axisymmetric modes with parity^] +1 "axial- led 
hybrids" if £ a and h a p receive contributions only from 

axial terms with / = 1, 3, 5, . . . and 
polar terms with I — 2, 4, 6, . . ., 



5 The family of modes for which £ Q and h a p receive contributions only from polar terms with 
Z = 0,2,4,... and axial terms with Z = 1,3,5,... would also have parity +1 and could be designated 
"polar-led hybrids." However, these modes require a more general theorem to establish their 
character. 
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and we designate axisymmetric modes with parity —1 "polar- led hybrids" if £ Q and 
h a /3 receive contributions only from 

polar terms with I = 1, 3, 5, . . . and 
axial terms with Z = 2, 4, 6, . . .. 

We prove the theorem separately for each parity class in Appendix C. 



3.4 Boundary Conditions and Explicit Solutions 



A physically reasonable solution (£ a , h a p) to the perturbation equations ( |154| ), (|156| ) 



(|157|) and ( |212| )- (|2T4|) , must be regular everywhere in the spacetime. Of course, the 



fluid variables Wi(r), V/(r) and Ui(r) (for all I) have support only inside the star, 
r G [0, R\. The metric functions Hij(r) will also have support only inside the star 
(for all I), since they are directly proportional to Wi(r) by Eq. ( |154j ). The metric 
functions hi(r), on the other hand, satisfy a nontrivial differential equation, (|156|) , 



in the exterior spacetime and will, therefore, have support on the whole domain 
r G [0, oo]. Let us now consider the boundary and matching conditions that our 
solutions must satisfy. 

At the surface of the star, r = R, the perturbed pressure, Ap, must vanish. (This 
is how one defines the surface of the perturbed star.) The lagrangian change in the 
pressure is given by Eq. fl2D|), 

Ap=- 1 - 1 pq^Ag afi . (216) 



Making use of Eq. ( |194j ) and the equilibrium equations ( |170| ) and ( |172| ), we find that 
at r = R 

° = ^= R^R M °2M ) ^ W '^ Y " eM < 217) 
where M = M(R) is the gravitational mass of the equilibrium star and satisfies 
2M < R. 

For the equations of state we consider^ the energy density e(r) either goes to a 
constant or vanishes at the surface of the star in the manner, 

k 

r) 

This restriction can be dropped if the boundary condition Ap(r = R) = is replaced by 



(■(/•) - ( i -4^1 



Ah(r = R) = 0, with h the comoving enthalpy defined below in Sect. 4.1.1, Eq. (299). 
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(for some constant k). If e(R) ^ 0, then Eq. (|217|) requires that Wi(R) = for all 
I. Otherwise, one finds from Eq. ( |157| ) that Vi(r) will diverge at the surface unless 
Wi(R) = 0. Thus, the boundary condition, 



Wi(R) = (all I) 



(218) 



must be satisfied at the surface of the star. (By Eq. (|154|) , this also implies that 
H\^i{r) vanishes at the surface of the star). 

In the exterior vacuum spacetime, r > R, we have only to satisfy the single 
equation[| ( |156| ) for all I, which becomes 



r " 

\ + 



(2-P-l) c2X 



or 



1 



2M „ 
— ) h i 



1(1 + 1) 



4Mn 



0. 



hi = Q, 



(219) 
(220) 



where we have used e~ 2A = (1 — 2Mq/t) for r > R. 

Since this exterior equation does not couple hi(r) having different values of I, we 
can find its solution explicitly. The solution that is regular at spatial infinity can be 
written 



h t (r) = J2hs 



s=0 



R\ l+S 
r J 



(221) 



If we substitute this series expansion into Eq. (|220|) , we find the following recursion 
relation for the expansion coefficients, 



(2M Q \ (l + s-2)(l + s + l) 
V R 



(222) 



s(2l + s + l) 

with hifi an arbitrary normalization constant. We, therefore, have the full solution 
to zeroth order in Q of the perturbation equations in the exterior spacetime. 

This exterior solution must be matched at the surface of the star to the interior 
solution for hi(r). One requires that the solutions be continuous at the surface, 



Jim [hi(R -e)- h(R + e)\ = 0, 



(223) 



for all /, and that the Wronskian of the interior and exterior solutions vanish at 
r = R, i.e. that 



lim \hi(R - e)h{(R + e)-h{(R- e)hi(R + e)] = 0, 

6— >0 



(224) 



7 This equation was first written down by Regge and Wheeler |65j in the context of Schwarzschild 
perturbations. 
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for all I. 

Thus, in solving the perturbation equations to zeroth order in Q we need only 
work in the interior of the star (as in the newtonian case). Inside the star, the 
perturbation (£, a ,h a/3 ) must satisfy the full set of coupled equations ( |154| ), (|156|) , 



(|157|) and (|212| )-( [2T4] ) for all /, subject to the boundary and matching conditions 
( |2T8|) , (|223|) and Q22§. 



Finally, we note that since we are working in linearized perturbation theory there 
is a scale invariance to the equations. If (£ a , h a p) is a solution to the perturbation 
equations then (K^ a , Kh a p) is also a solution, for constant K. Thus, in order to 
find a particular mode of oscillation we must impose some normalization condition 
in addition to the boundary and matching conditions just discussed. We choose the 
condition that 

U m (r = R) = 1 for axial- hybrids, or that (225) 
U m+ i(r = R) = 1 for polar-hybrids. 

3.4.1 The purely axial solutions 

In Sect. 3.3 we saw that if a mode of a slowly rotating star has a stationary non-radial 
perturbation as its spherical limit, then it is generically a hybrid mode with mixed 
axial and polar angular behaviour. However, we have also seen in Sect. 2.3.1 that 
newtonian stars retain a vestigial set of purely axial modes - the so-called "classical 
r-modes" - whose angular behaviour is a purely axial harmonic having / = m. Let 
us now address the question of whether or not such r-mode solutions exist in the 
relativistic models. 

For relativistic stars, Kojima has recently derived an equation governing 



purely axial perturbations to lowest order in the star's angular velocity. Based on 
this equation, he has argued for the existence of a continuous spectrum of modes, 
and his argument has been made precise in a recent paper of Beyer and Kokkotas 
j/J. Beyer and Kokkotas, however, also point out that the continuous spectrum they 
find may be an artifact of the vanishing of the imaginary part of the frequency in 
the slow rotation limit. (Or, more broadly, it may be an artifact of the slow rotation 
approximation. ) 

In addition, Kojima and Hosonuma have studied the mixing of axial and polar 
perturbations to order Q 2 in rotating relativistic stars, again finding a continuous 
mode spectrum. Their calculation uses the Cowling approximation (which ignores 
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all metric perturbations) and assumes an ordering of the perturbation variables in 
powers of Q which forbids the mixing of axial and polar terms at zeroth order. 

In contrast to these results, we do not find a continuous spectrum of purely 
axial modes for isentropic stars, but rather, a discrete spectrum of axial-led hybrids. 
Indeed, we show below that there are only two purely axial modes (and their complex 
conjugates) in isentropic rotating relativistic stars. Both are discrete stationary 
modes having spherical harmonic index 1 = 1. (One mode has m = and the other 
has m = 1.) These modes are the generalization to slowly rotating stars of the I = 1 
axial modes discussed by Campolattaro and Thorne []TTJ] and the generalization to 
relativistic stars of the I = 1 axial modes discussed in Sect. 2.3.1. In particular, 
none of the newtonian r-modes having I = m > 2 retain their purely axial character 
in isentropic relativistic stars. 

As in Sect 2.3.1, let us write down the equations governing an axial mode be- 
longing to a pure spherical harmonic of index I. In other words, let us assume that 
hi(r) and Ui(r) (for some particular value of I) are the only non- vanishing coefficients 
in the spherical harmonic expansions of the lagrangian displacement ( p.77| ) and the 
perturbed metric ( |1 79| ) . The set of equations that have to be satisfied are the ze- 
roth order (spherical) equations ( |154j ), ( |156| ) and ( 157 ); the order Q conservation of 
circulation equations ( 212 )-( |214| ) and the matching conditions at the surface of the 
star (12231) and 



With hi(r) and E/j(r) the only non-vanishing perturbation variables, Eqs. ( 154 ) 
and (|157 ) vanish identically, while Eq. ( |156j ) remains unchanged, 



h'; - (i/ + \')h[ + 

Eq. (|212j) becomes, 



(2-Z 2 



2 2 



hi 



r 



and Eq. (PI3|) with I 

= IQi+xQi+2 



= [1(1 + l)/cfi(/ii + Ui) - 2miuUi] , 

1 + 2, I — > I and I —> I — 2 gives the equations 8 

(l + l 



-2d r (. 



ue 2u Ui 



-d r ( 



2- —2u 

r uoe 



Ui 



3 Alternatively, one can get these equations from the coefficients of Y™ and YJ™ 2 in Eq. 



(226) 
(227) 

(228) 
0£ 
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(229) 



+2d r (ue-^Ut) ((/ + 1)Q? - lQf +1 ) 



+ z(z + i)(^ +1 + g?-i 



(/ + i)Q,_ig, 



2<9 r (uoe^Ui) + -dr ( 



r 2 u;e 



2^^ 



(230) 



respectively. Recall that we need only work with two of the three equations ( 212 )- 
(|214|) since they are linearly dependent as a result of Eq. fl204j ) . 

By Thm. (0), we know that a non-axisymmetric (m > 0) mode must have I = m 
as its lowest value of / and that an axisymmetric (m = 0) mode must have I = 1 
as its lowest value of /. (Hence, in the present context of pure a spherical harmonic 
these are also the only allowed values of /.) We consider each case separately 



The case m = and 1 = 1. 



We seek a solution to Eqs. (|226|) -( [2~3"0| ) with m = and 1 = 1. The fluid perturbation 
turns out to have a particularly simple form in this case. From the definition of Qi, 
Eq. pi, we find that 



QU = 0, Q 2 = l - and Qhi = ^- 



(231) 



These imply that Eq. (|230|) is trivially satisfied, while Eqs. (|228D and ([229]) both 
become, 

= d r (ue-^Ui) - ^d r (r 2 uoe- 2v ) U t (232) 

with solution 

Ut(r) = K ir 2 } (233) 

for some constant K\. 



The exterior solution (|221| ) also takes on a particularly simple form in the present 



case. When 1 = 1, the recursion relation ( |222j ) terminates and the exterior solution 
simply becomes 



hi{r) 



(234) 



for some constant K2. 



ss 



Eq. ( |227|) can be satisfied for m = if either k = or hi = — t/j = —K\r 2 . In 
the latter case, however, the matching conditions ( |223| ) and at the surface of 

the star would require that K2 = —K\R 3 = 0. Thus, a non-trivial solution will exist 
if and only if 

k = 0. (235) 

(Since m — 0, this implies that the frequency <r also vanishes.) 

Finally, we consider Eq. ( |226| ). Defining f(r) = hi(r)/r 2 , it is not difficult to 
show that 



3 *. 



2 2 

V+ao + 4 



(236) 



4_-(i/+A) r/ 



r e 



/ 4 



/') + A')/. 



But this is simply the operator appearing in Hartle's [29| equation, Eq. (|1 74 ) ! 
Indeed, defining^ 



n 
j 

oj(r) 



-i-Y 

2 Utt; 
3 y 

Cl — Cj(r 



K 2 
f(r) 



3_\ 2 hi(r) 
Air 



(237) 

(238) 

(239) 
(240) 



Eq. (|226|) becomes precisely Hartle's equation - here governing a change Q in the 
uniform angular velocity Q of the star and the associated changes J and u in the 
angular momentum J of the spacetime and the frame-dragging metric variable u, 
respectively. 

Just as in the newtonian case (Sect 2.3.1), we find that the only allowed purely 
axial mode with m = is a stationary perturbation in the fluid velocity of the form, 



5u c 



(241) 



Associated with this fluid perturbation is a change in the frame-dragging metric 

variable uo — > uo + u given by a solution to Eq. ( |174|) for u = Q — Cj. Observe that 

9 The exact form of these expressions assumes a standard definition of the spherical harmonics 
Y™. (See, e.g. Jackson |3(|, p. 99.) In particular, we are using Y° = ^S/An cos#. 
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this change, Cl, in the angular velocity of the star is uniform - a fact which follows 
from the radial dependence of the perturbed velocity field, Eq. ( |233| ). In a star that 



is already rotating uniformly, a perturbation inducing differential rotation would vi- 
olate conservation of circulation. Thus, the equation that enforces conservation of 
circulation - the curl of the perturbed Euler equation - with components ( |212|) - (|2T4]) , 
restricts the radial behaviour of 5u a . In the non-rotating star, however, such a re- 
striction does not apply (See Sect. 3.1.4 and Campolattaro and Thorne [JT/TJ ) . Since 



Eqs. ( |212j )-( E?T4]) are of order Q, they vanish in the spherical limit and impose no 
restriction on the form of Ui(r). In freely specifying Ui(r), in this case, one is spec- 
ifying the form of differential rotation about the 2-axis. The corresponding metric 
perturbation (the frame dragging term) is then determined by Hartle's equation in 
the form of Eq. 



The case I — m > 0. 

We now seek solutions to Eqs. (|226|)-fl2"3U|) with I = m. The fluid perturbation, 



again, turns out to have a particularly simple form. In this case, Eq. (|227|) becomes, 

(m + l)K,Q(h m + U m ) = 2uU m . (242) 
From the definition of Qi, Eq. (^), we find that 

Ql = 0, and Ql +1 = J + - . (243) 



These, again, imply that Eq. ( |230p is trivially satisfied, while Eqs. ( |228| ) and ( |229| ) 
both become, 

= 2d r (ue- 2 »U m ) - ^^d r (r 2 uje- 2 ») U m (244) 

where we have used Eq. ( |242| ) to substitute for nQ(h m + U m ) in Eq. ( [229| ). 
This last equation has solution 

U m {r) = K 1 r m+1 u^ m ~ 1) e' (m ^ 1)u (245) 

for some constant K\. Then, Eq. (|242|) implies, 



h m (r) = K 1 



2lo 

1 



(m + 



rm+1 _|( m _i) e _ (m _i) I , ^ 246 ^ 
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Finally, consider Eq. ( |226| ). Letting / = uj^e and substituting for U m (r) and 
h m (r) the expressions ( |245|) and ( |246| ) , respectively, one can show that, 



(m — 1)< 



+ 



(2-^- TO ) eK _g (i/ + V) _2 



Auo' 



(m + l)Ktt \f r 



h r 



(247) 



+ 



2uj 



(m + 1)k£1 



f" , x ff'Y 2(m+ 1) ( f 



- u/ 



A') 



,2A 



If m = 1 this equation is obviously satisfied. However, with some effort, one can 
also show that it is satisfied only if m — 1. For m > 1, therefore, the system of Eqs. 
(|226|) - (|230|) is overdetermined and no solutions exist. The purely axial modes with 
I = m > 2 do not exist in isentropic relativistic stars. 

This result contradicts the claims by Kojima [|38| and Kojima and Hosonuma 
P§| that a continuous spectrum of purely axial modes exists in isentropic relativistic 
stars. 

Kojima |38 bases his conclusion on Eqs. (|226|) and (|227|) (which he derives us- 
ing slightly different notation) and he does not distinguish between the isentropic 
and non-isentropic cases. In non-isentropic stars, Eqs. ( |226| ) and ( |227| ) are, in- 
deed, the only equations^] governing axial perturbations to order Q. Kojima P%| 
derives the striking result that these equations can be combined into a single "mas- 
ter equation" whose highest derivative term has a frequency-dependent coefficient. 
The vanishing of this coefficient at a particular radius when the frequency is real 
gives rise to the continuous spectrum. However, because the frequency is complex 
(with imaginary part higher order in Q) we conjecture that a higher-order calculation 
for non-isentropic models will find only a discrete set of modes - those modes that 



10 In isentropic stars Eqs. ( 228 )-( [230|) have the form shown and involve only hi and Ui. In non- 
isentropic stars, however, these equations are modified in such a way that they give a coupling to 
order f2 polar variables such as the perturbed pressure and density. 
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generalize to relativistic stars the newtonian r-modes first studied by Papalouizou 
and Pringle |60| (see also Provost et al. [Q, Saio ||| and Smeyers and Martens 

For the case of isentropic stars, on the other hand, it is no longer true that Eqs. 
(|226|) and ( |227|) are the only equations to be satisfied. As we have just seen, these 
are joined by Eqs. ( [228| )-( |230|) and comprise an overdetermined system for modes 
with m > 2. The result is not a continuous spectrum of purely axial modes, but no 
such modes at all! 

Kojima and Hosonuma perform a higher order calculation and they do distin- 
guish between what we are calling the isentropic and non-isentropic cases. (They use 
the terms "barotropic" and "non-barotropic".) However, they also make the simpli- 
fying assumption that the metric perturbation is small and can be ignored altogether 
(the Cowling approximation). Based on this assumption they, again, claim to find a 
continuous spectrum of pure r-modes. However, we have seen in the isentropic case 
that the full perturbation equations (which include contributions from the perturbed 
metric) forbid the existence of such r-modes. Again, in the non-isentropic case, we 
expect that a higher order calculation that includes the metric perturbation will find 
only a discrete r-mode spectrum. 

Finally, let us return to the solution (|24 5|) - (|246|) with I = m = 1, 



U m {r) 



K ir 2 
( & 



1 \r l 



(248) 
(249) 



This solution must satisfy the matching conditions 
solution with I = 1, Eq. (1234]) , 



h m (r) = — 



K 2 _ 

r 



21) and (E24|) to the exterior 



(250) 



Eq. 



gives, 



= K t R 3 



'u(R) 



- 1 



while ( |224|) gives, 

= KiR 3 
Adding these equations, we find 



'a) (it!) 
k£1 



1 +R 



uo'{R) 



+ K 2 



(251) 



(252) 



= u{R) + -Rw'(R) - Kfi 



(253) 
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and using Eq. ( |176|) to substitute for u>(R) and u'(R) we find, 

= Q - (254) 

or 

k = 1. (255) 

Thus, the purely axial solution with I = m = 1 has a co-rotating frequency equal 
to the angular velocity of the star. Since kQ = o + mfl, this implies that a = 0, 
i.e., that the mode is stationary as seen by an inertial observer. This mode is the 
generalization to slowly rotating stars of the I = m = 1 mode found in spherical stars 
by Campolattaro and Thorne |IT| , and the generalization to relativistic stars of the 



I = m = 1 mode found in Sect 2.3.1. It represents uniform rotation about an axis 
perpendicular to the rotational axis of the star. We note, again, that only uniform 
rotation is an acceptable perturbation, as was the case with the m — 0, I — 1 mode. 

In newtonian isentropic stars there remained a large set of purely axial modes with 
I = m\ the I = m = 2 mode being the one expected to dominate the gravitational 
wave-driven spin-down of a hot, young neutron star. In relativistic stars, however, 
we see that all such pure r-modes with I = m > 2 are forbidden by the perturbation 
equations, and instead must be replaced by axial-led hybrids. Let us then turn to 
the problem of finding these important hybrid modes. 

3.4.2 Relativistic corrections to the "classical" r-modes in 
uniform density stars 

Before turning to the general problem of numerically solving for the hybrid modes 
of fully relativistic stars, let us look directly for the post-newtonian corrections to 
the I = m newtonian r-modes. The equilibrium structure of a slowly rotating star 
with uniform density is particularly simple (see Chandrasekhar and Miller |IJ|) and 
lends itself readily to such a post-newtonian analysis. 
For a spherically symmetric star with constant density, 

3M 

<r) = ^, (256) 
the equilibrium equations ( |169| )-( [T73| ) have the well-known exact solution inside the 
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star (r < R), 



p(r) 

M(r) 
e M»0 

-2X(r) 



(1 



2M \ 2 
R 



2Mr 



(tf 



m " [ r) 



1 _ 2Mo 
1 R 



3 
2 

1 - 



\ 2M /r\' 
R \RJ 

2M ( r \ 2 



1 _ 2Mo ^ 2 



1 A 2M \2 

2 V i2 



(257) 

(258) 
(259) 
(260) 



where M is the gravitational mass of the star and R is its radius. (See, e.g., Wald 



H Ch. 6.) 

To find the equilibrium solution corresponding to the slowly rotating star, we 
must also solve Hartle's ]2?| equation ( |174| ), 

= r 2 u" + [4 - r{y' + \')}ruj' - 4r(i/ + A')w 

where we may use the spherical solution to write 

r{y' + A') = 4vrr 2 (e + p)e 2X 



(261) 



(262) 
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-(l-3£)*} 



To simplify the problem, we expand our equilibrium solution in powers of (2M /R) 
and work only to linear ordeij^]. We will need the expressions, 



r{v 2 [r) [-R 



o 



2M 



and 



e~ 2v = 1 + 



3 1 (r x 2 
2 ~ 2 \R 



2M 



+ 



/ 2M 
V R 



(263) 
(264) 



Since we are also working to linear order in the star's angular velocity, we may 
set O = 1 without loss of generality. We write, 



(265) 



n This expansion will give us the first post-newtonian (1PN) corrections to the I = m newtonian 
r-modcs. 
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and solve Eq. ([261]) subject to the following boundary condition (Hartle p9[ ) at the 
surface of the star, 



1 = 

To order (2M /R) the solution is, 



u) r 



1-1 



CiJ + -it 



3r 2 \ /2M 



5# 2 / V R 



r=R 



o 



/2Mc 
V R 



(266) 



(267) 



With this explicit equilibrium solution in hand, we now consider the perturbation 
equations. We are required to solve Eqs. ( |154j ), (|156|) , ( |1571 ) and (|212|) - (|214|) subject 
to the boundary, matching and normalization conditions (|218|) , (|223|) , ( |224|) and 
( |225|) . We seek the post-newtonian corrections to the / = m newtonian r- modes 
discussed in Sect. 2.3.1, Eqs. ( ]78| ) and (P|), 

2 



(to + 1) 

r \ m+l 

R, 



(268) 



(269) 



Therefore, let us make the following ansatz for our perturbed solution inside the star, 



U m (r) 
h m {r) 
W m+1 (r) 
V m+1 (r) 
U rn+2 (r) 



m + l) 

m+l 



1 + Kl 



2M 



O 



y. \ m+l 

R, 



1 + u mfl 1 



R 2 



/'//Ml : /'///.J. ( -j^J 



2M 
2M 
2M 
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/2Mq 
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2M n ^ 2 



i? 



\ m+l 



R 



m+l 



u m+2,0 



R 2 

Vmfi + v m,l 
m + 3 /2M 



R 



O 



2ilV 2 



R 



O 



2M 



2Mn N 2 



/2Mq 



R 



(270) 
(271) 
(272) 
(273) 
(274) 
(275) 



where we have chosen the form of U m (r) so as to automatically satisfy the normaliza- 
tion condition ( 225 ) and we have chosen the form of W m+ i(r) so as to automatically 
satisfy the boundary condition ( |218p . Note that we have assumed that hi, Vy, Wy 
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and JJyi are of order (2Mq/R) 2 or higher for alH > m, I' > m + 1 and I" > m + 2. 
We will justify this ansatz by finding a self-consistent solution to the perturbation 
equations. 



Observe that the exterior solution (|22 1|) for h m (r) already has a natural expansion 



in powers of (2M /R) as a result of the recursion relation (|222|) 



h m (r) = h mfl [ 7 ) (—j + o^—) . 

The normalization constant, h m fi, is determined by the matching condition 



(276) 



<>"ni,0 

+ h 



m,l 



while ( [224D imposes the following condition on the interior solution, 



= h mfi \ -m(h mfl + h m>1 ) - [(m + l)h mfl + (m + 3)/i mi i] 



(277) 



(27£ 



or, 



= (2m + l)/i mj0 + (2m + 3)/i mjl (279) 

Turning now to the perturbation equations, we will write Eqs. ( |156[ ), ( |157| ), ( |212| ) 
and (|213|) to order (2Mq/R). Eq. ( [154] ) merely expresses Hij,{r) in terms of Wi(r), 
and we need not work with Eq. (|214j) since (|212| )-( pT4l ) are related by Eq. (|204|) . 
Hence, a complete set of perturbation equations, accurate to first order in (2M /R), 
is as follows. 

Eq. ([157D with I = m + 1 is, 



= (m + l)(m + 2)V m+1 - (rW m+1 )' 



Eq. (|156|) with I = m is, 



= r 2 h" m - m(m + l)h m - 4r(i/ + X)U„ 



(280) 



(281) 



Eq. (|2~I2D with / = m is, 



= m(m + l)nh m + m [(m + 1)k — 2ui\ U n 



(282) 



-2mQ m+ i [W m+ i + (m + 2)V m+1 ] 



Eq. ( gig) with Z = m + 2 is, 
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0= [(m + 2)(m + 3)/t-2m][Z m+2 



(283) 



+2(m + 3)Q m+2 [W m+ i - (m + l)V m+1 ] 
Eq. (pTBj) with Z = m is, 



= mKrh' m + m/te^r (e" 2iy Z7 m )' - j^™^ (uje- 2u U, " ' 



m(m + 2)e 2 ' 
(2m + 3)r 



(r 2 ue- 2u )'u m + 2(m + 3)Q m+1 Q m+2 [rZ7^ +2 + (m + 2) (7, 



y m+2 



-Qm+i {[(m + 2)« + 2m]rV' i+1 + 2m(m + 2)V m+1 - (m + 2)«tF m+ i} 

(284) 

where we have used the fact that = and Q m+1 = l/(2m + 3), and we have set 
f2 = 1. One can readily verify that all other non-trivial equations are satisfied by a 
solution to these; for example, Eq. ( [213] ) with Z = m + 2. 

We now substitute for the equilibrium quantities (z/ + A'), e~ 2v and Co using Eqs. 
( |263| ), ( |264| ) and ( [26 7| ), respectively. We also substitute for the perturbation variables 
k, U m , h m , W m+ i, V m+ i and Z7 m+2 using our ansatz, Eqs. ( |270| ) to (|275|) . Collecting 
powers of (2M /R), one finds that the zeroth order terms vanish identically as a 
consequence of the newtonian solution. At order (2M Q /R) we find the following set 
of equations. 

Eq. ( |280p becomes, 

m+1 







(m + 2 



(m + l)f m +i,o - w m+lfi 



(285) 



+ 



(m + l)(m + 2)u m+ i j i + (m + 4)u? m+ i j0 



7? 



m+3 



Eq. ( |281|) becomes, 



2(2m + 3)Zi mi i - 6 



R 



m+3 



(286) 
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Eq. (|282|) becomes, 
= 2m\ h mfl + Ki + 1- Q m +i 



w m +i,o + {m + 2)v m+h0 



R 



m+l 



(287) 



+2ml h m>1 - - + Q m+1 



Wm+1,0 - (m + 2)v m+1; i 



R 



m+3 



Eq. (|283|) becomes, 



= 2(m + 3)Q, 



m+2 



%+i.fl - (m + l)u m +i,o 



R 



m+l 



(4(2m + 3) n/ „. 

+{ — f — ^-Um+2,0 ~ 2(m + 3)Q m+2 

[ (m + l) 



Eq. ( |284| ) becomes, 



= < 2mh m + 2mKi + 2m 



R 



m+3 



(288) 



-2Q 



m+l 



[(m + 2) + m(2m + 3)]t> m+ i i0 



m 



(m + l 



■Wra+1,0 



1? 



m+l 



f 2m(m + 3) 4m(m + 2) m(m + 3) 

+ "( — 7 ! ~ — h m< i — 7— — ; 7777^— — — —^ U rn,0 ' 



m + 1 



(m + l)(2m + 3) 



m + 1 



m(m + 3) 2m(m + 2) 



5(2m + 3) 5(2m + 3) 



2(m + 3) (2m + 5)Q m+ iQ m+2 u m+ 2,o 



2Q 



m+l 



(m + 1, 



[(m + 2)(m + 3) + m(m + l)(2m + 5)]v m+1>1 

!/ r \ m+3 

(289) 

We must now solve these algebraic equations for the eight constants K\, M m> o, 
h m fi, h m ,ii w m +i,o, fm+1,1 and u m+ 2,o defined by our ansatz. Observe that of 

the four equations obtained from the coefficients of (r/R) m+1 , only two are linearly 
independent. (This is an example of the linear dependence in the perturbation 
equations expanded about r = discussed in detail in Appendix D.) Thus, we have 
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seven independent equations together with our matching condition ( |279|) for the eight 
unknown quantities. We find the following solution. 



K 



(m + 1) 



I _ 4(m-l)(2m + ll) (2M ^ + Q (2M„ 



5(2m + l)(2m + 5) V R 



R 



U m {r) 



m+l 



R 2 \ R 



R 



21 



(290) 
(291) 



h m {r) 



\ m+l 



R 



+ 



(2m + 1) (2m + 3) \RJ 



2M \ +(J f2M 



R 



R 



(292) 



W m+ x(r) = (m + l)(m + 2)K( 



m+l 



R 2 \ R 



2M n N 2 



(293) 



V m+1 (r) = K 



r\ m + 1 . Jr\ 2 /2iW Q \ /2M \ 

_) ( m + 2) _ (ro + 4)(-) (^2) +0 (^£) 



^ +0 f^) (294 ) 



^m+2(r) = 

where we have defined 



(2m + 3) \RJ V i? 
6(m - l)Q m+ i 



R 



K = 



and where 



KQ 



m+l 



24m(m + 2)(2m + 3) 



5(m + 2)(2m + 5) 



48(m+ l) 4 (m + 3) 2 



+ (2m + 3) 2 (2m + 5) 



m(m + 2) 2 - 4£ 



(295) 



(296) 



(297) 



Since our solution satisfies the full perturbation equations to order (2M /R), our 
ansatz was self-consistent. Thus, we have explicitly found the first post-newtonian 
corrections to the I = m newtonian r-modes of uniform density stars. 

The solution reveals the expected mixing of axial and polar terms in the spherical 
harmonic expansion of 5u a . All of the newtonian r-modes with m > 2 pick up 
both axial and polar corrections of order (2Mq/R). (When m = 1, the constant 
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K vanishes and we recover the purely axial solution with I — m — 1 described in 
Sect. 3.4.1.) In addition, we see that the newtonian r-mode frequency also picks 
up a small relativistic correction. To order (2M /R) the frequency decreases, but it 
is not clear whether this represents a tendency to stabilize or destabilize the mode. 
Finally, we note that the metric perturbation (whose radial behaviour is determined 
by the function h m ) is of the same order as the post-newtonian corrections to the 
fluid perturbation. Thus, there is no justification for the Cowling approximation in 
constructing the hybrid mode solutions. 

This analytic solution provides some useful insights into the effect of general 
relativity on the rotational modes of neutron stars. However, it is limited to uniform 
density stars, it is accurate only to first post-newtonian order and it applies only to 
those modes that are purely axial in the newtonian star. Thus, let us now examine 
the numerical methods that will allow us, in principle, to relax all three of these 
restrictions. 
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Chapter 4 

Relativistic Stars: Numerical 
Results 



4.1 Method of Solution 

We now consider the problem of solving numerically for the hybrid modes of slowly 
rotating relativistic stars. We follow the same approach used for the newtonian 
problem (Sect. 2.4). That is, we expand all quantities in regular power series about 
the center and surface of the star and solve an algebraic system for the coefficients 
of these expansions. 

In the newtonian case, the equilibrium solution did not play a large role in the 
perturbation equations, entering only into the perturbed mass conservation equation 
(Q). For the relativistic problem, however, equilibrium variables appear in all of the 
perturbation equations. Thus, we begin with a discussion of our numerical solution 
to the equilibrium equations ( |169|) - (|174|) . 



4.1.1 Numerical solution of the equilibrium equations 

As discussed in Sect. 3.2, we require our equilibrium solution to be that of a slowly 
rotating perfect fluid obeying a barotropic (one-parameter) equation of state, 

p = p{e). (298) 

For such models, it is convenient to define a comoving enthalpy, 

h{p) = f P . y - (299) 
v ' Jo e(p')+p> K ' 
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and to re-express the equilibrium equations ( |169|) - (|T74D such that h is the indepen- 
dent integration variable rather than r (Lindblom |46|]). We write these modified 
equations below. 

To integrate the equations in their usual form (for a given equation of state), one 
begins by fixing a quantity at the center of the star such as the central energy density 
e c or the central pressure p c . One then integrates Eqs. ( |169| )-( |TT4D from r = out 
to the surface of the star - i.e., to the radius, r = R, at which the pressure drops to 
zero. 

With the equations re-expressed so that h is the independent variable one pro- 
ceeds in a similar manner. One begins by fixing the central enthalpy h c . (This also 
fixes the central pressure and energy density as a result of the equation of state 
and Eq. ( |299| )) . One then integrates from h = h c out to the surface of the star, 
h = 0. This method, unlike the usual approach, has the advantage that the domain 
of integration is known from the outset. 

Using 



one can write Eq. ( |172| ) as 



or 



v (h) -v c = h c -h (302) 

where v c is an arbitrary integration constant. 

Lindblom |46j has re-expressed Eqs. (|170|) and (|171|) as, 

dr r(r - 2M) 



dh (M + 4vrr 3 p) 

and 



(303) 



dM 2 dr . . 

llh =i * re Th> (304) 

respectively. Similarly, we write Eq. ( |174j ) as the pair of first order equations 

^ = e (— f±. (305) 

dh dh K J 



and 



dh 



r . 



dr 

Th (306) 
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where we have defined 



/ 



-{v—v c +\) 



duj 
dr 



(307) 



As in the usual approach, the equations fl303j )-( j306l) are singular at the center 
of the star. Therefore, we start the numerical integration near h = h c using the 
following truncated power series solutions (see Lindblom [|6] for the first two of 
these expressions), 



r(h) 



3{h c - h) 



27r(e c + p c ) 



e c - 3p c + -e% 
5 



(K-h) 1 

(e c + 3p c ) J ' 



(30£ 



M(h) = ^e c r 3 (h)i [ l+'^-(h c -h)} 



3e, 



(309) 
(310) 



16n { 5 

f{h) = —{e c + p c )uj c r{h)ll + - 



6(2e c - 3p c 



+ 



(/i c - h) 



5(e c + 3p c ) {tc + Vo 
The central energy density e c , the central pressure p c and the constant 



de 
dh 



h=h c 



(311) 



(312) 



are determined from the equation of state, while uj c is a normalization constant which 
must be fixed arbitrarily to begin the integration. (We simply set ou c — 1.) 

Beginning with these initial conditions, we integrate Eqs. (|303|) -( [306l ) to the 
surface of the star, h = 0, using a standard Runge-Kutta algorithm (Press et al. 



[32|). The mass and radius of the star are given at this point by M = M(0) and 
R = r(0). The constant v c is then given in terms of these quantities by matching to 
the exterior solution at the surface of the star, 



-h c + u(0) = -h c + -\n[l 



2M 



(313) 



Finally, the angular velocity Q and the angular momentum J of the star are deter- 
mined by matching the interior solution to the exterior solution ( [L76| ). (See, e.g., 



Hartle and Thorne p0|.) At the surface of the star, h = 0, we have, 



(314) 
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and 



which gives 



and 



— = — = e-»°f(0) 
dr R* J{ j ' 



n 



1 du 



r=R 
1 



J=^R A e-^f(0). 

D 



(315) 
(316) 
(317) 



Once these quantities are determined, we may then set Q — 1 by renormalizing the 
solutions u(h) and f(h) (Hartle and Thorne ||30|| ) . 

For the special case of polytropic stellar models, one can explicitly relate the 



energy density and pressure to the enthalpy using Eq. ( |299| ). The equation of state 
for a polytrope is given by 



n+l 



V 
e 



Kp- 

p + np 



(318) 
(319) 



where n is the polytropic index, p is the rest-mass-energy density and K is a constant. 
Since K has dimensions (length) 2//n , it is convenient to work only with dimensionless 
quantities by making the redefinitions, 



V 

P 
e 



K~ n p 
K- n p 
K- n e. 



The equation of state then becomes 



P 
e 



and Eq. ( |299j ) gives 



h{p) 



In 



p n 

p + np 



dp' 
e(p') + p' 

1 + (n + l)p^ 



(320) 
(321) 



(322) 
(323) 
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or 



p{h) 
p(h) 
e(h) 



n+l 




1 + n 



'e h -l 
n + l 



(324) 
(325) 
(326) 



These expressions may then be used in the integration of Eqs. ( |303|) -( |306l ) with the 
initial values 



h 1 \ n 

e hc - 1 \ 



Pc 




1 + n 



1 



n 



n+l 



ne 



n+l 



(327) 
(328) 
(329) 



Finally, we comment on one last aspect of the numerical solution of the equilib- 
rium equations. In order to implement our numerical approach to the perturbation 
equations, we need the power series expansions of the equilibrium quantities that ap- 
pear in these equations. The equilibrium quantities we will require (see Sect. 4.1.2 
below) and the forms of their series expansions are as follows. About the center of 
the star, r = 0, we write 



r(e' + p') 
(e+p) 



e 2A 

TV 1 

rX' 

n 

E 
n 



2; 



r \ 

" 1 r) 



2 1 



E 



2/ 



2; 



and about the surface of the star, r = R, we write 



r(e' + p') 
(e + p) 



(330) 
(331) 
(332) 
(333) 
(334) 
(335) 

(336) 



where ' = d/dr and where we have defined the function 



fi(r) = re 2u (iue- 2u )'. (342) 

All of these quantities are expressible in terms of the functions r(h), M(h), u(h), 
00(h), f(h), p(h) and e(h) defined by our Runge-Kutta solution to the equilibrium 
equations. The coefficients of these series, 7T», vr^ etc., are extracted from our nu- 
merical solution using a polynomial fitting algorithm (Press et al. []62|). We note 



that although the polynomial fit accurately reproduces the various equilibrium func- 
tions defined above, it does not find the coefficients of the higher order terms in 
the fitting polynomials with great accuracy. Since it is the coefficients, themselves, 
that are used in the numerical solution of the perturbation equations, this limits the 
accuracy of our solutions. 

4.1.2 Numerical solution of the perturbation equations 

Having described our numerical method for the solution of the equilibrium equations, 



we now turn to the perturbation equations ( |154j ), ( |156| ), ( p.57|) and ( [212j )- (|214|) . These 
must be solved numerically subject to the boundary, matching and normalization 
conditions (|2T8| ), ( p3|) , ( g2g ) and (|225| ) at the surface of the star. 

Since Eq. ( |154[ ) merely expresses ifi,z( r ) in terms of Wi(r), we may ignore this 
metric variable and drop Eq. (|154|) from the system of equations to be solved, for 



all I. In addition, we need not work with Eq. (|214|) since (|212|) - (|214|) are related 



by Eq. ( |204j ) . Thus, a complete set of perturbation equations is provided by Eqs. 
( 156 ), ( |157| ), ( |212j ) and ( |213j ) which we will now re-express in a form more suitable 



for our numerical approach. 
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We write Eq. ( |156|) as 

= r 2 h'i - r(v' + X)rh! x - [2 + (2 - I 2 - l)e 2X + r(v' + A')l h 



(343) 



-4r(v' + \')U h 



Eq. ([157]) as 



= rW{ + 



1 + r(v' + A') + 



r(e' + p') 
(e + p) 



Eq. (HI) as 



= [1(1 + l)K,n(hi + t/i) - 2mc^] 



(344) 



(345) 



+ (/ + l)Qt [(2uj + n) Wi-t - 2(1 - l)uVi-i] 



-lQi +1 [(2u> + li) W l+1 + 2(1 + 2)cuV l+1 } 



and Eq. ( [213j ) as 



= (/ - 2)1(1 + l)Q,_iQ, \-2iorUi 2 + 2(Z - l)a;^_ 2 + (I - 3)//*7 z _ 



(346) 



+(/ + i)q* 



(/ - l)lKQ,rV{_ x - 2mlurV l '_ 1 - 2(1 - VjlnVtrv'V^ 



+2m(l - l)luVi-! + m{l - 3)Z/xV^_i 



(/ - l)Z/cfie 2A Wi_i + 4«fir(i/ + A')W, 



+Z(Z + 1)| mnVLr (h\ + [//) - 2m/dW (hi + [/;) 

+ ((I + 1)Q 2 - IQ 2 +1 ) (2QrU[ + IpUx) 



+ 



+ Z(Z + l)(Q2 +1 + Q2-l)] (2o) + /i)^} 
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-IQ 



i+i 



(I + !)(/ + 2)/dW ; ' 1 + 2m(l + l)QrV{ +1 - 2(1 + !)(/ + 2)/dWV m 



+2m(l + !)(/ + 2)QV l+1 + m(l + !)(/ + 4)//^+i 



-(/ + !)(/ + 2)«fie 2A Wi + i + 4«flr(i/ + A')Wi+i 



+Z(Z + !)(/ + S)Q l+l Q l+2 2ojrU' l+2 + 2(Z + 2)uU l+2 + (I + A)fiU, 



1+2 



where we have used the definition (|342|) of the function n(r). 

Eqs. ( |343| )-( |3"46|) comprise a system of ordinary differential equations for the 
variables hi', Uy, Wy and Vy (for all I'). Together with the boundary and matching 
conditions at the surface of the star, these equations form a non-linear eigenvalue 
problem for the parameter k, where kQ is the mode frequency in the rotating frame. 

To solve for the eigenvalues we proceed exactly as in the newtonian case (see Sect. 
2.4). We first ensure that the boundary and matching conditions are satisfied by 
expanding hi>(r), Ui>(r), Wi>(r) and Vj'(r) (for all I') in regular power series about the 
surface and center of the star. (We present these expansions explicitly in Appendix 
D, Eqs. (|403|) - (|4ToD ). Substituting these series into the differential equations (|343|) - 



( 346 ) results in a set of algebraic equations for the expansion coefficients. These 
algebraic equations may be solved for arbitrary values of k using standard matrix 
inversion methods. For arbitrary values of k, however, the series solutions about 
the center of the star will not agree with those about the surface of the star. The 
requirement that the series agree at some matching point, < r$ < R, then becomes 
the condition that restricts the possible values of the eigenvalue, Kq- 

We begin by replacing all of the equilibrium quantities in Eqs. ( |343D -( |3"46| ) with 
their series expansions, (|330p-(|3"¥TD and by replacing our perturbation variables hi>(r), 



Ui'(r), Wv(r) and Vv(r) with their series expansions, (^03|) -( f4d0| ). The result is a set 
of algebraic equations for the expansion coefficients hi, hk etc., which we present 
explicitly in Appendix D, Eqs. ( |417| )-( ]42"4| ). We then write down the matching 
conditions at the point r equating the series expansions about r = to the series 
expansions about r = R. (These matching conditions are also presented explicitly 
in Appendix D, Eqs. (|412|) - (|416|) .) The result is a linear algebraic system which we 
may represent schematically as 

Ax = 0. (347) 
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In this equation, A is a matrix which depends non-linearly on the parameter k, and 
x is a vector whose components are the unknown coefficients in the series expansions 
for hi'(r), Ui'(r), Wi>(r) and Vp{r). 

To satisfy equation (|347|) we must find those values of k for which the matrix 
A is singular, i.e., we must find the zeroes of the determinant of A. We truncate 
the spherical harmonic expansions of 5u a and h a p at some maximum index Z max and 
we truncate the radial series expansions about r = and r = R at some maximum 
powers z max and /c max , respectively. 

The resulting finite matrix is band diagonal. To find the zeroes of its determinant 
we use the same routines from the LAPACK linear algebra libraries (Anderson et al. 
[|I|) and root finding techniques that we used in the newtonian calculation. 

The eigenfunctions associated with these eigenvalues are determined by the per- 
turbation equations only up to normalization. Given a particular eigenvalue, we 
find its eigenfunction by replacing one of the equations in the system ( |347| ) with the 
normalization condition ( |225| ). Since this eliminates one of the rows of the singular 
matrix A in favor of the normalization equation, the result is an algebraic system of 
the form 

Ax = b, (348) 

where A is now a non-singular matrix and b is a known column vector. We solve this 
system for the vector x using routines from LAPACK and reconstruct the various 
series expansions from this solution vector of coefficients. 



4.2 Eigenvalues and Eigenfunctions 

We have computed the eigenvalues and eigenfunctions for uniform density stars in 
full general relativity. If Mo is the gravitational mass of the spacetime and R the 
coordinate radius of the star, then the dimensionless constant (2M /R) is a useful 
measure of the strength of relativistic effects. We have studied a set of axial- and 
polar-hybrids over a range of values of this constant from the nearly newtonian to 
the relativistic regimes, 10~ 6 ^ (2M /R) ^ 0.2. 

Although the code is fully relativistic and was written to handle any polytropic 
equation of state, it does not yet accurately find the modes for compressible stars 
(polytropic index n ^ 0) or for highly relativistic uniform density stars (2Mq/R ^ 
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0.2). The problem appears to be related to the fitting of the equilibrium variables 
to their power series expansions, and the difficulty in accurately computing the 
coefficients of the higher order terms in the fitting polynomials (see Sect. 4.1.1). 
Work is in progress to address this problem and extend this study to the broader 
class of stellar models. However, we do not expect qualitative differences between 
our current results and the results of such a study. 

The structure of the hybrid mode spectrum in relativistic stars appears to be 
identical to that in newtonian stars. We find that for each of the hybrid modes 
considered, there exists a family of relativistic modes parametrized by (2M /R), 
whose limit as (2Mq/R) — > is the corresponding newtonian mode. 

By means of the post-newtonian solution of Sect. 3.4.2, we have already verified 
our claim that the newtonian r-modes with / = m > 2 do not exist in relativistic 
stars as purely axial modes. We may now make further use of this analytic solution 
to test the accuracy of our relativistic code in the small (2Mq/R) regime. We do so 
by examining the difference between the frequencies of corresponding modes in the 
newtonian and relativistic stars. If <7o is the frequency of a mode in the newtonian 
star with eigenvalue kq and a is the frequency of the corresponding relativistic mode 
with eigenvalue k, then this difference is given by, 

a - a = (ft - «o)fl. (349) 

In Fig. [13| we display function of (2M /R) for the modes whose newtonian 

limit is a pure r-mode with 2 < I = m < 5. The curves show the post-newtonian 
solution given by Eq. ( |290| ), while the symbols show the fully relativistic numerical 



solution. At (2M /R) = 10% the fractional change in the frequency from the new- 
tonian value is approximately 4% for these modes. Observe that not only do the 
analytic and numerical solutions agree in the small (2Mq/R) regime, but that the 
code scales correctly with (2M /R). 

This agreement with our analytic solution gives us confidence that our code is 
able to find the relativistic modes. Thus, we may now explore those modes for 



which we have not worked out a post-newtonian solution. In Fig. [IJ we display 
K — no as a function of (2Mo/R) for a set of modes whose newtonian limits are 
axial- and polar- led hybrids with m = 2. The modes whose frequency corrections 
are shown correspond to the first six entries in the m = 2 column of Table Q At 
(2Mo/R) = 10% the fractional change in the frequency from the newtonian value 
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is, again, approximately 4%. Fig. [14] also reveals the feature that the frequencies of 
modes with k,q < increase in the relativistic star, while the frequencies of modes 
with k > decrease. So the magnitude of the frequency always decreases as the 
star becomes more relativistic. It is natural that general relativity will have such 
an effect for two reasons. Gravitational redshift will tend to decrease the fluid os- 
cillation frequencies measured by an inertial observer at infinity (i.e., frequencies 
measured with respect to the Killing time parameter, t). Also, "the magnitude of 
the centrifugal force is determined not by the angular velocity Q of the fluid relative 
to a distant observer but by its angular velocity relative to the local inertial frame, 
u)(r)." (Hartle and Thorne p0||.) Thus, the centrifugal and Coriolis forces diminish 
- and the modes oscillate less rapidly - as the dragging of inertial frames becomes 
more pronounced. It is not clear what effect these small frequency shifts will have 
(if any) on the stability of the modes. 

We now turn to a discussion of the eigenfunctions. From Sect. 2.3.1, we know 
that the radial behaviour of the I = m newtonian r-modes is given by U m (r) = 
(r/R) m+l , and in Sect. 3.4.2 we computed their post-newtonian corrections. The 
mode expected to dominate the gravitational- wave instability in a hot, young neutron 
star is the newtonian r-mode with I = m = 2. In Fig. |15|, we display the functions 
Ui(r), Wi(r), and Vi(r) for I < 4 associated with this mode in a uniform density 
star with (2M /R) = 0.2. These components of the eigenfunction were computed 
numerically using the fully relativistic code. The 20% corrections to the structure 
of the equilibrium star induce only 1% corrections to the character of this mode; 
thus, the function U2{r) is barely distinguishable from its newtonian form, (r/R) 3 , 
and the coefficients shown of the other axial and polar terms have been multiplied 
by a factor of 100 to make them visible on the scale of the plot. These functions 
are normalized so that ^(r) = 1 at the surface of the star, r = R. In Fig. [16] we 
show the metric functions hi(r) for I < 6 for the same mode. The vertical scale is 
set by the normalization of U2{r), and this implies that the metric perturbation is 
at most a 4% correction to the equilibrium metric. The fact that Z^O") dominates 
the perturbed metric is the statement that this mode couples strongly to current 
quadrupole radiation. 

Figs. |17| - |19| are a series of plots displaying the m = 2 axial-led hybrid mode 
whose newtonian eigenvalue in the uniform density star is k = 0.466901. (See Figs. 

and Table |3|.) Fig. |I7| shows the functions Ui(r), Wi(r), and Vi(r) for this mode 
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as calculated by the relativistic code in the newtonian regime, (2Mq/R) = 10~ 6 . As 
expected, these functions agree with the newtonian forms displayed in Figs. 
up to the difference in the normalization conditions used in the newtonian (|89|) and 
relativistic (|225|) calculations and corrections of order 10~ 6 . In Fig. [18| we display 
the functions Ui(r), Wj(r), and Vi(r) for I < 6 for the same mode, but now with 
(2Mq/R) = 0.1, and in Fig. [19] we display the corresponding metric functions hi(r) 
for I < 6. Observe that U2, W3, V3 and U4 are barely distinguishable from their 
newtonian forms and that the mode has acquired relativistic corrections of order 



1%. Fig. ^ reveals the interesting feature that the metric function /12O") nearly 
vanishes in the exterior spacetime. Recall the discussion at the end of Sect. 2.6 of 
the vanishing of the / = 2 current multipole for this particular mode in the uniform 
density newtonian star; we saw that the only nonzero current multipole is that with 
I = 4. Here, the relativistic calculation reveals explicitly that the perturbed metric 
in the exterior spacetime is dominated by h^r), and that h,2(r) is smaller by two 
orders of magnitude in the exterior spacetime. 

Finally, we present a series of plots displaying a polar-led hybrid mode and its 
relativistic corrections. Figs. ^0|-|22| show the m = 1 polar-led hybrid mode whose 
newtonian eigenvalue in the uniform density star is k,q = 1.509941. (See Figs. f§-[| 
and Table |^ Fig. ^0] shows the functions Ui(r), Wi(r), and VJ(r) for this mode as 
calculated by the relativistic code in the newtonian regime, (2M /R) = 10~ 6 . As 
expected, these functions agree with the newtonian forms displayed in Figs. [|-|6|, 
up to the difference in the normalization conditions used in the newtonian (|89|) and 
relativistic ( 225 ) calculations and corrections of order 10~ 6 . In Fig. |21| we display 
the functions Ui(r), Wi(r), and Vi(r) for I < 4 for the same mode, but now with 
(2M /R) = 0.05, and in Fig. we display the corresponding metric functions 
hi(r) for I < 4. Observe that W\, V% and JJi are barely distinguishable from their 
newtonian forms and that the mode has acquired relativistic corrections of order 



1%. Fig. 22 shows that h 2 (r) dominates the metric perturbation. However, since 
the mode is stable, this strong coupling to current quadrupole radiation will serve 
only to damp the mode rapidly. 
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Figure 13: The difference between the relativistic and newtonian eigenvalues (k — 
Kq) as the uniform density star becomes increasingly relativistic. These frequency 
corrections are shown for the modes whose newtonian limit is a pure r-mode with 
2 < I = m < 5. The curves were calculated using the analytic post-newtonian 
expression ( |290| ), while the symbols were calculated numerically. 
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Figure 14: The difference between the relativistic and newtonian eigenvalues (re — 
Kq) as the uniform density star becomes increasingly relativistic. These frequency 
corrections are shown for a number of both axial- and polar-led hybrid modes with 
m = 2 (see Table |5|). All of the data points were calculated numerically. 



114 



0.5 



-0.5 



U 2 (r) 

W 3 (r) x 10 2 
V 3 (r) x 10 2 
U 4 (r) x 10 2 



\ r 
\ \ 
\ 

V 



-1 



0.2 



0.4 0.6 

(r/ R) 



0.8 



Figure 15: Coefficients Ui(r), Wi(r), and Vi(r) with / < 4 of the spherical harmonic 
expansion ( |177[ ) for the m = 2 axial-led hybrid mode whose newtonian limit is the 
pure axial r-mode with I = 2 and co moving frequency k,qQ = 2Q/3. The mode 
is shown in the uniform density star with (2Mq/R) = 0.2, for which its comoving 
frequency has shifted to kQ = 0.625fi. The vertical scale is set by the normalization 
of U2(r) to unity at the surface of the star. Observe that there are both axial and 
polar relativistic corrections of order 1%. The coefficients of expansion (|177 ) with 
I > 4 are of order 0.01% and smaller and are not shown. 
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Figure 16: Coefficients hi(r)(= ho j i{r)) with I < 6 of the spherical harmonic expan- 
sion Q179Q of the perturbed metric for the same mode shown in Fig. [15]. The vertical 
scale is the same as that of Fig. [15] and is set by the normalization of t^tj)- Observe 
that, as expected, h,2(r) dominates the perturbed metric, which implies that this 
mode couples strongly to current quadrupole radiation. 
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Figure 17: Coefficients Ui(r), Wi(r), and Vi(r) with / < 4 of the spherical harmonic 
expansion ( |177| ) for the m = 2 axial-led hybrid mode whose newtonian limit has 
comoving frequency k,qQ = 0.466901fi. The mode is shown in the uniform density 
star with (2Mg/R) = 10 -6 to check for agreement with the newtonian calculation. 
As expected, these functions agree with the newtonian forms displayed in Figs. 0-|9|, 
up to the difference in the normalization conditions used in the newtonian (|89|) and 
relativistic fl225| ) calculations and corrections of order 10~ 6 . 
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Figure 18: Coefficients Ui(r), Wi(r), and Vi(r) with I < 6 of the spherical harmonic 
expansion ( |177| ) for the same mode shown in Fig. [L7| but now with (2M Q /R) =0.1. 
The comoving frequency has shifted to kQ = 0.443f2, and the mode has acquired 
relativistic corrections of order 1%. Accordingly, the functions U2, W3, V3 and U4 are 
indistinguishable from those shown in Fig. [17|. The coefficients of expansion (|177|) 
with I > 6 are of order 0.01% and smaller and are not shown. 
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Figure 19: Coefficients hi(r)(= h i(r)) with I < 6 of the spherical harmonic expan- 
sion ( |1 79| ) for the same mode shown in Fig. 
(2Mo/R) = 0.1. The vertical scale is the same as that of Fig. [H 
normalization of L^r). Observe that h,2(r) nearly vanishes in the exterior spacetime 
and that the metric perturbation is dominated by h^r). This is the expected result 
that the mode couples strongly to the I = 4 current multipole only (see Sect 2.6). 



T8| in the uniform density star with 
and is set by the 
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Figure 20: Coefficients Ui(r), Wi(r), and Vi(r) with / < 2 of the spherical harmonic 
expansion ( |177| ) for the m = 1 polar-led hybrid mode whose newtonian limit has 
comoving frequency k,qQ = 1.5099411). The mode is shown in the uniform density 
star with (2Mg/R) = 10 -6 to check for agreement with the newtonian calculation. 
As expected, these functions agree with the newtonian forms displayed in Figs. |]-|6|, 
up to the difference in the normalization conditions used in the newtonian (|89|) and 
relativistic fl225| ) calculations and corrections of order 10~ 6 . 
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Figure 21: Coefficients Ui(r), Wi(r), and Vi(r) with / < 4 of the spherical harmonic 
expansion ( 177|) for the same mode shown in Fig. [20] but now with (2M /R) = 0.05. 
The comoving frequency has shifted to kQ = 1.4802fi, and the mode has acquired 
relativistic corrections of order 1%. Accordingly, the functions W\, V\ and U2 are 
indistinguishable from those shown in Fig. ^0|. The coefficients of expansion (|177|) 
with I > 6 are of order 0.01% and smaller and are not shown. 
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Figure 22: Coefficients hi(r)(= ho j i{r)) with I < 4 of the spherical harmonic expan- 
sion ( |179| ) for the same mode shown in Fig. [21] in the uniform density star with 
(2Mq/R) = 0.05. The vertical scale is the same as that of Fig. ^1] and is set by the 
normalization of t^fV). The metric perturbation is dominated by h2(r), which im- 
plies that this stable mode will be rapidly damped by current quadrupole radiation. 
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Appendix A 



Proof of Theorem \t\ 



A.l Axial-Led Hybrids with m > 

Let I be the smallest value of I' for which Uy 7^ in the spherical harmonic expansion 
(|50|) of the perturbed velocity field 6v a . The axial parity of 5v a , (— 1) , and the 
vanishing of Y™ for / < m implies I > m. That the mode is axial-led means Wy = 
and Vy = for I' < I. We show by contradiction that I = m. 

Suppose I > m + 1. From equation fl65f ), J g' r YJ* m <if2 = 0, we have 
1 



-«/(/ + 1) 



m 



U l = lQ l+1 [W l+1 + (l + 2)V l+1 ] 



and from equation (|66| ) with Z replaced by I — 1, / q Y{^dVt = 0, we have 



Q m (J + 2)1? +1 + 



m+-«(Z + l) 



U[ + m(l + l) 



U, 



(350) 



(351) 



These two equations, together imply that 

V[ + -U, 

r 



0. 



or 



U, = Kt~\ 



which is singular at r = 0. 



A. 2 Axial-Led Hybrids with m = 



Let m = and let I be the smallest value of I' for which t/j/ 7^ in the spherical 
harmonic expansion (EDI) of the perturbed velocity field 5v a . Since V a l^° = 0, the 
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mode vanishes unless I > 1. That the mode is axial- led means Wy 
for /' < I. We show by contradiction that 1 = 1. 
Suppose I > 2. Then f qVY&dtl = becomes, 



and V v = 



or 



which is singular at r = 0. 



r 



0. 



U, = Kr~\ 



(352) 



A. 3 Polar-Led Hybrids with m > 

Let I be the smallest value of I' for which W\> 7^ or Vy 7^ in the spherical harmonic 
expansion (|50f) of the perturbed velocity field <5f a . The polar parity of 5v a , (—1)', 
and the vanishing of Y™ for I < m implies I > m. That the mode is polar-led means 
Uy = for /' < /. We show by contradiction that / = m. 
Suppose / > m + 1. Then / q r Y l *^dQ = becomes 



Wi + {l + l)Vi = 0, 



(353) 



and / (fYffidto = becomes, 



-k(1 + l)+m 



V, 



V/ + m(l + l)^ 



2 r 



(354) 



+ (l + 2)Q 



1+1 



U' l+l + (/ + !) 



U l+1 



These two equations, together imply that 



~-k(1 + l) + m 



Vi 



+ (l + 2)Q. 



1+1 



ul +1 + (/ + 1 



u l+1 



or 



-k(1 + 1) + m 



V t + {1 + 2)Q l+1 U l+1 = Kr 



-(i+i) 



which is singular at r = 0. 
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A. 4 Polar-Led Hybrids with m = 



Let m = and let I be the smallest value of I' for which Wy ^ or Vy ^ in the 
spherical harmonic expansion (|50D of the perturbed velocity field 5v a . When I = 
the mode is automatically polar-led; thus we need only consider the case I > 1. That 
the mode is polar- led means TJy = for V < I. We show by contradiction that 1 = 1. 
Suppose I > 2. Then / q r Y l *\dVL = becomes 



(355) 



and / q^Y^dfl = becomes, 



-k{1 + 1) 



V, 



, w l 



+ {l + 2)Q 



i+i 



U[ +l + (I + 1 



Ui+i 



. = 



(356) 



These two equations, together imply that 



--k(! + i) 



v{ + (i + i)— 

r 



+ + 2)Q l+1 



uu + ii + i) 1 ^ 



or 



--k(1 + m + + 2)Q l+1 U l+1 = Kr-V +l \ 



which is singular at r = 0. 
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Appendix B 

Algebraic Equations: Newtonian 



In this appendix, we make use of the following definitions: 



a t = ~Km+(l + l)Qf-lQt +1 
bl = m 2 - i(i + 1) (l - ^ - QiVa 



-«/(/ + 1) 



For reference, we repeat the definitions (^) and (p0[ ) 

_ (cr + mfi) 

(I + m)(l — m) 



Qi 



(21 - 1)(2Z + 1) 



(357) 
(358) 
(359) 



(360) 
(361) 



B.l Axial-Led Hybrids 



For I = m, m + 2, m + 4, . . . the regular series expansions 1 about the center of the 
star, r = 0, are 



V m+j+ i(r) 



R 



i odd 



r \ 



-i odd 



(362) 
(363) 



1 We present the form of the series expansions for Ui(r) for reference; however, we do not need 
these series since we eliminate the TJ\ (r) using equation (pq) . 
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/ r \ m+j oo , v i 

U m+J (r) = [-) E«i,i(^J ( 364 ) 

i odd 

where j = 0, 2, 4, . . .. 

The regular series expansions about r = R, which satisfy the boundary condition 
Ap = are 

W m+j+ i(r) = 53iS i+1>fc ^l-_j (365) 

Kn+ J+ i(r) = E fy+i,* 1 - d ( 366 ) 
fc=o V ^ y 

OO / r v k 

U m+j (r) = ^2u jtk {^l - - j (367) 



k=0 v 

where j = 0, 2, 4, . . .. 

These series expansions must agree in the interior of the star. We impose the match- 
ing condition that the series (|362|) -( |363| ) truncated at i max be equal at the point r = r 
to the corresponding series ( |365| )-( P6"6| ) truncated at /c max . That is, 

/ r \ Tn-\-j *raax / r \ i fcffiax / r \ k 

0=1?) E »i+M ( ?) " E ( 1 - £ ) (368) 



-ft/ , \ -ft / 71 \ -ft/ 



z odd 



r i) e^-eM^) <*») 



i=l x ' fc=0 

z odd 



When we substitute (M -(363|) and (H) into (H), the coefficient of (r/R) m+ ^ +t in 
the resulting equation is 

i-2 

= (to + j + i + + E w j+i,i-s-t (370) 

3=1 
s odd 

-(m + j + l)(m + j + 2)v j+1>i 



Similarly, when we substitute ( |365| )-( |3"6T)| ) and (WJt) into ([84]), the coefficient of [1 — 
(r/R)] k in the resulting equation is 

k 

= (A; + 1) [wj+i,k ~ Wj+i,k+i] + E ~ ^- 2 ) *j+i,fe-«+i ( 371 ) 

s=0 



-(to + j + 2)(m + j + l)%+i,fe 



128 



where we have defined 7r_2 = = Wj+ifl. 

When we use (|65| ) to eliminate the Ui(r) from ( ]67| ) and then substitute for the Wi±i{r) 
and Vi±i(r) using ( |362| )- (|363|) , the coefficient of (r / ' B) m+ i +l in the resulting equation 
is 

= (i + 1)(to + j - 2)(m + j - l)<2 m+i Q m+j _iQ m+i _ 2 c m+i c m+j+ 2 (372) 



x 



^i-3,i+4 - (m + j - 3)uj_ 3i 



i+4 



- < {i + l)(m + j - 2) 2 Q 2 m+j _ 1 c m+j + §«(m + j - l)c m+j _ 2 c 



+ (m + j + 1) [(m + j + i)a m+i + 6 m+j ] c m+j _ 2 



i+2 



+ 



|/c(to + j — l)(m + j + i) — (i + l)m 



Cm+j — 2 Cm+j 



+ (m + j + l)(m + j - 1) [(rn + j + i)a m+j + b m+j ] c m+i . 



■(z + l)(m + j)(m + j - 2) 2 <& +j _ 1 



xQm+j c m+j+2 v j-l,i+2 



|/c(m + j + 2)c m+j c m+J+2 



+(to + j) [(to + j + i)a m +j + 6 m+ j] c m+i+2 



-(2m + 2j + i + 2)(m + j + 3) 2 Q 2 m+j+2 c m+J 



x Qm+j+l c m+j~2 w j+l,i 
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+ < (m + j + 2)(m + j) [(m + j + i)a 



\K(m + j + 2)(m + j + i) + m(2m + 2j + z + 2) 



Cm+jCm+j+2 



+ (2m + 2j + i + 2)(m + j + 3) 2 (m + j + l)Q 2 m+j+2 c, 



+j+2^m+j 



x Qm+j+lC m +j-2 v j+l, 



+ (2m + 2j + i + 2)(m + j + 3)(m + j + 2)Q m+i+3 

Qm+j +2 Qm+j +1 



^ ('m~j - ~ j 



w i+3ii _ 2 + (m + j + 4)uj +3j . 



i-2 



When we use (|65|) to eliminate the {7i(r) from (|67f) and then substitute for the Wi±i{r) 
and V/±i(r) using (|365|) - (|366|) , the coefficient of [1 — (r / ' R)] k in the resulting equation 
is 



= -(m +j-k-l)(m + j-l)(m + j- 2)Q m+j Q m+j _ 1 Q m+j _ 2 (373) 



X C m _|_jC m _)_j_)_2 



%-3,fe - (m + j - 3){; i _ 3>fc 



(A; + l)(m + j - l)(m + j - 2)Q m+j (5 m+J -_i(5 m+j _2C m +jC m+J - + 2 



x 



w. 



-3,fc+i — (m + i — 3)%_ 3ifc+ i 



(m + j - fc - l)(m + j - 2) 2 Q 2 m+j _ 1 c m+j 
— \n{rn + j — l)c m+ j_2 c m+j 



-(to + j + 1) (b m+j + ka m+j ) 



*Qm+j c m+j+2 w j-l,k 



(m + j- 2) 2 Q 2 n+j 

— 1 Cm+j 

+ (m + j + l)a m+i c m+j _ 2 f 



x(fc + l)Qm+jC m+i+2 Wj-i jk+1 



(m + j — k — l)(m + j — 2) 2 (m + j)<5^ +J _iC m+j 



\nk{m + j — 1) + m(m + j — k — 1) 



Crn+j— 2Cm-\-j 



+ {m + j + l)(m + j - 1) (6 m+i + ka m+j ) 

Wn+j-2 



xQm+j c m+j+2 v j-l,k 



c m+j 



+ 



m — \K{m + j — 1) 



Cm+j— 2Cm+j 



-(m +j + l)(m + j- l)a m+j c m+j - 2 



X{k + l)Q m +i c m+'j+2 "^jr'-l.fc+l 



(m + j) (6 m+j + fca m+J ) c m+j+2 



+ |«;(m + j + 2)c m+ jC m+ j + 2 



— (m + j + k + 2)(m + j + 3) 2 Q 2 m+j+2 c m+j 



x Qm+j+lCm+j-2 Wj + \ t k 
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+ < — (to + j)a m +jC m +j+2 + (m + j + 3) 2 Q 2 m+j+2 c m+j 



+ { (m + j + 2)(m + j) {bm+j + ka m+j ) Cm+j 



+2 



m(m + j + k + 2) + ^nk(m + j + 2) 



Cm+j Cm+j +2 



+ (m + j + k + 2)(m + j + 3) 2 (m + j + l)Q 2 m+j+2 c m+j 



X-Qm+j+l c m+j—2 Vj+l,k 



+ \ - (m + j + 2) (m + j)a m+i c m+j+2 



|k(to + j + 2) + m 



Cm+j C-m+j+2 



-(m + j + 3) 2 (m + j + l)Q 2 „ + , +2 c m+j 



x(/c + l)Q m+ j + iC m+ j^2 



+ (to + j + fc + 2) (to + j + 3) (to + j + 2)Q m+j+3 Q m+j+2 Qm+j+i 



X Cm+j—2Cm+j 



(k + \){m + j + 3)(m + j + 2)Q m+j+3 Q m+j+2 Qm+j+lCm+j-2Cm+j 



w j+3tk+1 + (m + j + 4)v j+3jk+1 



The equations (|368|) through (|373|) make up the algebraic system (|88|) for eigen- 
values of the axial-led hybrid modes. One truncates the angular and radial series 
expansions at indices j max , i max and /c max and constructs the matrix A by keeping 
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the appropriate number of equations for the number of unknown coefficients Wj + \ : i, 
Vj+i t i, Wj + i t k and Vj+i^- In following this procedure, however, one must be aware of 
the following subtlety in the equations. 
For each q = j + i the set of equations 



T0| ) with i = l and j = q — 1, and 
72]) for all i = 1, 3, . . . , q with j = q — i 



can be shown to be linearly dependent for arbitrary k and for any equilibrium stellar 
model. For example, taking the simplest case of q = 1, one can show that equation 
( |370| ) with 2 = 1 and j = becomes 

= (m + 2) Ki-(m + lKJ 

while equation ( [372| ) with i = 1 and j = becomes 

0=| |«(m + 2)c m c m+2 + m [(m + l)a m + b m ] 

-(2m + 3)(m + 3) 2 <5^ +2 c m |Q m+ iC m _ 2 - (m + l)vi,i] • 

This problem can be solved by eliminating one of these equations from the subset 
for each q (for example, equation fl372| ) with i = 1). Thus, to properly construct the 
algebraic system fl88|) we use, for all j = 0, 2, . . . , j max , the equations 

(HD 



( P70D with z = 1,3, . . . ,i max 

( P7ID with fc = 0,l,...,fc max -l 

(|372| ) with i = 3,5, . . . ,z max 

(0731) with fc = 0,l,...,fc max - 1. 



B.2 Polar-Led Hybrids 

For I = m, m + 2, m + 4, . . . the regular series expansions^ about the center of the 
star, r = 0, are 



2 We present the form of the series expansions for Ui(r) for reference; however, we do not need 
these series since we eliminate the Ui(r) using equation (p5[). 
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W m+j {r) 



Vm+j(r) 



R 



m+j °° 



r V 

r) 



E 

! = 

i even 
OO 

E 

i even 



3,i 



R 



r) 5 u ^ [r 



(374) 
(375) 
(376) 



where j = 0, 2, 4, . . .. 

The regular series expansions about r = R, which satisfy the boundary condition 
Ap = are 



W m+J (r) = E^i,4 1_ ^J 
°° / r \ k 



U m+j+1 (r) = ~ ft) 



(377) 
(378) 
(379) 



where j = 0, 2, 4, . . .. 

These series expansions must agree in the interior of the star. We impose the match- 
ing condition that the series (|374|) - (|375|) truncated at z max be equal at the point r = r Q 
to the corresponding series (|377|) -( |37"8"D truncated at /c max . That is, 



8=0 
i even 



E ^ ( r i) - E - 'i) 



l ft max 



m+j 'ma 



2 "-max 



r 



(380) 
(381) 



When we substitute (p?^)-CT) and (§§ into (§§, the coefficient of (r/R) m+j+i in 
the resulting equation is 



i-2 



= (m + j + i + l)wj ti + n s+i 



w 



j,i-a-2 



(m + j)(m + j + l)vj,i (382) 



s=0 
s even 
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Similarly, when we substitute ([377|) -( |378|) and ( jSTj) into fl84]), the coefficient of [1 
(r/R)] k in the resulting equation is 



= (k + 1) [w j>k - Wj,k+i) + ~ ^-2) w i|fe _ s+ i (383) 



s=0 



-(m + j)(m + j + l)v jjk 
where we have defined 7r_ 2 = = Wj j0 . 

When we use (|65|) to eliminate the Ui(r) from (|66|) and then substitute for the Wi±i{r) 
and Vy-i(r) using ( |374| )- ([375|) , the coefficient of (r / R) m+ i +l in the resulting equation 
is 



= -im(m + j - l)Q m+j Q m+j -iC m+j+1 



(384) 



Wj-2,1+2 ~ (m + j - 2)Vj-. 2 ,i+2 



+ < [(i + l)m- -K(m + j - l)(m + j + i)](m + j - l)Q 2 m+j c m+j+1 



+ 



[m 



m+j Qm+j+1 



-urn 
2 



C-m+j — 1 C771+ j + 1 



[m(2m + 2j + i + 2) + -«(m + j + 2)(m + j + i)\ 



x ( m + J + 2)g2 



+ <[(i + l)m - -ft(m + j - l)(m + j + «)] 



x (m + j - l)(m + j + l)Q 2 m+i c 



m+j c m+j+l 



m + (m + j + i)a m+j 



Cm+j — 1 Cm+j+ 1 



m(2m + 2j + i + 2) + -K,(m + j + 2)(m + j + i) 



x (m + j){m + j + 2)Q 2 m+j+1 c m+j - 1 \v jt 



+ Q m+j+2 Q m+ j + i [m(m + j + i) + m(m + j + l)(2m + 2 j + % + 2)] 
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^ Cm+j-l 



Wj+2,i-2 + (m + j + 3)v j+2 ,i-2 



When we use (|65|) to eliminate the C/j(r) from (|66|) and then substitute for the (r ) 
and V;±i(r) using (|377|) - (|378|) , the coefficient of [1 — (r / R)] k in the resulting equation 
is 



= m(m +j-l)(m + j- k)Q m ^Q m ^ x c m ^ nx 



(385) 



x 



Wj-2,k - (m + j - 2)vj-2,k 



+ (k + l)m(m + j - l)Q m +j Qm+j— l^m+j+l 



Wj-2,fc+i - (m + j - 2)u,-_ 2) fc+i 



+ j -[(g^ + rre)(m + j — 1) — fem](m + j - l)Q 2 m+j c m+j+1 



+ 



KUl + k ( 1 Qm+j Qm+j+1 



£-171+ j — 1 C-m+ j + 1 



[(-k/c + m)(m + j + 2) + /cm] 



Cm+j — l ( 



Wj,k 



(k + l)l[m- -K(m + j - l)](m + j - l)Q 2 m+j c m+3+1 



+ (l Qm+j Qm+j+lj c m+j-l c m+j+l 



[m + -K(m + j + 2)] 
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x 



(m + j + 2)Q 2 m+j+1 c m+j ^w jjk+1 



+ j — \{-nk + m)(m + j — 1) — km] 



x(m + j-l)(m + j + l)Q 2 m+i c 



m+j c m+j+l 



- (m 2 + ka m+j ^ 



+ [(- K k + m){m + j + 2) + /cm] 



x (m + j')(m + j + 2)Q^ +j+1 c m+ j_i|u jifc 



+ + 1)< -[m- -K(m + j - 1)] 



x 



(m + j - l)(m + j + l)Q^ +i c m+i+1 



-|- Q, m +jC m +j—iC m +j+i 



— [m + -K,(m + j + 2)] 



x (m + j)(m + J + 2)Qm +j+ iC m +j-i|%fc+i 



+ m(m + j + 2)(m + j + k + l)(5 m+J+ 2<5m+j+ic m +j-i 



x 



Wj+2,k + (m + j + 3)v j+2 ,k 



(k + l)m(m + j + 2)Q m+j+2 Q 

m+j+l Cm+j — 1 
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w j+2 ,k+i + (m + j + 3)v j+ 2 lk +i 



The equations ( |380| ) through ( |385| ) make up the algebraic system (|88|) for eigen- 
values of the polar-led hybrid modes. As in the case of the axial-led hybrids, one 
truncates the angular and radial series expansions at indices j max , i max and fc max and 
constructs the matrix A by keeping the appropriate number of equations for the 



number of unknown coefficients Wjj, Vj^, Wj^ and Vj t k- 

We, again, find that certain subsets of these equations are linearly dependent for 
arbitrary k and for any equilibrium stellar model. For all j, it can be shown that 
both equation ( |382| ) with % = and equation ( |384|) with i = are proportional to 

= K',0 - ( m + j)Vjfl] ■ 



This problem can, again, be solved by eliminating, for example, equation (384) 
with i = for all j. Thus, to properly construct the algebraic system ( |58| ) we use, 
for all j = 0, 2, . . . , j max , the equations 



(HD 



( gH) with i = 0, 2, . . . , i 

([383]) with fc = 0,l, 

( gg§ with i = 2, 4, 

(BSHT) with k — 0,1, 



k 

5 "'max 



t — 1 

5 "'max 
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Appendix C 



Proof of Theorem 



C.l Axial- led hybrids with m > 0. 

Let I be the smallest value of I' for which Uy ^ in the spherical harmonic expansion 



( 177 ) of the displacement vector £ Q , or for which hy = hoj> 7^ in the spherical har- 
monic expansion ( |179|) of the metric perturbation h a p- The axial parity of (£ a , h a p), 
(— and the vanishing of Y™ for I < m implies I > m. That the mode is axial- 
led means Wy = 0, Vy = and H\y = for I' < I. We show by contradiction that 
I = m. 

Suppose l>m + l. From Eq. (|TJ), J Aue^Y^dQ = 0, we have 

= 1(1 + + C/,) - 2mc^ (386) 



-/Q m [^<9 r (r^e" 2 ") + 2(1 + 2)uV l+ i 
and from Eq. (gig ) with Z -»• Z - 1, / Aw^^tZfi = 0, we have 

= -4 (Z + l)«fid r [e^ih + + 2m<9 r (ue^Ut) + ^±^<9 r (r 2 cue-^) £/, 



-Q 



1+1 



±<9 r (rW 2 ^) W l+1 + 2(1 + 2)d r (ue- 2 »V l+l 



(387) 



Together these give, 



= 2d r (ue-^Ut) + — d r (r 2 uje~ 2v ) U t 



1(/+2) TT 



' f ,-.„-2i>\ 2 



(388) 
(389) 
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or, 



U t = K (Qe-*')- 1 *™ r~ l 
(for some constant K) which is singular as r —>■ 0. 



(390) 



C.2 Axial- led hybrids with m = 0. 



Let m = and let I be the smallest value of I' for which Uy 7^ in the spherical 
harmonic expansion ( |177| ) of the displacement vector £ a , or for which hy = /io,z' 7^ 
in the spherical harmonic expansion (179) of the metric perturbation h a p. Since 



VaFg = 0, the mode vanishes unless I > 1. That the mode is axial- led means 
W\> = 0, Vi> = and ii^j/ = for I' < I. We show by contradiction that 1 = 1. 
Suppose I > 2. Then Eq. (pT5[) with 1^1-2, J /\uj rd Y*%dVt = 0, becomes 



= 2<9 r (be" 2 ^) + 4<9r ( 



r 2 ooe- 2 n Ur 



2- -2u 

r ue 



l ( - -2v\ 2 

r [uje 



Ui 



(391) 
(392) 



or, 



(393) 



(for some constant K) which is singular as r —>■ 0. 



C.3 Polar-led hybrids with m > 0. 



Let / be the smallest value of V for which Wy 7^ or Vy 7^ in the spherical 



harmonic expansion ( |1 77| ) of the displacement vector £ a , or for which Hiy 7^ in 
the spherical harmonic expansion (|1 79| ) of the metric perturbation h a p. The polar 
parity of h a p), (—1)', and the vanishing of Y™ for I < m implies / > m. That 
the mode is polar- led means Uy = and hy = for I' < I. We show by contradiction 
that I = m when m > and that I = 1 when m = 0. 

Suppose I > m + 1. From Eq. ( gTj ) with Z -> / - 1, / Au 9[p Y*™dQ = 0, we have 



= (/ - l)Qi 



e "\l ( 

O r 



r 2 ue 



-2v 



Wi + 2(1 + l)QVi 



(394) 
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Substituting for Vi using Eq. ( |157| ), we find 



~d r (r 2 uje~ 2v ) Wi + 2ue 



2u< 



e + p) 



{t + p)e i - v+X) rWi 



2 ( 



2- — 2u\ 2 

r we ' 



-i(i-2) e"> +A ) 



r 2 (e + p) 



with solution, 



Wi — K (ue- 2u ) 



.1 p-(^+A) 



-(HI) 



(395) 
(396) 



(397) 



(for some constant iT) which is singular as r — > 0. 

When m = this argument fails to establish that I cannot be equal to 1, because 
Eq. ( [394 ) is trivially satisfied for I — 1 as a result of the overall I — 1 factor. Instead, 
the argument proves that I cannot be greater than 1 in this case and therefore that 
I = 1. 
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Appendix D 

Algebraic Equations: Relativistic 



In this appendix, we make use of the following definitions. 

ai = (/ + 1)Q 2 - IQ 2 1+1 (398) 
k ee m 2 - 1(1 + 1) (l - Q 2 - Q 2 +1 ) (399) 



and we repeat the definitions 



K EE ^±^1 (400) 



(/ + m)(Z — m) 



_(2Z - 1)(2/ + 1) 
We will also make use of the definition 



(401) 



f o fc < o 

e(fc) ee J . 402 

U A;>0 

The regular power series expansions of the perturbation variables about the center 
of the star, r = 0, are 



Hr) = EMo > (403) 



=o 



r \ l+l+2i 



Ui(r) = E«m(^J > ( 404 ) 



i=0 

^ / I 2 & 



=0 



W,(r) = E W V d . ( 405 ) 



i=0 



r \ Z+2i 



V^(r) = E w w(jjJ , (406) 
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while about the surface of the star, r = R, they are 

hi(r) 
Ufa) 

VAr) 



00 ~ / r \ k 



r) 



00 , 

J2ui,k (1 



(407) 
(408) 
(409) 
(410) 



The boundary condition ( [218] ) is automatically satisfied by the form of the surface 
expansion of Wi(r), Eq. ( |409|) . On the other hand, The matching condition (|224|) 
is not automatically satisfied by the form of the surface expansion of hi(r), Eq. 
( 407 ). Instead, condition ( |224| ) places the following non-trivial restriction on the 
series expansion fl407|) , 



= hi. 



.s=0 



00 
s=0 



(411) 



where the constants h^ s are given by the recursion relation fl222j ) up to normalization. 
The normalization factor h^Q is then fixed by the matching condition (|223|) once the 
interior solution is known. 

The series expansions about r = must agree with those about r = R everywhere 
in the interior of the star. To ensure this agreement we impose the matching condition 
that the series (^03|) -( [40"6| ) truncated at i max be equal at the point r = r to the 
corresponding series ( f407] )-( [4T0|) truncated at /c max . That is, 

l+l+2i fema 



o = £M^ 



£ hk (1 - ^ 

k=0 v K 



^max 



*max / y \ &max / j* \h 



^max / ^ \ 



l+2i k n 



(412) 

(413) 
(414) 
(415) 
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Furthermore, since the function hi(r) obeys a second order differential equation, we 
must also impose a matching condition on its derivative h[(r); namely, 



l + 2i fcmax / r \ fc-1 

) + E 

i=0 k=l 

We now consider the perturbation equations (|343|) -( p46|) . We substitute for the 



imax / rp \ 1 + 2% ^max _ / ~, \ k — 1 

Q= y E{l + l + 2i)h hi (j^ + E fc v(i-£) • ( 41 6) 



equilibrium variables in these equations their power series expansions ( |33Q) -( j3~4~l"|) . 
We substitute for the perturbation variables in these equations their power series 
expansions ( |403| )-( f4T0|) . Then, by applying straightforward rules for the multiplica- 
tion of power series, we extract the series expansions of the perturbation equations, 
themselves. The requirement that the coefficients of these expansions vanish inde- 
pendently then gives us the algebraic equations for the unknown constants h^i, u^i, 
wi t i, Vi } i, hi^, ui t k, wi t k and vi t k, (for all allowed values of /, i and k) discussed in Sect. 
4.1.2. 

When we substitute the series (p0|)-(g35|) and fl403D -(g06|) into Eq. Q343D the 
resulting series expansion about r = is 



i=i 



° = Ei [(l + 2i)(l + 2i + l)-l(l + l))h hi (417) 
-E [fo+i + x i+i)( l + 2i - 2 J + l) + (l 2 + l- 2)E j+1 ] h lti - s -x 



3=0 



i- 1 } f r \ l+1+2i 



When we substitute the series (p0|)-(p5|) and (p3|)-(M into Eq. (pg|) the 
resulting series expansion about r = is 

00 ( »-i 1 / r \ l+2i 

= E|( z+2i+1 ) w M+X>.7+i+Vi+^^ • ( 41 §) 

When we substitute the series (p0|)-(p5|) and ( g03[) -(g0g) into Eq. fl345l) the 



resulting series expansion about r = is 



i-l 
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= ^< l(l + l)K,hi t i + [l(l + l)K-2rmj ]ui ti -2mJ2ujj +1 ui >i -j- 1 (419) 

i=0 I j=0 



+{l + l)Qi 



i-l 



j=0 



+At i+ i u;i_i,o - 2(/ - 1)^3^ ^_i 5 i_ i+ i 

j=0 



-IQ 



i+i 



i-l 



j=0 



+2(1 + 2)^u; i ^ +M _ j 

j=0 



r \ l+l+2i 

r) 



When we substitute the series (p0|)-(g35|) and fl403D -(g06|) into Eq. Q346|) the 
resulting series expansion about r = is 

oo ^ i— 1 

= E { " 4 (* + 1 )^o^-2, J +i + £[-4(i - J>i+i + ( l - «i-2,i-j(420) 



8=0 



j=0 



+(/ - 3)/i m u,_ 2 ,o)(Z - 2)/(/ + l)Qi_iQ, 
+ { [(/ - 1)1(1 + 2i + 1)k - 4m/(i + l)w ] 



+^[m(7 - 3)//i i+ i - 2(7 - - 4m/(z - v^ lt 

j=0 



+ [ m(Z - 3)/yU i+ i - 2(Z - l)/fi;z/ i+ i 



-(/ - l) 2 Zre£ i+ i + 4(1 - 1)k(^+i + Ai+i] u,_i,o 



-(/ - l)kwn >i+ i 



i-l 



+^[4k(^ +1 + Aj-+i) - (Z - l)iK%+i] Wi-i,^}^ + 1)Q, 



i=o 
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i-l 



+ | mn{l + 1 + 2i) hi yi - 2mK^^ + i Ztyi-j-i 

j=Q 



+ 



mn(l + 1 + 2i) + 2u [(l + 1 + 2i)a t + k] 



i-l 



+J2 [ -2mnv j+1 + 2u j+1 [(l + 2i - 2j - l)aj + k) 



3=0 



+H j+1 (2a l + b t ) ui,i-j-i}l(l + 1) 
-{ (Z + 1)[(Z + 2)(Z + 2i + 1)k + 2m(2Z + 2i + 3)cj ] 



+ 1)[ m(Z + 4)/i j+1 - 2(Z + 2)«i/ i+ i 

i=o 

+2m(2Z + 2z — 2j + l)uj+x] -wz+i^-j-i 



-(Z + l)(Z + 2)«Wi + i,. 



i-l 

+X^[ 4k ( z/ j'+ 1 + Vh) - (Z + 1)(Z + 2)K,Ej + i] wi+i t i-j-i }lQ 

j=0 



+6(z - 1){ 2(2/ + 2i + 3)a;o«i + 2,i-i 



i-2 



+]T[ 2(2/ + 2i - 2j + l)uj j+l 

3=0 

+(Z + 4)/i i+ i] M«+2,i-i- 2 } 



x/(/ + i)(/ + 3)g i+1 g ;+2 



r \ Z+l+2i 



When we substitute the series (5551)- (Bgg) and (M )-CT) into Eq. (^43|) the 
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resulting series expansion about r = R is 

oo ( 

= E (? + 1) (? + 2) /^ +2 + (g + 1) (z> + A - 2g) k, q+ i (421) 

5=0 I 

+ [ ( q + l)( q -2)-(l 2 + l-2)E 



-(q + 2)(£ + A ) + g(vi + Ai)] h l>q 



q-2 



-J2l(l 2 + l~ 2)E j+1 + (q-j + l)(u j+1 + X j+1 ) 
3=0 



(q-j- i)(^-+2 + \+2)]hi. 



g-i-i 



-Q(q - 1)[(/ 2 + I - 2)E q + 2(u q + A,)] h lfi 



When we substitute the series (p6|)-(^4l|) and (gB?D -(gIED into Eq. Q344D the 
resulting series expansion about r = R is 



= Ei - ? - !) «Wi + (9 - 1)(tto + *>o + A + q + 1) (422) 

5=0 { 



r \ q 



"Efe+i + ^+1 + Aj +1 )wi, ff _j_i - /(/ + 1) fJj^j (l - -^J 



When we substitute the series (g3|)-(g4|) and (M )-CT) into Eq. Q345|) the 
resulting series expansion about r = i? is 

oo ( 5—1 

= ^< /(/ + l)Khi, q + W + 1)« - 2mcj ] £t i;(? - 2m^tDj+i ■uj.g-j-i 

5=0 I j=0 



5-1 5 

E( 2 % + fij) ™l-l,q-j ~ 2 (7 - 1)E^J ^-1,5-i 
3=0 3=0 



-IQl+l 



5-1 5 

E( 2 % + Ai) tfj+i,9-j + 2 ( z + 2 )E c2} i 

3=0 3=0 



(423) 
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When we substitute the series (p6l)-(TO and fl407D -(M) into Eq. d346|) the 
resulting series expansion about r = R is 



= E { 2(g + l)cD ^_ 2>g+1 



(424) 



i=0 



q-l 

+ Et 2 (<? ~ + 2(1 -q + j - 1)% + (I - 3)/ij] «/-2,g-i 



+[2(Z - + (Z - 3)/2 9 ] w,_2,o}(/ - 2)Z(Z + l)Q i _iQ I 
-{ -(g + !)/[(/ - 1)k - 2mw ] Vt-i,q+i 



+Z[ (Z — l)gK + 2m(Z — q — l)u> + 2mquji 



+m(Z - 3)/i - 2(Z - l)«z5 ] vi-i >q 



g-2 



2m ( l ~ <? + + 2m (<? - J - 1 Vj+2 

3=0 

+m(l - 3)p,j+i - 2(1 - l)Ki> j+1 ] vi-i, q -j~i 



+0(q - l)Z[2m(Z - l)u q + m(l - 3)/i, - 2(Z - 1)kz? 9 ] v,_ 1)0 



9-1 ~ i 

+^[4/t(z> i + Aj) - (Z - l)lKEj] wi- hq -j\(l + 1)Q, 

j=0 
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+ { -mn(q + 1) hi, q +i + mn{q - 2z/ ) \ q 



i- 1 

-^2mKi> j+1 hi, q -j-i - (? + l){mn + 2a^ ) 

j=0 



+[mqK + 2(qa t + fy)u) - 2qa l u 1 + (2a/ + fy)// - 2mKU ] u l;l 

q-2 

+E [ 2 [(9 - J - + ~ 2 (<? - J - 1 )^j+2 

j=0 

+ (2a ; + - 2mK,u j+1 

+&(q - l)[26,Wg + (2a, + - 2m/tz> g ] u,, }Z(Z + 1) 



-{ -(I + l)(g + 1)[(Z + 2)k + 2mw ] 
+(Z + 1)[ (Z + 2)ok + 2m(Z + g + 2)£u 

-2mquo l + m(Z + 4)/2 - 2(Z + 2)kvq] 

9-2 

+^(Z + 1)[ 2m(Z + g - j + l)d) j+ i - 2m(q-j - l)Q j+2 

3=0 

+m(l + 4)^+1 - 2(Z + 2)Ki> j+ i] 
+0(g - 1)(Z + l)[2m(Z + 2)a), + m(Z + 4)/2, - 2(Z + 2)/ci/ g ] u,+i i0 

+E[ 4 «(^ + A i) - (i + + 2)«^-] «) I+ i i5 _ j }ZQ, +1 

j=0 
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+ { -2(9+1) 



<^0 ^/+2,g+l 



g-1 

+J2l 2 ( l + Q-J + 2)w 3 - - 2(q - j)Q j+1 + (I + 4)/^] uj+^-j 



+ [2(Z + 2)uj q + (/ + 4)/2,] m/ +2 ,o}/(/ + !)(/ + 3)Q l+1 Q l+2 



1 



r \ l+l+2i 



Eqs. (|411j) -( [i2"4T) make up the algebraic system ( [3471) determining the eigenvalues 
of the axial- and polar-led hybrid modes. As in the Newtonian case, one truncates the 
angular and radial series expansions at indices Z max , i max and /c max and constructs the 
matrix A by keeping the appropriate number of equations for the number of unknown 
coefficients h^i, ui t i, Wij, v^i, hi^, ui } k, w^k and Vi^. Just as in the Newtonian case, 
however, one must be aware of a certain linear dependence in the expansions about 
r = 0. With the set of relativistic equations we have chosen to work with, this linear 
dependence arises only for the axial-led hybrids and may be seen as follows. 

For a given choice of q G [0, 1,2,.. .] the set of equations, 



(I 



for all I 



( EEISD for all i 
(POD for all i 



m + 1, m + 3, . . . , m + 2q + 1 with i = 0; 
0,1, ... ,q with I = m + 2q — 2i; and 
0,1,..., q with / = m + 2q — 2i 



can be shown to be linearly dependent for arbitrary k and for any equilibrium stellar 



model. For example, taking the simplest case of q = 0, one finds that Eq. ( 418 ) with 
I — m + 1 and i = becomes, 



= w m+lfi - (m + l)i> m+ i,o, 



(425) 



Eq. ( [41 9| ) with Z = m and i = becomes, 



= m{m + l)«/i m)0 + m [(m + 1)k — 2co> ] w m ,o 



(426) 



-2mu Q m+1 [w m+ i fi + (m + 2)t> m+ i,o] , 
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Eq. ( [42 0| ) with I = m and i = becomes 



= m(m + 1)< m(m + l)nh m fi + m [(m + l)/t — 2cuo] u m $ 



(427) 



—mQ 



m+l 



(m + 1) [(m + l)(m + 2)« + 2m{2m + 3)o>q] Vm+1,0 



(m + l)(m + 2)Kw m+ i i0 



and it is not difficult to show that 



= (425) - m{m + 1){0 + Q m+1 [(m + 2)« + 2mcj ] (p27p}, (428) 



which is the claimed linear dependence. 

As in the Newtonian case, this problem can be taken care of by eliminating one 
of the linearly dependent equations for each q. To properly construct the algebraic 
system (|347| ) for the axial- led hybrid modes we use Eqs. (|411|) - (|424|) for all i except 
Eq. ([42 0|) with % = 0, for all allowed I. 
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